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PREFACE. 



The preparation of this volume was undertaken with the 
view of presenting 'in as simple and concise a manner as 
possible the fundamental principles of the Method of Least 
Squares. While it is believed that everything essential to 
the solution of all ordinary problems has been included, no 
attempt has been made to develop at length those special 
methods and forms that are so useful and almost necessary in 
case large numbers of observations of certain kinds, such, 
for instance, as those met with in geodetic and astronomical 
measurements, are to be adjusted. 

Frequent references throughout the text, and more particu- 
larly the list of works given on page v of the Appendix, will, 
however, enable the student to extend his studies in what- 
ever special direction his profession may require ; it being 
expected that this book will in such cases be looked upon 
merely as an introductory treatise. All of the works men- 
tioned have been freely consulted in the preparation of these 
pages, and the author desires in particular to acknowledge 
his indebtedness for many of the examples. 

"Dasx p. Baetlett. 
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THE METHOD OF LEAST SQUARES. 



CHAPTER I. 

GENERAL PRINCIPLES. 

1. In scientific investigations of all kinds it is frequently 
necessary to determine the values of certain quantities by 
means of actual measurements either with or without the aid 
of instruments. The observations may be made directly upon 
the values of the unknown quantities or upon certain functions 
of the unknowns. In the latter case the values of the required 
quantities must be obtained by computation from the observed 
values of the functions. In order to obtain more accurate 
values of the unknowns than would be given by a single 
measurement, or set of measurements, the observations are 
usually repeated either in the same way and under the same 
conditions or in a variety of different ways and under vary- 
ing conditions. 

Under these circumstances it will invariably be found that 
the different measurements give discordant results, the amount 
of the discrepancies varying with the character of the observa- 
tions ; and the question that now presents itself is how to 
determine from these discordant observations the true values 
of the required quantities. From the nature of the case, 
however, we can not expect to obtain our values with absolute 
accuracy; all that we can hope for is to obtain those values 
which are rendered most probable after all the observations 
are taken into account, and, further, to determine the degree 
of confidepce that can be placed in those values, 
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2. The attainment of the above results constitutes the 
primary object of the Method of Least Squares. The method 
is also employed in comparing the relative worth of different 
measurements of the same quantity, and in determining the 
equation of a curve which shall suitably represent the relation 
between two variables in cases where the exact law connecting 
them is not known. 

Also, before making any observations, we may employ the 
method to determine how precise the component measurements 
of a series must be in order to yield a required degree of pre- 
cision in the final result; or, conversely, to determine what 
the precision of the final result will be, knowing the precision 
attainable in the component measurements. This latter appli- 
cation of the method will be treated at length in the course 
on " The Precision of Measurements." 

3. Errors. The cause of the discrepancies between the 
results of our different observations is that every observation 
that is a measure is subject to error. These errors are of two 
kinds, — Constant or Systematic Errors and Accidental Errors. 

4. Constant Errors are errors which in all measures of the 
same quantity, made with the same care and under the same 
conditions, have the same magnitude, or whose presence and 
magnitude are due to some fixed cause. These constant 
errors may be of several classes, which are designated as 
follows: — 

First. Theoretical Errors, such as those due to the refrac- 
tion or aberration of light, the effect of a definite change in 
temperature or moisture on our standards of measurement, 
etc. As soon as their causes are known the magnitude of 
these errors may be calculated and their effect eliminated 
from the observations. 

Second. Instrumental Errors, such as errors of division of 
graduated scales, defects in micrometer screws, eccentricity of 
circles, etc. These errors will be discovered by an examina- 
tion of the instruments and their effects eliminated from the 
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observations, either by a particular method of using the 
instruments or by subsequent computation. 

Third. Personal Errors. These are due to personal pecu- 
liarities of an observer, who always answers a signal too soon 
or too late, always estimates a quantity smaller than it is, etc. 
The character and magnitude of these errors may be deter- 
mined by a study of the observer, his "Personal Equation" 
may be obtained, and his observations thus corrected for this 
source of error. 

5. Mistakes. Although of a somewhat different character, 
these should be considered in connection with constant errors. 
A mistake is made when a figure 3 is read for a figure 8, or 
when in reading a graduated circle which is numbered in both 
directions the angle is read 48° instead of the complementary 
angle 47°, etc. These mistakes are usually of such a charac- 
ter that they may be detected by an inspection of the observa- 
tions and a proper correction made. 

6. Accidental Errors are errors due to irregular causes, 
whose effect upon the observations is not determined by any 
circumstances peculiar to that particular set of measurements, 
and which cannot therefore be computed and allowed for 
beforehand. Such errors are those due to sudden changes in 
refraction owing to sudden and unobserved changes in tem- 
perature; unequal expansion of different parts of an instru- 
ment with change in temperature ; shaking of an instrument 
in the wind, etc. But most important of all are those errors 
which arise from imperfections in the sight, hearing, and other 
senses of the observer, which render it impossible for him to 
adjust and use his instruments with absolute accuracy. 

After a full investigation of the constant errors, the observer 
should diminish the accidental errors as much as possible, both 
in number and magnitude, by taking every precaution and care 
in the measurements themselves. The problem now remains 
to combine the observations so that the remaining accidental 
errors shall have the least probable effect upon the results, 
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and it is to bring about tbis combination of observations that 
we employ the Method of Least Squares. 

When no more observations are made than are sufficient to 
determine one value for each of the unknown quantities, we 
must accept these values as the most probable ones. But if 
additional observations are made leading to discordant results, 
we can not take any one of them as the correct value, and in 
fact, as already stated, we shall probably not be able to obtain 
the true values of the unknowns. All that we can do is to find 
values of the unknowns which shall remove the discrepancies 
between the different observations and which shall be those 
values that are rendered most probable by the existence of 
the observations themselves. 

On first thoughts it may seem that these accidental errors, 
being due to so many different and unknown causes, will be 
beyond the scope of mathematical investigation. Neverthe- 
less, the theory of probability requires that these errors shall 
follow in magnitude and frequency a law that is capable of 
exact mathematical expression, and experience confirms the 
correctness of this law. 

For more extended remarks on these subjects see — 

Holman, " Discussion of the Precision of Measurements," pp. 1-14. 
Merriman, " Text-Book of Least Squares," pp. 1-6. 
Chauvenet, " Spherical and Practical Astronomy," pp. 469-473. 
Wright, " Treatise on the Adjustment of Observations," pp. 11-18. 



LAWS OF ERRORS OF OBSERVATION. 

7. The derivation of the general laws of the occurrence of 
errors of observation, and of the processes for determining 
the most probable values of the unknown quantities, will be 
based upon the following 

Axiom. If a series of n direct observations. Mi, M^, . . . M^, 
is made upon the value of a quantity 31, all the observations 
being made with the same care and under the same circum- 
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Stances, the most probable value M^ of that quantity is the 
arithmetical mean of the observations. Or 

" n n 

^ 8. The Real Error (a;) of an observation is the difference 

between the observed value of the measured quantity and the 

real value. 

~~ 9. The Residual (w) of an observation is the difference 

between the observed value of the measured quantity and 

the value rendered most probable by the existence of the 

observations. 

10. Example. Eight observations are made upon the 
resistance of a coil of wire, the true resistance being 512. 
Find from these observations the most probable resistance, 
and also the real errors and residuals. 



Observations. 


Real Errors. 


Residuals. 


M. 


X 


V 


512.4 


+ •4 


+ .30 


512.2 


+ •2 


+ .10 


511.9 


-.1 


-.20 


512.3 


+ .8 


+ .20 


511.8 


-.2 


-.30 


512.3 


+ .3 


+ .20 


511.9 


-.1 


-.20 


512.0 


.0 


-.10 


Mean = 512.10 




%v= .00 



From the observations, then, we should say that the most 
probable resistance of the coil is 512.10. It will also be 
noticed that the sum of the residuals is zero. That this is a 
general result following from the assumption of the arith- 
metical mean as the most probable value may be proved as 
follows: If the observations are il/i, 31^, . . . Mn, the 
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arithmetical mean 3fo, and the residuals Vi, v^, . . . v^, then 
we have 

Vi = Mi-Mo, V2 = M2 — Mo, ...v„ = Mn — Mo. 

= 'SiM—'2,M, since ^= 

.-. Sv^O (2) 

11. Weighted Observations. The weight of an observa- 
tion expresses its relative worth compai-ed with other obser- 
vations. Thus, if six observations are made upon the value 
of a quantity, five of which give the same result, while the 
sixth differs, in combining these two different results to 
obtain the most probable value of the unknown, the first 
value ought to have five times the influence upon the final 
result that the second has, since it has taken five times as 
much labor and time to obtain it. Hence in general we may 
say — 

12. The Weight (p) of an observation may be considered 
as representing the number of times the observation has been 
repeated and the same result obtained. 

The weights assigned to observations may be due to a 
variety of causes, as difference in skill of observers, difference 
in the instruments used or the circumstances under which 
the observations are made, etc. But whatever the cause, the 
effect on the final values of weighting an observation will be 
the same as indicated in the preceding paragraph. 

13. Example. Suppose n observations, M-^, M^, . . . M^, 
of weights pi, jt>2, . . . p„, are made upon the value of a 
quantity M. To find the most probable value Mi, of the 
quantity. 

From the above interpretation of the meaning of weight, 
we may consider that the whole number of observations is 
Pi-\-J''2-\- ■ ■ • Pm or %p, and that the result M^ has been 
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obtained in p^ observations, M^ in p^ observations, etc., 
Therefore, by (1) 

^'= s^ = -W ^^ 

Ma is called the General Mean. 

If the residuals are Uj, v^, . . . u„, we have 

v-^=M^ — Mo, V2 = Mi — Mo, . . . v„= Ma — M» 

Pil^l +P2V2 + • • . Pn'On = PxMi -\-PiM^ + • . • PnMi — M>^P 

^pv=0~ (4) 

Which shows that in the case of direct observations of differ- 
ent weights the sum of the weighted residuals is zero. 

14. If the observations are not made directly upon the 
values of the required quantities, the method of adjusting 
the results so as to obtain the best possible values of the 
unknowns will depend upon the laws which govern the dis- 
tribution of the errors of these observations. It is found in 
practice that the accidental errors of observations follow cer- 
tain well defined laws, and what these are may best be seen 
by taking an actual example. 

15. Example. One thousand shots are fired at a target 
which is divided into a number of horizontal sections by lines 
one foot apart, the centre line of the target being in the 
middle of one of these spaces. The shots were distributed 
as follows : — 



7n Space. SJtota. 


In Space. 


Shots. 


In Space. 


Shots. 


_|_5^to+4^ 1 


+Uto+i 


190 


-2^ to -3^ 


79 


+4i«+3i 4 


+ i " - * 


212 


-3* « -4* 


16 


+ 3^ « -\-2i 10 


- i " -H 


204 


_4^ « _5^ 


2 


4-2* " +1* 89 


-n « -2i 


193 







In this case the errors are evidently the distances of the 
shots from the centre of the target. Further, as far as can 
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be judged from these one thousand shots, if another shot is 
fired the probability that this shot will fall between the 
lines 



4- 5^ and -1-4^18 


.001 


— i and —1^ is 


.204 


+^ " +3i" 


.004 


-1^ " —2^ " 


.193 


+ H " +2i« 


.010 


-2^ " -3^ « 


.079 


+ 2J « +1^«« 


.089 


_3i « -4^ " 


.016 


+ U " + i" 


.190 


-^ " -5J « 


.002 


+ i " - i " 


.212 







The sum of the above probabilities is unity, and, therefore, 
as far as the preceding shots show the lOOlst shot will cer- 
tainly hit the target. 

16. Now using as abscissas the distances of the horizontal 
lines from the centre of the target, and as ordinates the num- 
ber of shots falling in the corresponding spaces, we may 
construct the following figure : — 



Figure 1. 



And if the entire area of this figure is taken as unity then 
the area of each rectangle will denote the probability of a 
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shot, if fired, falling within the corresponding space of the 
target. 

The graphical representation of the accidental errors of 
observation will always give a figure similar to the above. 
Hence denoting errors by abscissas, and their frequency by 
ordinates, the law of error of any series of observations may 
be represented by a curve whose general form is determined 
by Figure 1. This curve is called the " Curve of Error," and 
is shown in Figure 2. 

r 



Figure 2. 




PDU 



In order that this curve may represent exactly the distri- 
bution of the errors in any given series of observations it 
ought to meet the axis of X at some definite distance to the 
right and left of the origin and coincide with the axis from 
there on, for in all actual observations there is a limit beyond 
which no errors occur. But as the exact point of meeting 
could not be determined for any given case, and as it would 
not be possible to obtain the equation of such a curve, we 
make it asymptotic to the axis of JC, taking care that the 
error thus introduced shall in any set of observations be so 
small as to be negligible. 



10 
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17. An inspection of Figures 1 and 2 will now exhibit 
some of the general laws of errors of observation and the 
corresponding properties of the curve of error. 



Laws of Error derived from an 
inspection of Figure 1. 



Representation of these laws by the 
Curve of Error. 



First. Small errors are more 
frequent than large ones. 

Second. Positive and negative 
errors of the same absolute mag- 
nitude are equally likely to occur. 

Third. The probability of the 
occurrence of very large errors is 
very small. 

Fourth. The frequency of any 
error depends upon the magnitude 
of that error. 



The maximum point of the 
curve is on the axis of Y. 

The curve is symmetrical with 
respect to the axis of Y. 

The curve is asymptotic to the 
axis of X. 

The equation of the curve will 
be of the form 



2/=<^(a:) 



(5) 



18. If, now, the total area between the curve and the axis 
of X be denoted by unity the probability that the error of 
any given observation will fall between the magnitudes x and 
x-\- dx will be represented by the area included between the 
curve, the axis of J^ and the ordinates of the curve at the 
errors x and x -|- dx, or by 



y dx = <j>(x) dx 



(6) 



And this probability will be known as soon as we find the 
form of the function <l>{x). 

19. The above expressions in (5) and (6) are the ones 
that we should use if we regard the curve of error as repre- 
senting the law of occurrence of errors of observation. If, 
however, we look upon the curve as expressing the law to 
which we must make the residuals conform, in order that the 
values of the unknown quantities obtained from them may be 
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the most probable values, we should replace x hy v and use 
the expressions 

y='i>{v) (7) 

and ydv=<t>{v)dv (8) 

for (5) and (6), respectively. 

THE EQUATION OF THE CURVE OF ERROR. 

20. Let n observations, all of the same weight, with 
results Ml, Ji^, . . . 31^, be made upon any function or 
functions of a number of unknown quantities Sj, Zn, . . . z^; 
and let the residuals of M^^, M^, ■ • ■ M„ be v^, v^, . . . v„, 
and the probability of the occurrence of these residuals be 
(f>(Vi) dv, <j)(V2) dv, . . . 4i(Vn) dv, respectively. Then the 
probability of the simultaneous occurrence of all these 
residuals will be 

P=: <^(«l) ^V,) . rr<l>(v„) {dv)n (9) 

.-. log F = log <^(t)i) -f . . . log <^(w„) + « log dv (10) 

Each different method that might be adopted for computing 
the values of the unknowns Sj, z^, . . . Zg would lead to a dif- 
ferent set of residuals Vi, v^, . . . v„; but obviously that set 
of values of z^, z^, . . . Zq should be considered the best which 
corresponds to the particular set of residuals Uj, v^^ ■ • • '"ni 
the probability of whose occurrence is greater than that of 
any other set. 

Therefore the most probable values of gj, Sj, . . . Zq are 
those that make P in (9), or log P in (10), a maximum. 

The values of z^, Sj, . . . Zq corresponding with this latter 
condition are those that satisfy equations (11). It may be 
noticed that these equations also express the preliminary 
conditions leading to a minimum value of logP, but the 
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nature of the problem is such that a maximum value of P 
evidently exists while a minimum does not, and it is there- 
fore unnecessary to investigate further the mathematical 
conditions for a maximum. Hence we have 

SlogP 1 d<l>(Vi) , 1 d «^(t>„) ^ 

I • • • Ai-j, ^ ;5«- 



(11) 



3z2 ~~<^(Vi) azj "'"■■"c^(y„) Sgj 



aiogp ^ 1 5.^(t>i) _^i*K)=o 



and for convenience we may put 
substituting in (11) we have 



^{v) (13) 






(14) 



^(M 35^4-^(^5-^+. ..^K)5^; = o 

These equations contain all the unknowns z^, z^, . . . Zq, and 
there are as many equations as unknowns, hence as soon as 
we find the form of the function i/'(«) we can solve these 
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equations for the most probable values of Zi, Zo, . . . Zg. Since 
we have considered the general case, and the above results 
are to hold true whatever the number of unknown quantities 
and the form of the functions observed, we may deduce the 
form of ip{v) by solving a special example. 

Example. Let n observations of equal weight be made 
upon the value of a single unknown Zi, with results Mx, M^, 
. . . 3Tn, and let the residuals be Vi, v^,... v„. Then the 
most probable value of Zj is given by 

Zl = 3fx — Viz= 3f2 — Vi=: . . . M„ — Vn 

differentiating with respect to Zi, 

dzi dzi ' ' ' dZi *■ ' 

substituting (a) in (14), changing all the signs, we have 

^(«0 + "ACM +■••</'("«) = (b) 

But in this case, as was shown in (2), 

Vx-\- V2-\-. . .Vn=0 (C) 

In order that (b) and (c) may both be true the functional 

symbol xj/ must indicate multiplication by a constant. That 

is, in general 

xl/[v)=zcv (15) 

Substituting this in (13) and (12), 

dd>(v) , dv 

1 a<^(y) dv 

^^^^^^^^ ^(^-ai- = "''5F 

Integrating, log <^( y ) = icv" -\- c' 

^(w) = eJ<»" + «' 
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Since yz=,j>(^v) is the equation of the curve of error, (7), 
we may therefore write it 



y 



= ke^<^'' 



But on examination of the curve, y is seen to be a decreas- 
ing function of v, and hence the exponent of e is essentially 
negative. Accordingly we will write our equation in the 
form 

yz=1iie-f>^^ (16) 

the values of k and h depending upon the character of the 
observations, but in all observations of the same kind and 
weight having the same values. 

This equation represents the law in accordance with which 
the residuals must be distributed in order that the best results 
may be obtained from our observations. But, as before 
mentioned, if we wish our curve to represent the most prob- 
able distribution of the real errors of observation we should 
write the equation in the form 

y=ke-^'^' (17) 

Hereafter we shall use without further remark either form 
of the equation according to the aspect in which we are 
considering our curve. 

An inspection of the above equation will show that it 
satisfies all the conditions noted in discussing the form of the 
curve of error in paragraph 17. 

21. It is important to notice that in all discussions in the 
Method of Least Squares the number of observations is sup- 
posed to be large and always greater than the number of 
unknown quantities. As will be illustrated later on, para- 
graph 91, whenever this is the case there is a remarkable 
agreement between the results obtained in practice and those 
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indicated by the theory. And even when the observations 
are few in number the method still affords the best means at 
our command for their adjustment, the results obtained 
merely having a smaller weight than they would have had if 
derived from a greater number of observations. 

THE METHOD OF LEAST SQUARES. 

22. We are now in a position to see whence comes the 
name " Least Squares." 

In paragraph 20 it was pointed out that whenever we 
make a series of observations, each observation of the set 
having the same weight, the most probable system of values 
of the unknown quantities will be that which corresponds 
with the set of residuals the probability of whose occurrence 
is a maximum. That is, the best set of values of the 
unknowns will be that which gives a maximum value to 

But from equation (16) this reduces to 

P = Fe-'^'C"!' + ^2^ + • • • ''«') (dv)" (18) 

Since the exponent of e in this expression is negative, 
evidently P will be a maximum when 

Vi'-]-Vi^-{- ■ ■ ■ Vn^^'Zv^ is a minimum. (19) 

Hence the adjustment of observations by the Method of 
Least Squares is based upon the principle that the most prob- 
able system of values of the unknowns is that which renders 
the sum of the squares of the residuals a minimum. Hence 

the name. 

The conditions for a maximum value of P were expressed 
in equations (14), and since it has been shown that the 
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function i^ means multiplication by a constant, those equa- 
tions reduce to the following, called 

NORMAL EQUATIONS. 






(20) 



3fl , SVi , BVn 

An inspection of these equations will show that they also 
express the conditions that will make the sum of the squares 
of the residuals a minimum. 

In the adjustment of observations the above are the funda- 
mental equations. In order to obtain the most probable 
values of the unknowns in any set of observations, all that is 
necessary is to form the Normal Equations for that set and 
solve them simultaneously. The examples already solved for 
direct observations are merely special cases of the above 
general solution. 



CHAPTER n. 

THE ADJUSTMENT OF OBSERVATIONS. 

INDIRECT OBSERVATIONS. 

23. In the determination of the values of quantities by 
means of observations the functions of the unknowns that it 
is necessary to observe may be of any form, but if they are 
not linear the normal equations derived from them are likely 
to be complicated and difficult, if not impossible, to solve. 
Hence if the observations are not upon linear functions of 
the unknowns, the first step will be to reduce them to 
equivalent linear expressions by transformations depending 
on the character of the observed functions; see paragraphs 
43 and 44. It will be necessary to consider, therefore, the 
method of adjusting observations on linear functions alone, 
and the procedure in cases of this kind may be illustrated by 
the following simple example. 

24. Example. -S'l, S^, S^, are three solids whose masses 
are required. Not having standard weights enough to obtain 
all these masses directly, by varying the distribution of the 
solids in the pans of the balance the following observations 
are made : — 

xS'i := <5^2 -)- 1.7 grams. 
Ss = 2.4 « 

If, now, the most probable values of Si, S^, S^, are repre- 
sented by 2i, 22, S3, and the corresponding residuals of the 
observations by Vi, V2, . ■ ■ v^, the values of these residuals in 
terms of Sj, Zj, Zg may be found from the above observations 



18 METHOD OF LEAST SQUARES. 

by transposing all terms to the first members of the equa- 
tions, and we obtain at once the following, called 

OBSERVATION EQUATIONS. 

Zl — Z2 — 1.7 := Wi 

Z3 — 2.4 = V2 

— Zl -f- Z2 + Z3 — 1.0 = W3 

Z2 — Z3 — 3.0 ^ ^4 



(A) 



Applying equations (20) as formulas we have the 
NORMAL EQUATIONS. 

(2l_22-1.7)+(-Zi+Z2+Z3-1.0)(-l) = 

(Zl - Z2 - 1.7)(- 1) + (- Z1 + 32+ 23 - 1-0) 

+ (Z2-Z3-3.0) = (B) 

(Z3-2.4)+(-Zi + Z2+Z3-1.0) 

-f (S2-33-3.0)(-l) = 
Simplified, these become 







2zi - 222 - 


S3- 


-0.7=0 


(a) 






— 2zi-|-3z2 




- 2.3 = 


0>) 






- Zl + 


3Z3- 


-0.4=0 


(c) 


3X(a) 




6zi — 6z2 — 


3Z3 


-2.1 = 




(c) 


- 


+ 


3Z3 


-0.4=0 






5zi — 6Z2 




— 2.5=0 




2X(b) 


- 


-4zi + 6z2 




— 4.6=0 






Zl 




— 7.1 = 












2i=7.1 


(d) 


substitute 


(d) 


in(b) 




Zj = 5.5 


(e) 


(( 


(d) 


in (c) 




Z3=2.5 


(f) 



An inspection of the work in this example will show that 
for the adjustment of observations of equal weight on linear 
functions of the unknowns we may derive the following : — 
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25. Sale. J^or each observation write an " Observation 
Equation " / then for each unknown form a " Normal 
Equation,'''' by r)\idtiplying the first member of each observa- 
tion equation by the coefficient of that unknown in that 
equation, adding the results and placing the sum equal to 
zero. Solve these equations simultaneously for the values 
of the unknoions. 

In solving for the most probable values of the unknowns 
the second members of the observation equations are very 
commonly written zero instead of v^, v„, . . . Vn- For this is 
the form in which the equations naturally appear, and if the 
observations were exact the residuals would actually all be 
zero. The method of solution is the same in either case. 

26. Observations of Unequal Weight. If the observations 
are not all of equal weight the same method will apply, 
except that in the formation of the normal equations each 
observation equation will be used the number of times denoted 
by its weight. Thus in the last example if the observations 
have the weights 4, 9, 1, 4, the normal equations will have 
the same form as in (B), page 18, but each part of each 
equation will be multi-plied by the weight of the observation 
equation from which it is derived. This will give the 

NORMAL EQUATIONS. 
4(3i-S2-1.7) + (-Zi + Z2 + Z3-l-0)(-l) = 
4(zi - s„_ 1.7)(- 1) -I- (- Zi + z, + 2, _ 1.0) 

-f 4(Za- 3,-3.0) = 
9(23 - 2.4) + (_ Si + z, + Z3- 1.0) 

-|-4(s2-Z3-3.0)(-l) = 
or reduced 

5zi -^ 5s2 — Z3— 5.8 :=0 (a) 

— 5zi + 9z2— 3z8— 6.2=0 (b) 

_ Sj_3zj4-14z3— 10.6=0 (c) 

the solution of which gives 

Zi=7.07 Zj=5.42 83=2.42 
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Further it will at once be seen that if />i, p^, ■ ■ ■ Pw are the 
weights of the corresponding observations, equations (20) 
take the general form: — 

WEIGHTED NORMAL EQUATIONS. 

^^'^^a^+^^^=a;.^+---^"*'"s^ = * (21) 



Hence for the formation of the normal equations in weighted 
observations on linear functions of the unknowns, we have 
the following: — 

27. Eule. I^or each observation write an " Observation 
Equation. " / then for each unknown form a '■'■Normal 
Equation,'''' hy multiplying the first member of each observa- 
tion equation by the coefficient of that unknown in that 
equation and by the weight of that equation, adding the 
results and placing the sum equal to zero. Solve these 
equations simultaneously. 

28. The same result will be obtained if we begin by 
multiplying each observation equation by the square root of 
its weight and then proceed according to the first rule 
(paragraph 25). 

This result illustrates the important principle that m,ulti- 
2^^yi»g « set of equations by the square roots of their weights 
reduces them all to equivalent equations of vieight unity. 

29. Helation between the Weight of an Observation and 
the Value of it. If in paragraph 20 the n observations 
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have weights 2hi Pii • ■ ■ Pm and the quantity h values 
Ai, Aj, . . . A„, then equation (18) becomes 

= kik^ . . . A;„e-(''i''''i' + ^^^''^^ +■■■ ^'«"') (dv)" (22) 

The most probable set of values of the unknowns is that 
which makes jP a maximum, and P is a maximum when 

hi^ Vi^ -\- h^^ v^^ -\- . . . TinVn is a minimum. (23) 

The conditions for a minimum value of this expression are 
the following, which are then for this case the 

NORMAL EQUATIONS. 



s, -r ■ ■ ■ nn Vn Q^^ 



^^"'^^di+^^"'^^di-^---^-'''-di=^ 



(24) 



But equations (21) are also the normal equations for 
this case. Hence (21) and (24) must be identical, and 

Pi : P2 ■ ■ ■ Pn^ hi' : W : ■ . . h„^ (25) 

That is, the square of h is proportional to the weight of 
the observation. Accordingly, since h increases in value as 
the quality of the observations is improved, it is called "The 
Measure of Precision." 
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Further, it follows from (23) and (25) that the most 
prohahle system of values of the unknowns will be that in 
which 

PiVi -{- x>2V2^ -\- . . .pn^n is a minimum. (26) 

And this is the most general form of statement of the 
principle of " Least Squares." The same principle is repre- 
sented in equations (21). 

30. Computation of Corrections. If large numbers occur 
in the observations it is better to compute the most probable 
corrections to apply to the observed values rather than the 
most probable values of the unknowns themselves. In this 
way we can often avoid a large amount of numerical work. 

31. £Jxample. Pi, P^i -f ai P*-, -Psi are five points whose 
altitudes above the mean level of the sea are to be determined 
from the following observations of difference of level. 

Pi = 573.08 P4-P2= 170.28 

P2-A= 2.60 ^4 — ^5 = 425.00 

A = 575.27 P5 = 319.91 

Pa -^2 = 167.33 P5 = 319.75 
Pi-Ps— 3.80 

An inspection of these observations shows that we may put 

Pi=:573 + Zi ^4=745+34 

P, = 575 + 22 P5 = 320 + 35 (A) 



Ps = 742 + 



Sb 



where Zj, Zj, Zj, z^, Zj, are small corrections whose most 
probable values are to be determined. We now have for 

OBSERVATION EQUATIONS. 
573 + zi- 573.08=0 or Zi-.08 = 

675 + Z2-573-Z1- 2.60=0 or Zj-Zj-. 60 = 
575 + Zj- 575.27=0 or z^— .27 = 



or 


2, _ gj _ .33 = 


or 


^4 - Zs — -80 = 


or 


Z4 - 2j - .28 = 


or 


24 — 25 =0 


or 


25+. 09=0 


or 


25 + .25=0 
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742-|-2s_575-Zj- 167.33 = 
745 + Z4 — 742 - 28 — 3.80 = 
745 + 24 — 575 - 2, - 170.28 = 
745 -j- 24 — 320 - 25 - 425.00 = 

32O4-Z5- 319.91 = 

320 + 35-319.75 = 

From these we now form the 

NORMAL EQUATIONS. 

22i- S2 + .52=0 

- Zi + 4z2- 23- Z4 - .26=0 

— Z2 + 223- Z4 + .47=0 

— Z2 - 23 + 3Z4- 25—1.08=0 

■^ Z4 + 3f;5+ .34 = 
and solving, 

2i=-.19; Z2 = .14; Z3=.05; Z4=.43; 25 = .03 

Substituting these in equations (A) we have for the most 
probable altitudes, 

Pi = 572.81 ^3 = 742.05 ^5 = 320.03 

^2 = 575.14 P4= 745.43 

If the original observation equations had been retained, 
the independent terms in the normal equations would have 
been 

570.48 240.26 163.53 599.08 214.66 

32. Significant Figures. The adjustment by the Method 
of Least Squares of observations which occur in practice, 
although not difficult, is apt to be long and laborious. Hence 
to reduce this labor as much as possible it is of great import- 
ance that careful attention should be given in the solutions 
to the proper use of significant figures. When in doubt, 
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liowever, as to the proper number of figures to retain it is 
better to keep too many rather than too few, as the superfluous 
figures can be rejected at the end of the computation ; while 
if too few are retained the results obtained from the compu- 
tations will be worthless. 

For a general discussion of the subject of significant figures 
see Holman's " Precision of Measurements," pages 76 to 84, 
but for the present the following rules will suffice for most 
cases. 

Bule 1. Tn casting off places of figures increase by 1 the 
last figure retained, when the following figure is 5 or over. 

Rule 2. Tn the jyrecision measure retain two significant 
figures. 

Rule 3. In any quantity retain enough significant figures 
to include the place in which the second significant figure 
of its ])recision measure occurs. 

Rule 4. When several quantities are to be added or 
subtracted, ax>ply Rule 3 to the least precise and keep only 
the corresponding figures in the other quantities. 

Rule 5. When several quantities are to be multiplied or 
divided into each other, find the piercentage precision of the 
least precise. If this is 

1 pjer cent or inore, use four significant figures. 
.1 " " " " " five " " 

.01 " " " " " six " " 

in all the worh. If the final result obtained in this way 
confiicts toith Mule 3, apply the latter. 

Rule 6. Wlien logarithms are used, retain as many 
places in the mantisscs as there are significant figures 
retained in the data under Rule 5. 

The application of these rules is not always possible in the 
course of the work, since the precision measures may not be 
known until the end of the computation. But as a general 
rule it is sufficient in direct observations to retain one more 
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place of figures than is given by the individual observations, 
and in indirect observations to retain two additional places. 

CONDITIONED OBSERVATIONS. 

33. Conditioned Observations are those in which the 
unknown quantities must be determined not only so as to 
satisfy as closely as possible the observation equations, but 
also so as to satisfy exactly certain other conditions. These 
conditions must be less in number than the unknown quantities, 
otherwise the unknowns could be determined from the con- 
ditions alone. 

The adjustment of observations of this class may be reduced 
to the method already used for unconditioned observations in 
the following way. 

The observations are represented by " Observation Equa- 
tions," and the conditions by certain other equations, called 
" Condition Equations." 

Between these two sets of equations we will eliminate as 
many unknowns as there are conditions. From the resulting 
equations, which will be the same in number as the observa- 
tions, we will form in the usual manner the " Normal 
Equations " for the remaining unknowns. Having solved 
these normal equations and substituted the results in the con- 
dition equations, we shall obtain the values of the unknowns 
first eliminated. 

All conditions of the problem are now fulfilled, for the 
condition equations are satisfied exactly and, moreover, 
according to the principle of Least Squares our results 
are those rendered most probable by the existence of the 
observations. 

As in the example last considered, it is often more 
advantageous to compute corrections to the observed values 
of the unknown quantities rather than the values of the 
quantities themselves. 
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34. Example. Find the most probable values of the angles 
of a quadrilateral from the observations, 



^=101° 13' 22" weight 3 

£= 93 49 17 " 2 

C= 87 5 39 " 2 

D= 77 52 40 " 1 



(A) 



^_|_J5-)-C+Z>=360° 0' 58" 

The condition to be satisfied is in this problem 

^_|_^-|.C+Z»=360° (B) 

Let 8i, Zj, Zg, Z4 be the most probable corrections to add 
to the observed values. This gives for 

OBSERVATION EQUATIONS 



Zi ::= weight 3 

22=0 " 2 

Z3=0 " 2 

z, = " 1 



(C) 



and for the 



CONDITION EQUATION 

Zi + Z2+Z3 + Z4 + 58=0 (D) 

Eliminating z^ between (D) and (C), the equations from 
which the normal equations are to be derived become 



Zi=: weight 3 
2.= " 2 

Z3=0 " 2 (^^ 

Zi + Z2 + Zj + 58 = " 1 



THE ADJUSTMENT OF OBSERVATIONS. 27 

Applying the rule in paragraph 27, these give the 

NORMAL EQUATIONS 

4s,+ z,+ Z3 + 58 = 
Si + Sz^-j- S3-}-58 = (F) 

zi-j- 22 + 328-1-58=0 

Solving, and substituting the results in equation (D), we 
find 

zi=- 8.29 Z3=_ 12.43 „ 

S2=- 12.43 Z4=- 24.85 ^ ' 

Applying these corrections to the observations (A), the 
most probable values of the angles are 



^=101° 13' 13".71 

-B= 93 49 4.57 

(7= 87 5 26.57 

D= 11 52 15.15 



(H) 



Note. In eliminating tmknowns between the observation 
and condition equations care must be taken that the obser- 
vation equations are not combined with each other or multi- 
plied by any quantity. For if this is done the loeights of 
the observation equations will be altered. (See §28.) 

35. In the above example it is evident that the corrections 
to be applied to the different observations are inversely as 
their weights. And, in general, when there is but one 
equation of condition, the observations expressing direct 
determinations of the unknowns, the corrections will be pro- 
portional to the coefficients of the unknowns in the equation 
of condition divided by the weights of the corresponding 
observations. A proof of this is given in paragraph 117. 
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The most common case is that in which these coefficients 
are all unity, as in the example just solved, and we may then 
derive the 

Rule. Find the difference between the theoretical and 
observed results and divide this correction among the 
observations in the inverse ratio of their weights. 

In the last example the sum of the observed angles exceeds 
360° by 58". Therefore the correction to be applied to A is 

1 

- 58 X ^ ^ ^ ^ = - 58 X y = - 8.29 

8- + ^ + ^+^ 

EMPIRICAL FORMULAS AND CONSTANTS. 

36. In the work so far considered the observations are 
supposed to be made either directly upon the values of the 
unknown quantities or upon some function of the unknowns 
whose form, and the constants entering into it, are definitely 
known. But another sort of problem frequently occurs, 
in which observations are made upon the values of a certain 
variable and the corresponding values of some function of it, 
the exact form of the function not being known. The object 
in this case is the determination of the most probable form of 
the function and the values of the constants involved ; that 
is, the derivation of the algebraic expression best representing 
the law connecting the variable and function. 

This expression may be looked upon as the equation of a 
curve, abscissas denoting values of the variable and ordinates 
values of the function, and for all values of the variable 
within the range of the observations we may determine from 
it the most probable values of the function corresponding. 
But except in special cases, where the number of observations 
is large, where the law connecting variable and function is 
well defined, and where the equation obtained is an accurate 
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representation of this law, it cannot be assumed to apply- 
beyond the range of the observations. And in no case would 
it be safe to make use of the curve very far beyond the limits 
of the observations. 

37. The Method of Least Squares will not assist in deter- 
mining the form of the function. This must be settled upon 
beforehand, either from theoretical considerations or by 
constructing a plot, using values of the variable as abscissas 
and of the function as ordinates, when the smooth curve 
drawn through the points thus obtained will indicate the 
form of equation to be used. 

It is to be observed that this is a method of trial and will 
not necessarily give the most probable form of the function ; 
and in fact we may not be able to obtain the form that would 
be absolutely best. Further, several forms of equation may 
be known which would represent well the plotted points. In 
such a case that should be considered the best in which the 
sum of the squares of the residuals is found to be the least. 

38. As soon as the form of the function is decided upon it 
should be reduced to the linear form, and the determination 
of the values of the constants involved is then a simple 
application of the preceding methods. 

As the " Observation Equations " in any given problem will 
all be of the same kind, it is usually advisable to write out 
the general form of the "Normal Equations" and arrange 
the computations in tabular form, while the retention of 
the proper number of significant figures is of particular 
importance in this work. 

39. A case that frequently occurs is that in which the 
quantity y is a constantly increasing function of the vari- 
able X, or where the plotted curve is approximately parabolic 
in form. Here the equation 

y^ A + Bx+Cx'+I)x'+ ... (27) 

may be taken to represent the relation between the variable 
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and the function. The larger the number of terms taken in 
the second member, the more accurately may the equation 
obtained be made to represent the results of the observations ; 
but the labor involved increases rapidly with increase in the 
number of terms, and if the plot shows a very nearly straight 
line the first two terms alone may suffice. 

40. Example. In measuring the velocity of the current 
of a river the following results were obtained : — 



Depths. 


Velocities 


X 


V 


1 


4.86 


2 


5.14 


3 


5.15 


4 


4.85 


5 


4.24 


6 


3.36 


7 


2.16 


8 


0.67 



The velocity at the surface is 4.250. Find the equation of 
a curve which will express the relation between x and V. 
Plotting the observations we find the curve 

V 
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This is approximately parabolic in form and passes through 
the fixed point (0, 4.25). Therefore the relation between x 
and Y may be expressed by the equation 



F=4.25 + ^a;-f Ca;2 



(A) 



and substituting in this the corresponding values of x and V 
as given by the observations, we shall have eight observation 
equations from which the most probable values of £ and C 
are to be computed. 

All of the observation equations being of the form (A) we 
have the 

NORMAL EQUATIONS 

(7 2a;*+^2a;' + 4.25 2x'- 2Fa;2= (a) 

C-^x'-^-B 2a;'+ 4.25 2a; — 2 Fa; = (b) 

For computing the coefficients in these equations it is most 
convenient to arrange the following table. 



X 


F 


Fa; 


x" 


Vx^ 


x" 


X* 


1 


4.86 


4.86 


1 


4.86 


1 


1 


2 


5.14 


10.28 


4 


20.56 


8 


16 


3 


5.15 


15.45 


9 


46.35 


27 


81 


4 


4.85 


19.40 


16 


77.60 


64 


256 


5 


4.24 


21.20 


25 


106.00 


125 


625 


6 


3.36 


20.16 


36 


120.96 


216 


1296 


7 


2.16 


15.12 


49 


105.84 


343 


2401 


8 


0.67 


5.36 


64 


42.88 


512 


4096 


36 




111.83 


204 


525.05 


1296 


8772 


2a! 




2 Fa; 


2a;^ 


2Fx» 


2a;» 


2x* 
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Substituting these results in (a) and (b) we have 

8772 C+ 1296 ^ + 341.95 = (c) 

1296C+ 204^+ 41.17 = (d) 

and solving, 

C—— .1493 B — .7465 (e) 

Therefore the required equation is 

V = 4.25 + .7465a; - .1493x2 (f ) 

Whenever the quantities in the observations are so large 
that the use of logarithms is desirable, these can best be put 
in the same columns directly over the natural numbers. 

41. It is to be remarked that in solutions like the above 
we assume, from our method of forming the observation and 
normal equations, that the observations on the values of the 
function are alone subject to error, the observations on the 
variable being supposed to be exact or to have errors so small 
as to be negligible. 

42. Periodic Phenomena. If as the variable increases the 
function passes through recurring values, that is, if y is a 
periodic function of x, some form of trigonometric equation 
would be the proper one to select. For instance, a good form 
to use is 

y = A-{-B8in ^x + C cos ^x (28) 

where A, B, C are the constants whose most probable values 
are to be found, and m is the number of units of x comprised 
in the entire cycle of values of y. The quantity m is to be 
determined from an inspection of the observations, and if it 
appears that several values of m might be used, all should be 
tried, and that which leads to the smallest sum for the squares 
of the residuals is to be considered the best. If the several 
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cycles are not similar and regular, additional terms involving 
multiples of x will have to be added to equation (28). 

43. Special Treatment of Exponential Equations. If the 
equation selected is not in the linear form as regards the 
unknown constants the general method of procedure is 
given in paragraph 44, but in some special cases a more 
simple reduction is possible. A quite common case is the 
following: — 

Suppose the relation between x and y is expressed by the 
equation 

y = fcx" (29) 

the problem being to obtain the best values of k and m. 
Taking the logarithms of both members of (29), 

log y = log k-\-m log x 

Denoting log k by k', this equation becomes 

m log x-\-k' — log y ^ 

which is in the linear form as regards the unknowns m and 
k', and the normal equations may now be formed in the usual 
manner. 

The most convenient way of determining whether equation 
(29) is a suitable one to select or not is to plot the corre- 
sponding values of x and y on logarithmic cross-section 
paper. If (29) is a proper equation to use the plotted points 
will lie upon a straight line. 

REDUCTION OF OBSERVATION EQUATIONS TO THE 
LINEAR FORM. 

44. When the observation equations are not linear as 
regards the unknowns, the only practicable method of pro- 
cedure, as already mentioned in paragraph 23, is to reduce 
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them to that form. This reduction may be effected as 
follows : — 

Let the observation equations be 

/i(Z„ Z„ . . . Zj) = J»f, 

(A) 



/„(Zl Z^,... Zg) = M. 

in which Z^, Z^, . . . Zg represent the unknown quantities and 
3Ii, 3/^, . . . Mn the observations, the functions being of 
known form. 

Let Zi', Zj', . . . Zg' be approximate values of Zj, Z^, . . . Zg, 
found by trial or the solution of a sufficient number of the 
observation equations, and let the most probable values of 
Zi, Z2, . . . Zg be 

-Zi' + si, Z,' + 3„...Zg' + 3g (B) 

Si, 82, . . . Sg being small corrections whose values are to be 
determined by the Method of Least Squares. The first 
observation equation in (A) may then be written, 

/i(Zi' + 8i, Z^' + s^, • • • Zg' + h) = M, 

Expanding the first member by Taylor's Theorem, denoting 
/i(^i') ^2) • ■ ■ Zq) by ^1, and neglecting terms containing 
powers of Zj, z^, . . . z^ higher than the first, this becomes 

If now we represent the coefficients of z^, z^, . . . Zg, by 
«!. hy ■ ■ ■ ^i) also ki — Ml by m^, and treat the other 
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equations in (A) in the same manner, our observation equa- 
tions will take the form 

«l3l + *lZ2 + • • • qi^q + Wi = Wi 
«22l + ^2^2 + ■ ■ • q^Zq + Wj = ^2 

(C) 

«n2i -|- 5„Z2 + . . . qnZq -f W2„ = V„ 

The second members reduce to the residuals since the 
quantities z^, z.,, . . . Zg represent merely the most probable 
values of the corrections. Equations (C) can now be solved 
in the usual manner and the most probable values of 
Zi, Zj, . . . Zq found by substituting the results in (IJ). 

45. Example. From the following observations find the 
most probable values of x and y. 

sin X -\- cos 2y = 1.5 
cos a; -|- 3 sin y = 1.7 (A) 

a;^ -f- 5y = 2.1 

By trial it is found that 42° and 18° are approximate 

values of x and y. Then in accordance with paragraph 44 

we put 

a; = 42°+zi y—l%°-\-z^ (B) 

and expanding the different functions, we find 



Ai = sin 42° + cos 36°, 


|^^ = cos42°, 


dk, 

3- = - 2 sin 36° 

dy 

= -1.18 


= 1.48 


= .74 


^2 = cos 42° + 3 sin 18°, 


^4^ = -sin 42°, 


f^ = 3cosl8» 
dy 

= 2.85 


= 1.67 


= -.67 


, /427r\^ 18^ 
*»-U80/+ 180' 


dks „ 42ff 
dx — 180' 


dk. 


= 2.11 


= 1.47 


= 5 


&! - Jfi = - .02, k^ 


- Jfs = - .03, 


k,-M, = M 


= mi, 


= OTj, 


= m. 



36 METHOD OF LEAST SQUABES. 

Therefore we have for 

OBSERVATION EQUATIONS 
.74gj _ I.I8Z2 - .02 = «i 
_ .67zi + 2.8532 — .03 = v^ (C) 

1.47zi -j- bz^ + .01 = -yj 

From these are obtained the 



NORMAL EQUATIONS 

3.I6Z1 + 4.572j + .020 = 
4.57zi 4- 34.5IZ2 - .012 = 



(D) 



Solving, 



.00845 Zi = .00147 (E) 



These results are in circular measure. Reducing to degree 
measure we have 

z, = -29'.0 Sj = 5'.l (F) 

Substituting these in (B), 



X : 



41° Sl'.O y = 18° 5'.1 (G) 



46. If the observations are conditioned, precisely the same 
method will be followed in reducing all the equations to the 
linear form. Tlie rest of the solution will then be as usual. 

If the values found for Zj, z.j, ■ ■ ■ Zg should turn out to be 
so large that the terms involving their second and higher 
powers can not be neglected as assumed, the process must be 
repeated using the values of Zi, Z^, . . . Zg first obtained as 
approximations. 

In a few cases the equations may be reduced to the linear 
form by some special artifice of a simple character, as in 
paragraph 43. In such cases the method of expansion by 
Taylor's Theorem should not be used. 



CHAPTER III. 

THE PRECISION OF OBSERVATIONS. 

47. The work so far considered has treated solely of the 
methods by means of which the most probable values of the 
unknown quantities may be determined from a series of 
observations. But in general something more than this is 
desired. We wish to know, if possible, how much reliance 
can be placed upon the results obtained, and how they com- 
pare in precision with other determinations of the same 
quantities. Preliminary to the discussion of the precision of 
our results it will be necessary to consider more fully the 
Curve of Error. 

48. The Constant h. The law of distribution of errors of 
observation has been shown to be represented by a curve 
whose equation is 

y = ke~^^^ 

In this, if a; = 0, y = k (30) 

Therefore the constant k represents the intercept of the 
Curve of Error on the axis of Y. It is not, however, an 
independent quantity but is determined by the value of A, 
as will now be shown. 

49. To Find k in Terms of h. Since, as shown in para- 
graph 18, the total area between the curve and the axis 
of X is denoted by unity, we have 

k C^^e-^'^' dx=l or k C'^e-^'^^ dx=\ (a) 
This may be written 






\-^^'^''hdx = - (b) 

2k ^ 
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Let t = hx, .-. dt ^ hdx. Also, when jb = oo, < = tx, 
and when a; = 0, < = 0. 



r e-^"" dt = r e-'l^^'' h dx 

^0 «/0 

Multiplying this equation by 



(c) 



we have 



[j^ e-^' ^« T = J" j^°° e-'^'d + ^^ h dx dh 



dh 



~2J<, 1 



dx 



-\-x^ 



tan ^ X 



^ a 



IT 

4 



e dt ■= — 



Substituting (31) in (c) and (b) 

Therefore equations (16) and (17) become 

Vir 



(31) 



(32) 



(33) 
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50. It has been shown that the quantity A is a "Measure 
of Precision" of the observations, and hence a determination 
of its value in each case would enable us to compare the 
relative reliability of different measurements. In practice, 
however, it is not found convenient to compute the value of 
h directly, and so this quantity is used only in developing 
the theory of the subject, while in the comparison of observa- 
tions certain other quantities now to be derived are used as 
precision measures. These latter quantities are called the 
"Mean of the Errors," the "Mean Error," and the "Probable 
Error," respectively, and may be computed directly from the 
observations. Further, as it will be found that they all bear 
a definite relation to A, the value of this quantity can be 
determined from them if desired. 

In the following discussions no distinction is made between 
positive and negative errors of the same numerical magnitude, 
and unless otherwise stated the observations are all of the 
same weight. 

THE MEAN OF THE ERRORS OR AVERAGE DEVIATION. 

51. The Average Deviation (a.d.) of an observation is 
the arithmetical mean of the errors all taken with the positive 
sign. 

Since from (6) the probability that the error of a single 
observation will fall between x and x -\- dx is 

4>(x) dx, 

if n observations are made, the number of errors falling 
between these limits is 

n (j>(x) dx. 

Hence the sum of all the errors of the observations is 
n C X 4,{x) dx or 2n C x <t>(x) dx. 
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Dividing this last expression through by n we have 



a.d. =: 2 I x<f>(x) dx = — r i 



e-^^^ X dx 



= fe-'''*' (-2h^x) dx 



or a.d. = — - (34) 



THE MEAN ERROR. 

52. The Mean Error (/a) of an observation is the square 
root of the arithmetical mean of the squares of the errors. 

The total number of errors being n, the number falling 
between x and x-\- dx is, as just shown, 

n (t>(x) dx, 

and the sum of the squares of these errors is 

n x° 4>(x) dx. 

Therefore the sum of the squares of all the errors is 

n I x^ <t>(^) dx 
••• 11" = — f" e-'i^' a!= dx (a) 

But as shown in paragraph 49, 

A r" e-A^x^ dx=l or r e-h"^' (fa; = ^ 
Differentiating (b) with respect to A, 



(b) 



.2A 



n/; 



r e-'><^'^^ x^ dx = -1^ (c) 
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end replacing the integral in (a) by its value as determined 
from equation (c), we have 



THE PROBABLE ERROR. 

53. The Probable Error (r) of an observation is an error 
such that one-half the errors of the series are greater than it 
and the other half less than it. Or it is an error of such a 
magnitude that the probability of making an error greater 
than it in any given observation is just equal to the probability 
of making one less than it, both probabilities being one-half. 

The probability that the error of an observation will fall 
between x and x-\-dx being <l>(x) dx, the probability 
that the error will fall between the limits r and — r is 



P = r,l>(x)dx = ^ fe-^'^^'dx (36) 



If r is the probable error, P is one-half, or 



- Ce-f^^^dx = 1- (37) 



ajd from this definite integral r is to be found. 

Let t = hx, .: dt = hdx. Also when x^r, v,'e have 
t =^ hr, and when a; = 0, t=zO. Substituting these results 
in (37), we have 

- C^e-^dt = i (38) 
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Denote hr by p. Then by interpolation in a table of 
values of this integral, the value of hr in (38) is found to 
be 

p = hr = .47694 (39) 

54. When t is small the values of j e"'^ dt may be 

found by expanding e~' into a series and integrating the 
successive terms. Thus, by Maclaurin's Theorem, 

Jo Jo ^ [2_ [3_ ^ 

3 5[2_ 7[3_^ 



RELATIONS BETWEEN |l, r, «.d., ft, AND p. 

55. From (40) r = ^, and from (34) a.d. = , 

^ hV^ 



r ^z p a.d. v' 



TT 



r = .8453 a.d. (41) 



Also from (35) p. = , 

Av/2 



r = ppSl2 
or »• = .6745 |i (42) 
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The relation between the values of /i, r, a.d., and h may 
be conveniently expressed as follows : — 

1 r 

,iy/^ = - = - = a.rf. 0F (43) 



or arranging in 


tabular form, 








H- 


r 


a.d. 


/» 


1.0000 


1.4826 


1.2533 


r 


0.6745 


1.0000 


0.8453 


a.d. 


0.7979 


1.1829 


1.0000 



From this it will be seen that 

/x > a.d. > r (44) 

56. Further, since in (25) it was shown that p oc h'^, it 
follows at once from (43) that 



1 1 1 

a.d.^ "^ ^^ "^ r"" 



P ^ I-n "^ —^ "^ ~. "^ ^' (*5) 



That is, the weights of different determinations of a 
quantity vary inversely as the squares of their Mean Errors, 
their Probable Errors, or their Average Deviations. 

It is to be observed, however, that the determination of the 
relative weights of quantities from a comparison of their 
precision measures according to (45) applies only when the 
quantities are of the same kind and subject to the same con- 
stant errors, if any of the latter exist. (See § 98.) The 
applications of (45) are numerous and important. 

57. Example A. Suppose n direct observations, all of 
the same weight, be made upon a quantity, and that the 
probable error of a single observation is r. Then since the 



1 




r 


— 


or 


'•o = -n: 


n 




\Jn 
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weight of the arithmetical mean is n, its probable error r^ 
will be given by 

2 1 «• 

(46) 

Or in general, suppose 8 is any precision measure of an 
observation of weight p, and suppose Po is the weight of a 
second similar quantity or observation, then the correspond- 
ing precision measure So of the latter will be 

8„ = 8^^ (47) 

The case of most common occurrence is that in which 
jo = 1, and then we have 

80=7^ (48) 

Example B. A line is measured five times and the aver- 
age deviation of the mean (A. D.) found to be .016 feet. 
How many additional measurements are necessary in order 
that the A. D. of the mean may be reduced to .004 feet ? 

Let X be the total number of observations required. Then 

a; : 5 = .000256 : .000016 
a; = 80 

Consequently the number of additional measurements 
required is 75. 

Example C. In two determinations of the quantity X 
there were obtained 

Zi = 427.320 ± 0.040, Zj = 427.30 ± 0.16 

Find their relative weights, and the most probable value 
of Z and its probable error. 
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Note. The above is the method commonly employed to 

denote that the probable errors of the observations are 0.040 

and 0.16. 

Pi 16^ 16 

From (45) — = = — 

From (3), the most probable value of i will be given by 

Z„ = 427 + '' X -f^ + ■'' 
= 427.319 
and the weight of Zq being 17, by (48), 

.16 
To = -— = .039 

VI 7 
Therefore we should write the result 

Zo = 427.319 ± .039 



REPRESENTATION OF Jl, a.d., AND r ON THE CURVE OF 

ERROR. 

58. To find the points of inflection of the Curve of Error 
we have 

ax 

dx' 

d'^y 
For a point of inflection, y— j = 0, |or 

2 A^ A e-'''*' (2 A^ a;^ - 1) = 

X = = M by (35). 

hsl2 '^ ^ ^ ' 
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That is, the Mean Error is represented by the abscissa of 
the point of inflection of the Curve of Error. See OM in 
figure 2, page 9. 

Next, for the abscissa of the centre of gravity of the area 
to the right of OY, we have 

I y X dx 

X,=:OD = ^ 

^ y dx 

2 h 

e-'''^ X dx 



\/7 Jo 



for, (§ 18), f'y dx=^ 
Integrating, 



ajo =: — — = a.d. by (34), 

h\7r 



Finally, if an ordinate PJ" be drawn so as to bisect the 
area to the right of the origin between the curve and the 
axis of X, the Probable Error will be represented by the dis- 
tance of this ordinate from the axis of Y. For the proba- 
bility of the occurrence of an error less than the amount OP 
is then equal to the probability of the occurrence of an error 
greater 1;han OP. This being the case, OP is the probable 
error by definition. 



CHAPTEE IV. 

COMPUTATION OF THE PRECISION MEASURES. 

DIRECT OBSERVATIONS. 

59. Observations of Equal Weight. Given n direct obser- 
vations all of the same weight on a single quantity M, to 
find the Mean and Probable Errors and Average Deviation 
of a single observation and of the Arithmetical Mean. 

Let the observations be Mi, M^, . . . M^. 

" " arithmetical mean be M^. 

" " real errors be iCj, x^, . . . a:„. 

" " residuals be Vy, Wj) • • • '"n- 

Denote the mean and probable errors and average devia- 
tion of a single observation by fi, r, and a.d., respectively, 
and the corresponding quantities for the arithmetical mean 
by fioj '"o) and A.D. Then by definition, 

Wi = JIfi — ibTo, Wj = Jlfj — Jl/o. • • • «» = -^ — Ma 



and 






If Mo represented the true value of M, the residuals 
would be the same as the real errors, and we should have 



=V' 



n 



and if n is large this formula is practically exact. But when 
n is small a more accurate expression is necessary. To 
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obtain this let M^ + x^ be the true value of 31. Therefore 

ail = Jfi - {Mo + Xo) — Vi — Xo 

Kj = J/j — (Mo + Xo) = V^ — Xo 
Xn = Mn - {Mo + Xo) = Vn ~ Xo 

Squaring, adding, and dividing by n 

^ = /^= = - (2«= - 2 Xo^V + 7W) 

= ^ + Xo' for by (2) 2» = 

The value of Xg is not known and can not be found 
exactly, but it is approximately equal to the mean error of 

Mo, that is, by (48), to /aq = Ji.. 

Substituting this value in the above we have 

'^ n ^ n 

{n - 1),.= = Su' (49) 

(50) 

^y(^'> ^" = V,^fzT) (51) 

by (42) r = .6745 \/-^ (52) 



and 



V- 


= 


^ n—l 




V-0 


^ n(n — 
.6745 \/- 

.6745 \] - 


1) 


r 


»t- 1 


To 


tin. _ n 



(53) 
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60. In order to avoid the use of the squares of the 
residuals we may proceed as follows : From (49) 

S«» = (n - 1) — 
On the average, the values of the residuals will then be 



In - I 
In— 1 



\n - \ 



!B» 



Adding and dividing by n, neglecting the signs of the 
residuals, 

n \ n n \ n 

a.d. = — (54) 

y/n (n — I) 

by (48) A.D. = (55) 

n \n — 1 

.8453 Sv ,-„. 

by (41) r = (56) 

\/n {n — I) 

.8153 S« ,.,, 

and n = (57) 

n Vn — 1 

The mean errors may also be computed from the above by 
using the table of equivalents in paragraph 55, but this is 
not customary, formulas (50) and (51) being used for this 
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purpose. Results derived from (50) are to be regarded as 
more accurate than those obtained from (54), the latter 
being a second approximation. 

61. Example. From the following measurements on the 
length of a base line find the most probable length and the 
values of the various precision measures. 



M 


V 


v' 


455.35 


.02 


.0004 


.35 


.02 


4 


.20 


- .13 


169 


.05 


- .28 


784 


.75 


.42 


1764 


.40 


.07 


49 


.10 


- .23 


529 


.30 


-.03 


9 


.50 


.17 


289 


.30 


- .03 


9 


3.30 


+ .70 


.3610 


455.830 


-.70 




M, 


Sw= 


2y" 



By (50) and (51) 
^/.3610 

r = .6745 ^ 

By (54) and (55) 

1.40 
a.d. = 



s/90 
r = .8453 a.d. 



.20 
.13 

.15 
.13 



V'lO 
To = .6745 ,.„ 



.063 
.042 



A.D. = 



1.40 



30 



.047 
Ta = .8453 A.J). == .040 
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and we should write for the most probable length of the 
base line 

Mo = 455.330 ± .042 

62. Observations of Unequal Weight. Using the same 
notation as above, with slight modifications, we will 

Let 3fo represent the General Mean. 
" i'u Pi-, ■ • • Pn be the weights. 
" a.d-xi a.d.^, . . . a.d.^ be the average deviations. 
" fii, fu, . . . /j^ be the mean errors. 
" ''15 ^21 • ■ • ''ii be the probable errors. 
" a.d., fi, r, and v refer to observations of weight unity. 



etc. 



If the " Observation Equations " are formed for this case 
they will be 

Ml — Mo = Vi, M2— Mo = v^, . . . M„ — Mo = Vn. 

And, as was shown in paragraph 28, if these equations are 
each multiplied by the square root of the weight of the corre- 
sponding observation, they will all be reduced to equivalent 
equations of weight unity. On performing this operation 
it will be seen that the residuals of the new equations 
become 



Then by (48) 






Mo = 




SPk 



fl 



'fPv, ^s V^, • • • ■"n V^ 



And evidently to these reduced observations the formulas 
of paragraph 59 apply. 
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Therefore 



= J ^^^' r = .6745 V- (58) 

V n — I 



. ^ J _1P^ r, = .6745 It, 



(59) 



^^ = i/ ^ f^^' ^ n = .6745 iJio (60) 

Also, by a method similar to that used in paragraph 60, it 
may be shown that 

a.d. = —^^^=- r = .%ma.d. (61) 

S n{n — 1) 

a.d.^ = ^^^^ — ^^ ^ g^gj ^^^ (g2) 
vi>i**(n - 1) 



4.1>. = ^^^ — n = .8453 A.D. (63) 
\l'S,pn{n — 1) 



Formulas (61), (62) and (63) are not of much value in 
practice unless all the weights are perfect squares, for other- 
wise no labor is saved in the computations, since the square 
root of each weight will have to be determined, and the results 
obtained can not be considered as reliable as those given by 
formulas (58), (59) and (60). 
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63. Example. Given a series of observations on 31, the 
difference in longitude between two stations. To find Mo, 
11^ fio, a-d., A.B., r, r^, r^. 



M 


P 


pM 


V 


v\jp 


y^ 


pv^ 


V 4' 30" 


4 


120 


1.6 


3.2 


2.6 


10.4 


41 


1 


41 


9.4 


9.4 


88.4 


88.4 


43 


1 


43 


11.4 


11.4 


130.0 


130.0 


37 


9 


333 


5.4 


16.2 


29.1 


261.9 


48 


4 


192 


16.4 


32.8 


269.0 


1076.0 


34 


16 


544 


2.4 


9.6 


5.7 


91.2 


25 


9 


225 


6.6 


19.8 


43.6 


392.4 


46 


1 


46 


14.4 


14.4 


207.4 


207.4 


28 


25 


700 


3.6 


18.0 


13.0 


325.0 


24 


4 


96 


7.6 


15.2 


57.7 


230.8 


1° 4' 31".6 


74 


2340 




150.0 




2813.5 


Mo 


2/> 


Sjom 




2 v\fp 




%pvi^ 



%pM 2340" 

By (3) Mo = =f::- = 1° 4' + 



-s/' 



2813.5 



9 



ip ' 74 

= 1° 4' 31".6 



= 17.7 ,10 = -^. = 2.1 

V'74 



r = .6745 u = 11.9 



150 

a.d. = -^ = 15.8 

V90 

, = .6745 ^0 = 1-4 



r< = - = 4.0 
* 3 

A.I). = ^= 1.8 

V/74 

.8453 A.I>. = 1.5 



Mo = r 4' 31". 6 ± 1".4 
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FUNCTIONS OF OBSERVED QUANTITIES. 

64. Theorem. Given any number of quantities and their 
Mean and Probable Errors and Average Deviations, to find 
the Mean and Probable Errors and Average Deviation of any 
function of the quantities. 

Let the quantities be M^, M^, • . ■ Mg. 
" mean errors be /ij, /t,, . . . fig. 
" function be Mz=f(Mi,Mi, . .. Mg). 

" mean error of Jf be £J. 
" probable error of M be H. 
" average deviation of M be D. 

The derivation of the general formula will be simplified if 
we consider first a few special forms of functions. 

65. Case I. Suppose M = Mi ± Mi- 

The number of observations from which Mi and M^, and 
hence ju,i and fi^, have been determined is not necessarily 
known, but we may assume that for each quantity it is any 
large number n, and that the real errors of the observations 
are 

for Ml, xi', Xi", xi'", . . . 

-'"2> •''2 ) •''2 ) •''2 > • • • 

Then the real errors of M, computed from the separate 
observations on My and M^ will be 

iCi i: ajg , ajj it iCg , Xi "x. x^ ... 

n 

_ Siei'± 2 2g!ia!i, + Saia' 
n 

or W = \ii' + 11,^ (64) 

since in the most probable case the term Sa^asj will dis- 
appear, as there most likely will be as many positive as 
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negative products of the same absolute magnitude of the 
form Xx Xj. 

By successive applications of the above the same principle 
may be extended to cover the algebraic sum of any number 
of quantities. So that if 

M =z M,±M,±... Mg 
then JE^ = jij^ + ji,^ + . . . ji^^ = S|Ji'' (65) 

Since probable errors and average deviations differ from 
mean errors merely by a constant factor, we shall likewise 
have 

B^ = r,'' + r,'' + . . . r/ = 2*-= (66) 

J>2 = a.d.^^ + a.d.^^ + . . . = l^a.d.^ (67) 

66. Example. Given the telegraphic longitude results, 

A. TO. sec. sec. 

(A) Cambridge west of Greenwich, 4 44 30.99 ± 0.23 

(B) Omaha " " Cambridge, 1 39 15.04 ± 0.06 

(C) Springfield east " Omaha, 25 8.69 ± 0.11 

Find X, the longitude of Springfield, and its probable error. 

X = ^ + i? — C = 5* 58"* 87'.34 ± 0'.26 
for by (66) H = v/(.23'' + .06^ + .IP) = .26 

67. Case II. Suppose 31 = aiJfi 

Using the same notation as in paragraph 64, the real 
errors of M will be 

E' = ^^^^ = a,W 
or B = ttiVi (68) 
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68. Case III. Suppose Jf = aiMi ± a^Mi. ■ ■ . a^Mq 
By combining (65) and (68) 

also n^ = l.aW^ (69) 

and J>= = 'Za'a.d.^ 

69. Example. The length of a bar at 20° Centigrade is 
found to be 75".0041 ± 0".0037, and its expansion per 
centigrade degree 0".0036 ± 0".0018. What is the length 
of the bar at 56° Fahrenheit ? 

20° C. = - X 20 + 32 = 68° F. 
5 

5 5 

The expansion for 1° F. will be - X .0036 ± - X -0018 

y y 

or .0020 ± .0010 

L = 75.0041 - 12 X .0020 
= 74.9801 



Ji = V-0037^+ (12 X -OOl)" 

= .013 
i = 74".980 ± 0".013 

70. Case IV. Suppose M = f {M„ M^, . . . Mq). 

Let Jfi = «! + mi, M^ = a^ -{- m^, . . . Mq = aq -\- nig, 
where a^, a^, . . . a^ are arbitrarily assumed quantities very 
nearly equal to M^, M^, . . . Mg, so that m^, m^, . . . nig are 
so small that their second and higher powers may be neg- 
lected. Then the errors of M^, If^, . . . Mg may be con- 
sidered to be in wij, m^, . . . nig, and hence /ixj, ^, . . . ^g 
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may be regarded as the mean errors of the quantities 
wii, mj, ... wig. We now have 

-^ = / («i + »»i, ttj + Wj, . . . aj + m,) 

Expanding this expression by Taylor's Theorem, and 
denoting /(«!, a^, . . . Og) by ^', we have 

M = M' 

dM' , 5Jf' , dM' 

-(- negligible terms in the second and 
higher powers of wij, m^, . . . m^. 

Then the mean error of M will be the same as the mean 
error of terms (A), and by (69) this is given by 

or what is practically the same thing, by 






71. Example A. Two sides of a right triangle are meas- 
ured with results 

a = 49.53 ± 0.59 
* = 50.38 ± 0.93 
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Find the length of the hypotenuse c and its probable 
error. In this example 



M= ^Mi'+ M^^ = sja^ + b^ = 70.65 

BM __ a BM b 

BM^~ y/a2_|_ 62' BM^ ~ yja^ J^ ¥ 

By (70) R^= V^^, + '^"^^^ 

■=. .59' + .93^^ 

|_70.65j ^ l_70.65j 

R — .78 
and c = 70.65 ± 0.78 

Example B. If the probable error of x is r, what is 
the probable error of the common logarithm of a; ? In this 
case 

M = logio Ml = logio X 

BM B logio «■ logio e 

B3fi Bx X 



Example C. If the weight of x is p, what is the weight 
Po of sin X ? 

Denoting the mean errors of x and sin a; by yu, and E, 
respectively, we have, by (45) and (70), 

•^„ o r^ sin a;"]' 

E^ = li!^ \ = iL^ cos= X 
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and - = 4 

P 
Pa ^ = p sec' X 

cos^a: 

72. Equation (70), which expresses the law of propaga- 
tion of error in functions of observed quantities, is one of the 
most important in the whole theory of the Method of Least 
Squares. Upon it in particular is based the discussion of the 
" Precision of Measurements." This subject treats, in the 
first place, of the methods of finding the precision of a 
quantity obtained by computation from a series of measured 
quantities ; and, in the second place, it investigates the pre- 
cision with which the component measurements of a series 
must be made in order to obtain a required degree of pre- 
cision in the final result. The following simple example will 
illustrate the character of the solutions : — 

73. Example. In the determination of a current by a 
tangent galvanometer we have 

H 

I z=z 10 — tan <^ 
G 

where I is the current in amperes, H the horizontal com- 
ponent of the earth's magnetic force, G the galvanometer 
constant, and <^ the angle of deflection. Given the errors 
8i, §2, hi va H, G and tan </>, to find the error A in I. 

By (70) 

100 100 ^» „ „ , 100.B"' „ 

= — 8i^ tan^ ,^ -f —^^ 8,^ tan« 4, + —^ 8,= (a) 
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1 no TT^ 
Dividing this equation by I^ = tan^ <^ we have 

That is, the square of the percentage error in I is equal 
to the sum of the squares of the percentage errors in S, G, 
and tan (f>- Hence if 

IT is determined within .4 per cent, 
tan <^ " " " .1 " " 



then 



— = v' .16 + .04 + .01 = .46 per cent. (c) 

Next, suppose the value of I is required to within .1 per 
cent. To find the necessary accuracy in the determinations 
of S, G, and tan <f> when the error in each of these quan- 
tities is to have the same influence upon the total error. 

From (b) we shall now have 



8i 


= 


.001 

= .000577 

f3 




8i 


= .00058^ 




8. 


= .00058 G 




83 


= .00058 tan <t> 



(e) 



It is comparatively easy to obtain the necessary accuracy 
in the measurements of G and tan 1^, but diflScult in the 
case of If. 

For additional work of this kind see Holman's " Precision 
of Measurements." 
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74. Combination of Functions of the Same Variables. 
It is to be noticed that equation (70) applies only when M 
is a function of independent quantities. If iHf], M^, . . . Mq 
are merely different functions of the same quantities we 
must proceed as follows : — - 

Let Ml = (j, («!, Sj, . . . z„) 

M^ = xj, (2i, S2, . . . Si) 

M = f {Ml, M,) 

If any single observations of Sj, z^y • ■ ■ ^k ^^^ subject to 
errors Xj, x.2, . . . x^, the corresponding errors in Mi and 
M2 will be 

for Ml, Xi = aiCCi -|- a^x^ + . • . aja:* (a) 

" Jfj, X^ = a/aji -f- aj'xj + . . . ai'a;* (b) 

Where Oi, a„, . . . a^ are the differential coefficients of 
Ml, and a/, «„'... a^.' the differential coefficients of M^ 
with respect to Zj, z^, . . . Zt- The corresponding error in 
M will then be 

X = AXi + ^'JT, (c) 

Where A and ^' are the differential coefficients of M 
with respect to Mi and M2. Substituting in (c) from (a) 
and (b), 

X = (Aui -\- A'a'i) Xi -f (Aa^ + A'a'^) x^ . . . 

= aXi -\- ^X^ -\- . . . kXjc 

Then if the number of observations or values of X be 
denoted by n, 

E^ — 



^X" 
n 


= 




+ P' 


2a;2= 
n 


+ 




+ 


^ '^ + 






= , 


^W 


+ p^l^' 


+ .., 


. xVv 


z 



^» = aVi^ + )3»^^ + . . . XVt^ (d) 
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since in the most probable case the product terms will 
cancel out. 

Expanding (d) we have 

E^ = {Aa^ + ^'a/)^ ix.^"- + {Aa^ + A'a^y mj' + • • • 

= A\a,W + «2W ■ • •) + ^'\<h!W + «»'W • • .) 
+ 2AA'{a^a^t,^^-\-a^a^^^-\-...) (71) 

75. Example. As a very simple problem take 

ilfi = 22„ Jfj = Szi, tx.^ = 0.1, 
and JIf = J!fx + M^. 

Then ^ = 1, ^' = 1, «i = 2, ai' = 3. 
By (71) ^2 _ 4 ^ 01 + 9 X .01 + 2 X 2 X 3 X .01 = .25 
or E = 0.5 

In this particular example the result may be found directly 
from (68) by substituting at first in M the values of Mi 
and M^. 

Thus M = 23i + 3si = 53, 

E — bfii = 0.5 

If il/i and M^ had been independent quantities, by (64) 
or (69) we should have had 

E = ^(2 X 0.1)2+ (3 y 0.1)= 
= 0.36 

INDIRECT OBSERVATIONS. 

76. The determination of the precision measures of the 
unknown quantities in case the observations are indirect 
involves a knowledge of the weights of the unknowns, and 
consequently the method of computing these weights must 
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first be demonstrated. It will be assumed at present that all 
the observations are of weight unity. 

FIRST METHOD OF COMPUTING THE WEIGHTS. 
Let the observations be 

ajZi + h^z^ + • • • yaSg + ^"2 = J^i 

«nSl + K'^i + • • • 2rfiq + ^"n = ^n 

in which M^, M^, . . . Mn denote the actual observations, 
and 2i, Sj, ... Sg the most probable values of the unknown 
quantities. Let 

^1 — Jlfl = Wl, ^2 — Ji^ = OT2, . . . ^n — ii^ = Wl„ 

Then the above equations give rise to the 

OBSERVATION EQUATIONS 
«i3i + h'^i + • • • ^i^s -\- m^ = Vi 
a^Zi -\- b^Zi + • • • 2'2Zg + mj = Uj 

(A) 

On«l + K^". + • • • irfii + "*» = ■"» 

By the rule, paragraph 25, we now form the 

NORMAL EQUATIONS 
Zi S a'' + 22 2 a6 4- . . . Zg 2 ag' + 2 am = 
21 2 a6 + S2 2 6^^ + . . . Zg 2 5? + 2 im = 



gjSaj+ 2,26?+ ... 232?" + 2?m 



= 
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Multiply the first of (B) by Qi, the second of (B) by 
§2, . . . and add the results. Then 



si(§i2a=+ Q.^ab + . 


.. Qq2aq) 


+ s, (^iSaJ + C226' + . 


.. Qg^bq) 


+ MQ,%aq+ Q.-Zbq^ . 


■ ■ Q,^q') 


-|- §i2 am -\- Q2%bm -]- . . 


. Qq%qm 


= 





(C) 



Let Qi, Q2, ... Qq be determined so that the coefficient 
of Zi in (C) shall be unity, and the coefficients of Zj, 
Sj, ... Zg each equal to zero. That is, let 



Ci2 a= + §jS aJ + . . 


, . Qgl,aq - 


1 = 




e,2a&+ Q,-%b^ + .. 


. Qq^bq 


= 


(D) 


$,2a?+ e^SS? + .. 


■■ QqM' 


= 




Then (C) becomes 









Sv + Qi^am + Qi^bm + • • • Qq'^qm = (E) 

Equations (D) may be derived from the normal equations 
(B), if in them we replace z„ z^, ... s, by Q^, Q^, . . . Qg, 
'S,am by — 1, and 2 5m, 2 cm, . . . 2 2'm by zero. Hence 
the solution of the normal equations with these changes will 
give the values of Qi, Q^, ... Qg. 

TO SHOW THAT Qi IS THE RECIPROCAL OF THE WEIGHT OF Xi. 
Expanding the coefficients of (E) we have 
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+ C2(<^i»Wi + JjWla + . . . *n»n„) 

4- Qq{Qi>^l + ?a»"j + • • • ?nW„) 
= 



or collecting the coeflBcients of mj, m^, . . . rrin 

Sl + »Wl( Ql^i + 62^1 + ■ • • Qq2\) 

4- mj( ^ittj + §262 + • • • Cgg'i) 

= 



For convenience in writing we will let 



Oj 



= 6l«2 + 62*2 + 



Qq2l 
Qqii 



(E') 



a„ = Cl«n + C2*» + • ■ • Qq2n 
Zl 4- <»i»ni 4- Ojmj + • • • <»n«»n = 

Multiply the first of (F) by a^, the second by 
and add the results. Then 

2 oa = CiS «^ + ^22 aft + . . . Qq% aq 
= 1 by first of (D) 

Multiply the first of (F) by *i, the second by 
and add the results. Then 



(F) 
(G) 



(H) 



K 
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= by second of (D) 

Likewise S ca = (H') 



S^a = 

Multiply the first of (F) by aj, the second by oj, . ■ . 
and add the results. This gives 

2a» = CiS Oa + ^jS Sa + . . . §,,2 qo. 

= Q^ by (H) and (H') (I) 

Let jj. be the mean error of an observation of weight unity. 
Let /i.2, be the mean error of Sj. 
Let Pz^ be the weight of Zj. 

The mean errors of wii, Wj, . . . m„ are the same as the 
mean errors of M^, M^, . . . Mn and each is accordingly 
equal to /i. Therefore from (G), by (69) 

Mzi" = «iV + "sV + • • • o«V 
= /i,^2a^ 

= q,f? by (I) (J) 

But by (48) ^^^ = ^ 

Comparing this with (J) we see at once that 

Therefore, for the First Method of computing the weights 
we have the following : — 
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77. Rule I. In the normal equation for Zj write — 1 
for the absolute term % am, and in the other equations 
zero for each of the absolute terms 2 bm, 2 cm, ... 2 qm,. 
The value of 2i found from these equations, is the recipro- 
cal of the weight of the value of z^ obtained by the solution 
of the normal equations. 

To find the weights of z^, z^, . . . Zg, proceed in a similar 
way, forming a corresponding set of equations for each 
unknown. 

SECOND METHOD OF COMPUTING THE WEIGHTS. 

78. Write equations (B) of paragraph 76 in the form 

Si 2 a° + ^2 2 aft + ■ • • 83 2 a? -|- 2 am = A 
zi 2 aJ + Z2 2 S*" + . . . 2„ 2 5? + 2 5m = £ 

(L) 

Zi'S.aq -\- Zi'S,bq -\- . . . Zgtq^ -^ %qm = Q 

Then in the solution by the preceding method, equation 
(E) becomes 

3i + CiS am -|- Qi% bm -\- . . . Qq% qm 

in which, as was proved in (K), Qi is the reciprocal of the 
weight of Zi. Whatever method of elimination is employed 
in the solution of the normal equations, the coefficient of A 
in the value of z^ must necessarily be always the same. 
Hence we have 

79. Bnle II. Write A, JB, . . . Q instead of zero in the 
second members of the normal equatiojis and carry out 
their solution in any convenient way. Then the most 
probable values of z^, z^, . . . z are given by those terms 
in the results which are independent of A, B. . . . Q. 
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The weight of z^ is the reciprocal of the coefficient of A 
in the value of z^. The weight of Zj is the reciprocal of 
the coefficient of £ in the value of Zj, etc., etc. 

THIRD METHOD OF COMPUTING THE WEIGHTS. 

80. From the second, third, . . . equations of (L) find 
the values of Zj, s„ ... Zg in terms of Zj and substitute 
in the first of (L) without reduction. Then the first of 
(L) becomes 

iJzi = T -\- A -\- terms m B, C, . . . Q 

Where T is the sum of all the numerical quantities result- 
ing from the substitutions. Dividing through by M, 

T A 

Zi = - + - + terms in B, C, . . . Q (N) 
Ji JtC 

T 
in which — is the most probable value of Sj, and, as was 
B 



shown 


in deriving 


the second 


method, 






A 
B 


= QiA 






1 

Ir 


= <2. 




••• 


B = 


= P^ 


From this follows at once 





(O) 



81. Bule III. Substitute in the normal equation for z^ 
the values of Sj, z,, . . . Zg in terms of Zj as found from 
the remaining equations. Then before freeing of fractions 
or introducing any reduction factor, the coefficient of z^ 
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in this equation is the weight of the value of z^ obtained 
in the solution. 

To find the weights of Sj, S3, . . . Zq, proceed in a similar 
way with the normal equations for each of these uiiknowns. 

For the solution of an example by the three different 
methods see paragraph 84. 

THE MEAN ERROR OF AN OBSERVATION. 

82. The next step will be to derive /x, the mean error 
of an observation of weight unity. In the following demon- 
stration the equations referred to by letters are those in 
paragraphs 76 to 81. 

Let the real values of Zj, Sj, ... z^ be 

Sl + «1, ^2 + SBa, • • • 2g + «« 

and substituting in (A) we have 

«i(3i + ail) + Si(s2 + «2) + • • • ?i(2s + ««) + wi = Ai 
aa(3i 4- ^i) + M^^ + a;^) + • • • ?2(«5 + ^g) + WJ^ = A, 

(P) 

«n(Si + a^) + *«(«2 + a!j) + . . . 2„(Sg -f x,) + m„ = A„ 

where Aj, Aj, . . . A„ are the real errors of My^, M^, . . . Mn, 
or of nil, "*2» • • • "*»■ 

Multiply the first of (P) by ai, the second by a^, . . . 
and add the results. This gives 

Si S a^ + z^'S.ab -\- ... Zg^aq -\- "2, am 
-f a;i2a' 4- ^sSaS -f . . . a^Sag- = SaA 

But by the first of (B) the first line in this equation is 
equal to zero, and therefore 
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ail 2 a^ + aja 2 aJ + . . . Kg 2 ay — 5 aA = 
Also ail 2 aft + ajj 2 5^ + . . . a;, 2 J? — 2 *A = 

(Q) 

Xj^taq -\- Xi'S^bq -\- . . . Xq'^q^ — 2 2'A = 

These being of the same form as the normal equations 
(B), the value of x-^ resulting from their solution will be 
of the same form as that of Zj in the solution of those equa- 
tions, with only the substitution of — A for m. From 
(G) we shall therefore have 

ail — ajAi — ajAj — . . . u„A„ = (R) 

Multiply the first of (P) by Wj, the second of (P) by 
1)2, .. . and add the results. Then 

(Si + a^)2aw + (Zj + Xi)%hv + . . . (z, -f- Xq)-%qv 
-\- 'S,mv = 2 Aw 

and multiplying the first of (A) by Oj, the second by 
ttj, . . . and adding the results 

2 aw = Zi 2 a^ + Zj 2 a6 4- . . . Zj 2 ag 4- 2 am 
= by the first of (B). 

Also 2J« =0 

2?w = 

Substituting these values in the above, we find that 

2mw = 2 Aw (S) 

Now multiply the first of (A) by Wi, the second by 
Wj, . . . and add the results. This shows that 
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z^'S.av -\- z^'ibv -\- . . . Zq%qv -\- 'Xmv =z 2 u^ 

and as above, 

2 a« = 2 6u . . . = 2 ?« = 

Combining this result with (S) we have 

•2,mv = %v^ = 2 Aw (T) 

Next multiply the first of (A) by A,, the second by 
Aj, . . . , add the results and compare with (T). This gives 

Si 2 aA + S2 2 JA -(- ... Zg 2 g'A -|- 2 wiA 

= 2 Ay — -Zv^ (U) 

And finally, multiplying the first of (P) by Ai, the second 
by A2, ... and adding the results, we have 

Si 2 «A + zj 2 6a + ... z, 2 3'A + 2 wiA 
-f a;i2aA + a;2 2iA + . . . a:52g'A = 2 A'' 

Therefore from (U) 

%v^ + a!i2aA + X22*A + . . • Xg•S.q^ = 2A^ (V) 

_^ _ 2A^ _ 2^ a;i 2 aA + ■ ■ . tgg 2 gA 

■'• f^ ~ n ~ n '^ n ^ ' 

We must now find the mean values of the terms a;, 2 aA, 
ij 2 *A, ... a;, 2 q^. Expanding 2 aA, 

2 aA ^ ajAi -|- ajAj + • • • a„A„ 
from (R) Xi = aiAj + OjAj + • • • a„A„ 

Multiplying, 

ail 2 aA = OiaiAi" -f- OifH^^ + • • • '"••"n'^n' 
-\- terms in AtA,, AiA„ . . . 
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In the most probable case these product terms vanish, and 
substituting for Ai^ Aj^ . . . A„2 the mean value /x^ we 
have 

iCi 5 aA = /i^ S aa 

= f^' by(H) 

Similarly x, 2 dA = ix^ 



a;, 2 g'A = /j!^ 

From (W) then, q being the number of unknowns, 

, _ 2r>' ^ g^' 
»i n 



11 = V (72) 



n — q 



^ = V,:,!^ (-) 



r, = .674511, = .6745 V^ — r (^i) 



83. By a method similar to that used in paragraph 60, we 
may derive 

a.d. = — (75) 

\] n{n — q) 

a.d.z = -=:^:=z=r (76) 

\l p^n(n - q) 

84. Example. In illustration of the above processes 
take the example in paragraph 24, where we found for 



COMPUTATION OF THE PRECISION MEASUBE8. 73 
OBSERVATION EQUATIONS 



2l — S2 — 1.7 = Vi 

Sj — 2.4 = Vi 
Sj — Zs — 3.0 = U4 



(A) 



and for 



NORMAL EQUATIONS 






2 2i — 2 Sj — Zj _ 0.7 


= 


(a) 


- 2 Zi 4- 3 Zj - 2.3 


= 


(b) 


— Zi _|_ 3 gj _ 0.4 


= 


(c) 



SOLUTION FOR THE WEIGHTS BY THE FIRST METHOD. 

For finding the weight of Zj the above normal equations 
would be written 

2 Zi - 2 Z2 — Z3 — 1 = (a') 

- 2 Si + 3 Zj =0 (b') 

— 3, + 3 s, =0 (c') 

Solve for Zj 

3 X (a') 6zi — 6Z2 — 8Z3-3 = 

(C) 

2 X (V) 



- zi + 3 Z3 


= 


5 Si — 6 Z2 

- 4 zi + 6 Z2 


-3 = 
= 



3 = 
1 



For finding the weight of s^ the equations are 



(d') 



2 Zi - 2 Zj - Z3 =0 (a") 

-2 8i+8zj — 1 = (b") 

- si + 3 z, =0 (c") 
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Solve for Zj 

3 X (a") 6 Zi - 6 Zj - 3 z, = 

(c") - zi + 3 gg = 

5 Si — 6 Zj =0 



5 ' 
Substitute in (b") Zj = - . . p^, z= - (d") 
For finding the weight of Sg the equations are 



2 2, - 2 Zj — Z4 =0 
- 2 si + 3 Zj =0 
_ zi + 3 z, - 1 = 


(a'") 
(b'") 
(C") 


Solve for z. 




4 
from (b'") 2 Zj = - Zj 







substitute in (a'") 2 Zj — 3 Zj =0 

2 X (C") — 2zi + 6z3-2 = 

3 Zs - 2 = 

2 3 

z. = - .. p,, = - (d'") 



SOLUTION BY THE SECOND METHOD. 

The normal equations will now be modified so as to appear 
in the following form : — 

2 Zi - 2 Zj - Z3 - 0.7 = ^ (a) 

- 2 Zi + 3 Zj - 2.3 = ^ (b) 

- Zi + 3 z, - 0.4 = C (c) 
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Solving, 

3 X (a) 6 zi - 6 Sj - 3 2, - 2.1 = 3 ^ 



- «1 






+ 


3Z3 


— 


0.4 : 


= 


C 


5 Zj 


+ 


6 Zj 
622 


— . 


2.5 
4.6 





3 A 
IB 


+ 


C 



zi — 7.1 = 3 ^ + 2 ^ 4- C (8) 

Zi = 7.1 and p,^ = - (d) 

Substituting (8) in (b) 

3 Zj = 16.5 -|-6^ + 5jB + 2C 

3 

Zj = 5.5 and p^ = - (e) 



Substituting (8) in (c) 

3z, = 7.5 + 3^-}-2-B+2(7 

3 

z, = 2.5 and p^ = - (f) 

SOLUTION BY THE THIRD METHOD. 
The normal equations are now taken in their original form. 

2 Zi - 2 Zj - z, - 0.7 = (a) 

— 2 Zi + 3 Zj - 2.3 = (b) 

- Zi -f 3 Z3 — 0.4 = (c) 
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To obtain z^ and its weight we proceed as follows : 

from (c) 2, = 1^ + ^ 

, ,, , 2 zi , 2.3 

from (b) a, = _ 4. _ 

Substitute in (a) 

4 4.6 Zi .4 „„ 

2si 2i L___0.7 = 

3 3 3 3 

z, 7.1 

Collecting terms, — — — = 
o o 

zi = 7.1 and p,, = - (d) 

For Zj 

3 X (a) + (c) 5 Zi - 6 Zj - 2.5 = 

6 

Zi = 7 3j + 0.5 
o 

12 
Substitute in (b) Zj — 1.0 + 3 Sj — 2.3 = 

3 
Collecting terms, - Zj — 3.3 = 
5 

3 
Zj = 5.5 and p^ = - (e) 


For Zj 

3 X (a) + 2 X (b) 2 Zi — 3 Zs - 6.7 = 

3 , 6.7 

o A 7 

Substitute in (c) z^ — 1- 3 Sg — 0.4 = 

^ 2 
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3 7.5 

Collecting terms, - Zj — — ^ 

3 

a, = 2.5 and p^ = - (f) 

It is evident that the three methods give identically the 
same results and that the work is about the same in each case. 
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Substituting the values found for z^, z^ and Sj in the 
observations equations (A), we have 

7.1 _ 5.5 - 1.7 = ui = - .1 .01 = V 

2.5 — 2.4 = «2 = + .1 .01 = V 

- 7.1 + 5.5 + 2.5 - 1.0 = 1)8 = - .1 .01 = V 

5.5 - 2.5 - 3.0 = 1)4 = .0 .00 = v^ 



.03 = ^v" 



In this example w := 4, y = 3. 
By (72) M = \/y = . 



17 



By (74) r = .6745/u, = .12 

By (73) ^, = \J— = .30 r„ = .20 

3 





'-. = \/J = ■'' 






^ = v/f = ■» 




By (75) 


•3 

a.d. = ^ = 


.15 



= .15 



= ,10 



= .8453 a.d. = .13 
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85. Observations of Unequal Weights. If the observations 
are not all of the same weight the formulas and operations 
are merely modified in the usual manner and equations (72) 
and (75) take the more general form 



y n — q ^ 

a.d. = — (78) 

v' n{n — q) 



CONDITIONED OBSERVATIONS. 

86. Suppose there are given n observations, n' condi- 
tions and q unknown quantities. Then by paragraph 33, 
the method of solution is to eliminate n' unknowns between 
the " Observation " and " Condition " equations, leaving 
q — n' independent unknowns in the " Normal " equations. 
Consequently formula (77) now applies to this case and it 
would be written 



. = v/' 



'Zpv^ 



n — q -\- n' 



(79) 



also ,.. = v/ . ^^""1 '. m 

^ Pz(n — q -\- n) ^ ' 

The weights of the q — n' unknown quantities can be 
found by any one of the three methods already given and 
then the mean error of each unknown may be computed by 
using formula (80). If the mean errors of the n' quantities 
that were first eliminated is wanted, their weights must be 
determined by eliminating a different set of n' quantities 
from the original observation equations and solving the neces- 
sary sets of equations for these weights. The first method of 
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solution for the weights would perhaps be best here as the 
actual values of the unknowns have already been found. 
87. Example. Given the 

OBSERVATION EQUATIONS 



(A) 



Zl 


4- 


2j - 3.0 = 




weight 1 


2l 




- S3 + 1.5 = 
S2 - 2.2 = 




« 4 
3 


Si 




+ S3 _ 8.4 = 




2 


and the 














CONDITION EQUATION 








»8 — «!! — 0.5 = 








(B) 

To find the most probable values of Sj, Sj, and S3, and 
also their mean errors. 

Eliminating Z3 between (A) and (B) there remains 



gj 4- s, _ 3.0 


= 


weight 1 


Si - sj + 1.0 


= 


« 4 


s^ - 2.2 


= 


" 3 


Si + 22 - 2.9 


= 


« 2 



From these we have the 



NORMAL EQUATIONS 

7 Si - sj - 4.8 = 
_ 2i -f 10 Za - 19.4 = 



(C) 



Solving, 



sj = 0.98 p^, = 6.9 

s^ = 2.04 , p^ = 9.9 (D) 

from (B) 33 = 2.54 



(E) 
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Now eliminating Zj between (A) and (B), we find 

2i + 28 — 3.5 = weight 1 

z^ — Zs ^ 1.4 =0 « 4 

23 - 2.7 = "3 

gj _j_ gj _ 3.4 = « 2 

From these we derive the new set of 

NORMAL EQUATIONS 

7 Si - Z3 — 4.3 = 
- Zi + 10 Za — 24.4 = 

and solving for z, 

Za = 2.54 jo,3 = 9.9 (G) 

Substituting the values of z^, Z2 and 23 in equations (A), 
we have the residuals 



(F) 



V 


v^ 


pv^ 


.02 


.0004 


.0004 


.06 


.0036 


.0144 


.16 


.0256 


.0768 


.12 


.0144 


.0288 



.1204 = %pv^ 
Since in this example n = 4, w' = 1, q = 3, 



By (79) 



_ ./ .1204 _ 



/*2l = 



M2i. = 



\/9.9 



/.1204 
' 2 


= .25 


v/6.9 


= .09 


/* 


= .08 


^9.9 


/^ 


= .08 
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The first significant figure in the mean errors is so large 
that it is not worth while to retain the second place as usual. 

If desired, the probable errors and average deviations can 
now be computed by the usual formulas. 

88. In case the observations are made directly upon the 
values of several quantities subject to certain conditions, we 
have n ^ q, and equation (79) reduces to 



^ = v/^" (s» 



n 



from which the mean error of any observation may at once 
be computed from its weight. 

89. Example. Taking the example in paragraph 34 on 
the measurement of the angles of a quadrilateral, we had for 



OBSERVATION EQUATIONS 

3i = weight 3 

s, = "2 

z, = "2 

Z4 = "1 



(A) 



for the 



and for the 



CONDITION EQUATION 

-f s, + s, + 2, + 58 = (B) 

NORMAL EQUATIONS 

4 si + 82 + 2s + 58 = 
2i -j- 3 z, -I- z, + 58 = (C) 

gj -I- Zj -i- 8 Zs + 58 = 

Solving these equations for the values of the unknown 
quantities and also for their weights, we have 

Zi = — 8.3 p^, = 3.5 

z, = - 12.4 p^ = 2.5 

s, = - 12.4 p^ = 2.5 
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and forming a new set of normals containing e^, we find, on 
solving for that quantity, 

a^ = — 24.9 jOj. = 1.75 

These values of s,, z^, Zg, Zi are also in this case the 
residuals of the observations, and therefore to compute the 
precision measures we have 



V 


v' 


pv' 


8.3 


68.9 


206.7 


12.4 


153.8 


307.6 


12.4 


153.8 


307.6 


24.9 


620.0 


620.0 



= 1 



1441.9 = 2i»»' 



By (81) ^ = y/liM:i = 38 

/**.= -=:= 20 r,. = 13 

V'3.5 

/»*• = /**.= -=: =24 r^, = r^ = 16 

V2.5 

Mz. = -== = 29 r„ = 20 

V'1.75 

We shall accordingly write for the most probable values of 
the angles of the quadrilateral 



A 


= 101° 


13' 


14" ± 


13" 


B 


= 93 


49 


5 ± 


16 


C 


= 87 


5 


27 ± 


16 


B 


= 77 


52 


15 ± 


20 



In the original solution in paragraph 34 it is evident that 
the results were carried out to a greater number of places of 
significant figures than the character of the observations 
warranted. 



CHAPTER V. 

MISCELLANEOUS THEOREMS. 

THE DISTRIBUTION OF ERRORS. 

90. Having developed the processes for the adjustment of 
observations according to the Method of Least Squares, it 
will now be interesting to show how closely the distribution 
of errors found in actual practice corresponds to the theo- 
retical distribution upon which our methods of solution are 
based. 

By formula (36), the probability that the error of a single 
observation will be numerically less than a is 



P = M C\-hH-^dx (82) 

\l TT Jo 



Let t ■= hx, .•. dt z= h dx. Also when a; = 0, 

« = ; and when x z= a, t z= ha =:=: p-. Sub- 

r 
stituting in (82), 

i> = 4= r'^'^dt (83; 

Values of P for values of the argument - are given in 

Table I. Also for any series of observations this quantity 
P will represent the fraction of the entire number which 
should have errors less than the amount a. Hence if P is 
multiplied by the whole number of observations the result 
will be the number of errors which should be less than the 
limit a. 



84 
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91. Example. Forty measurements on the diameter of 
Saturn's ring were made by Bessel, with the following 
results : — 



M V 


M 


V 


M 


V 


M V 


38".91 - .40 


39",35 


+ .04 


39".41 


+ .10 


39".02 -.29 


39 .32 +.01 


39 .25 


— .06 


39 .40 


+ .09 


39 .01 - .30 


38 .93 - .38 


39 .14 


-.17 


39 .36 


+ .05 


38 .86 - .45 


39 .31 .00 


39 .47 


+ .16 


39 .20 


-.11 


39 .51 +.20 


39 .17 — .14 


39 .29 


-.02 


39 .42 


+ .11 


39 .21 - .10 


39 .04 — .27 


39 .32 


+ .01 


39 .30 


-.01 


39 .17 - .14 


39 .57 +.26 


39 .40 


+ .09 


39 .41 


+ .10 


39 .60 +.29 


39 .46 +.15 


39 .33 


+ .02 


39 .43 


+ .12 


39 .54 +-23 


39 .30 - .01 


39 .28 


-.03 


39 .43 


+ .12 


39 .45 + .14 


39 .03 - .28 


39 .62 


+ .31 


39 .36 


+ .05 


39 .72 +.41 



From these the most probable value of the diameter is 
found to be 

J) = 39".308 ± 0".022 

the probable error of a single observation being r = 0".136. 
Compare the theoretical and actual distribution of errors 

between 0".00 and 0".05 



over 



.05 


(C 


.10 


.10 


(C 


.20 


.20 


it 


.30 


.30 


u 


.40 


.40 







In the following table the first column gives the successive 
values of the limiting error a, the second column the values 

a 
of -, and the third column the corresponding values of I*- 
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The fourth column contains the differences between the suc- 
cessive values of P, and by multiplying each of these 
differences by 40, the number of observations, we have the 
quantities in column five, which are the numbers of errors 
that according to the theory should fall within the corre- 
sponding limits. Column six shows the actual number of 
residuals occurring between these limits. 



a 


a 
r 


P 


d 


n 


n' 


0.00 


0.000 


0.000 


0.196 

0.184 


8 

7 


9 
6 


0.05 


0.368 


0.196 


0.10 


0.735 


0.380 


0.299 


12 


12 


0.20 


1.471 


0.679 


0.184 


7 


8 


0.30 


2.206 


0.863 


0.090 


4 


3 


0.40 


2.942 


0.953 


0.047 


2 


2 


OO 


00 


1.000 









This is a close agreement considering that the number of 
observations is not very large. Also the number of errors 
greater than the probable error should be equal to the num- 
ber less than it. On counting the residuals we find twenty-one 
less than 0".136 and nineteen greater. 



THE REJECTION OF OBSERVATIONS. 

92. After a series of measurements have been made, it is 
frequently found that one or two of the observations differ 
widely from the others, and hence it becomes a matter of 
great importance to establish, if possible, some criterion by 
which we may determine whether such discordant observa- 
tions should be rejected or not. We are not concerned here 
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with the question of the detection of a mistake or constant 
error, which a consideration of the circumstances of the 
observations or of the instruments might reveal, but it is 
assumed that there is nothing whatever to guide us except 
the naere fact of the unusual size of the residuals of the 
observations under discussion. To reject an observation 
merely because it differs considerably from the others is 
entirely unjustifiable, while to retain it without any investi- 
gation is a neglect of the evidence furnished by the observa- 
tions themselves. 

The adoption of any rigid criterion based upon the magni- 
tude of the residuals is perhaps more satisfactory from a 
mathematical standpoint than from that of a practical observer, 
and some of the latter are of the opinion that no observation 
should be rejected entirely, even the most widely discordant 
ones being given a certain weight. In this latter case, 
however, the Theory of Probability will furnish a guide as to 
the proper weights to assign to the different observations. 

Of the various criteria that have been proposed, that 
developed by Pierce (see Chauvenet, page 558) is perhaps 
the most complete. The derivation and application of this 
criterion is, however, somewhat long and complicated, and 
for all ordinary cases the following simple methods will give 
practically as good results. 

93. Criterion for the Rejection of a Single Donbtfnl 
Observation. It was shown in (83) that in a series of n 
observations the number of errors numerically less than a 
should be nP, and therefore the number of errors greater 
than a should be 

n — nP = n(l — P) (84) 

If the value of the expression in (84) is less than one-half, 
the occurrence of an error of magnitude a will have a 
greater probability against it than for it, and hence the obser- 
vation corresponding may be rejected. 
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Accordingly the limit of rejection, a, of a single doubtful 
observation is obtained from the equation 



n (1 - P) = _ 
2 



or 



2 n 



2n 



(85) 



94. Example. Fifteen observations on the value of an 
angle are made. Ought any of the observations to be 
rejected ? 



M 


V 


w^ 


v' 


w' = 


v" 


«"2 


2° 23'.90 





.30 


.090 


- .41 


.168 


_ 


.33 


.109 


23.76 


— 


.44 


.194 


- .55 


.303 





.47 


.221 


25.21 


+ 


1.01 


1.020 


+ .90 


.810 








24.68 


+ 


.48 


.230 


+ .37 


.137 


+ 


.45 


.203 


23.96 


— 


.24 


.058 


- .35 


.123 




.27 


.073 


24.26 


+ 


.06 


.004 


-.05 


.003 


+ 


.03 


.001 


24.83 


+ 


.63 


.397 


+ .52 


.270 


+ 


.60 


.360 


24.07 


— 


.13 


.017 


-.24 


.058 




.16 


.026 


23.98 


— 


.22 


.048 


- .33 


.109 





.25 


.063 


24.14 


— 


.06 


.004 


-.17 


.029 





.09 


.008 


24.40 


+ 


.20 


.040 


+ .09 


.008 


+ 


.17 


.029 


24.38 


+ 


.18 


.032 


+ .07 


.005 


+ 


.15 


.023 


24.59 


+ 


.39 


.152 


-.28 


.078 


+ 


.36 


.130 


24.10 




.10 


.010 


+ .21 


.044 




.13 


.017 


22.80 


— 


1.40 


1.960 












2° 24 .20 






4.256 




2.145 






1.263 


Jtfo 






•Zv'' 




21)'^ 






■Zv"^ 
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Using all the observations we find 



Mo = 2° 24'.20 r = .6745 v/l^ = .37 

^ 14 

29 
By (85), P = _ = .967 

By Table I, - = 3.17 .-. a = 1.17 



As the residual — 1.40 is larger than a, we reject the 
last observation. 

From the remaining observations we now compute a new 
mean value and a new set of residuals. And we find 



/ 2.145 
M'o = 2° 24'.31 V = .6745 y = 


.27 


27 
By (85), P = — = .964 




By Table I, - = 3.11 .-. a = .84 





The third observation may accordingly be rejected. 
From the thirteen observations that remain we find 



/l 2fi^ 

M"o = 2° 24.-i;J r" = .6745 \/^—^ = 

V 10 



22 
12 



P = — = .962 
26 



- = 3.08 .-. a = .68 

r 



Therefore no more observations are to be rejected. 
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95. The Huge Error. In cases where the number of 
observations is not unusually large, a simple and safe criterion 
for the rejection of a doubtful observation is found in the use 
of the " Huge Error." 

This is an error of such a magnitude that 999 out of every 
1000 errors are less than it and only 1 as large as or greater 
than it. 

Therefore the probability that the error of any given 
observation will be less than the " Huge Error " is .999, and 
from Table I, 

when P = .999, — = 4.9 

r 

a = Huge Error = 4.9 *• 

= 3.3 |i (86) 

= 4.1 a.d. 

Then in any limited series of observations, if an error 
greater than the huge error is found, we should reject the 
observation corresponding. 

See also, Solman, page 30 ; Wright, page 131. 



CONSTANT ERRORS. 

96. Throughout our discussion of the methods of adjusting 
observations so as to obtain from them the most probable 
values of the unknown quantities, all constant errors are 
supposed to have been eliminated before the Method of 
Least Squares is applied in deducing the results. 

If this is not done, and each observation is subject to the 
same constant error, the final result will be affected by an 
equal amount, and in short, the Method of Least Squares is 
not capable of removing or reducing the effect of errors of 
this kind. All that is accomplished by the use of the method 
is to reduce to a minimum the effect of the Accidental Errors. 
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Hence it will be seen that although by increasing the 
number of observations of a given kind we may increase the 
precision, that is, reduce the probable error, of our final 
result as much as we choose, yet we do not in this way 
necessarily increase the accuracy of the determination. 

But if the unknowns can be determined in several ways, or 
under a variety of different circumstances, with various 
instruments, or by different observers, then it is most proba- 
ble that the constant errors of the different sets of measure- 
ments will be grouped about the true values of the unknowns 
according to the exponential law of error. Accordingly a 
combination of such observations will enable us to increase 
not only the precision, but also the accuracy of the final 
result, the constant errors of the different sets tending to 
cancel each other in the same way that the accidental errors 
of a single set do. 

It is for this reason that determinations of a quantity from 
observations made in a variety of ways are more valuable 
than those obtained merely from different sets of measure- 
ments of the same kind. 

97. The probability of the existence of a constant error 
may often be expressed in the following manner. 

Mxample. A standard 100 ohm coil is compared with 
a Wheatstone's bridge and the mean result found to be 
100.90 ± 0.20. To find the probability that there is an error 
in the bridge between -)- 0.30 and -|- 1.50 ohms. 

Suppose the result 100.90 ± 0.20 is treated as a single 
observation, and we find by an application of (83) the prob- 
ability that the error of this observation is numerically less 
than 0.60 ohms. 

Here a _ M ^ ^^^ ^ ^ ^^^ 

»• .20 

Hence, as far as is shown by the observations, the proba- 
bility that 100.90 ohms is within 0.60 ohms of the true 
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value is .957. But since it is known that the true resistance 
is 1 00 ohms, it follows that there is the same probability 
that there is a constant error in the bridge between -|- 0.30 
and -\- 1.50 ohms. 

98. Combination of Determinations having Different Con- 
stant Errors. In case two or more determinations of a 
quantity, together with their probable errors, are obtained, 
the method of combining them so as to secure the best final 
result was considered in paragraph 56, and in Example C, 
paragraph 57. But it was there assumed that all the results 
were subject to the same constant errors, while if this is not 
true the probable errors of the separate determinations bear 
no relation to their weights, and accordingly in such cases 
another process must be adopted. 

To determine whether the different measurements may 
fairly be considered to have the same constant errors we may 
proceed as follows : — 

Let the determinations of the quantity M be 

Jfi ± ri (a) 

and let the difference between these results be 

<Z = JIfi - Jf, (c) 

Then the probable error of d is by (66), 



R = s/r,' + r,' (d) 

If d is of such a magnitude that an accidental error as 
great as it may reasonably be expected, we may assume that 
the constant errors of (a) and (b) are the same, and pro- 
ceed as in Example C, paragraph 57. 

But if the probability of making two determinations which 
differ by the amount d is very small, we had best consider 
Ml and M^ to have the same weight, provided there are no 
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special reasons for regarding one better than the other. The 
final value ot M will then be the arithmetical mean of Mi 
and il/o, and its probable error will be found by (53). 

99. Exam]jle A. An angle is measured by a theodolite 
and by a transit with results 

By Theodolite, 24° 13' 36".0 ± 3".l 
By Transit, 24° 13' 24" ± 14" 

What is the most probable value of the angle and its prob- 
able error ? 

Referring to the preceding paragraph, 

r-^ = 8.1, rj =14, c? = 12, 

and the probable error of d is 



a = ^'i.V -\- 14^ = 14 

Then from Table I the probability that the accidental 
error of a determination will be at least as large as 12 is 
found from 

^ = ^ = .86 ... 1 - P = .57 

r 14 

That is, there is more than an even chance that two such 
determinations of the angle will differ by as much as 12. 

Hence it is fair to assume that the two determinations are 
not affected by constant errors of different magnitudes, and 
they would be combined as in Example C, paragraph 57. 

Example S. Suppose the zenith distance, M, of a star, 
observed at two different culminations, is found to be 

Ml = 14° 53' 12".10 ± 0".30 
3T^ = 14° 53' 14".30 ± 0".50 

What is the best final value ? 

Here d = 2.2, H = \/ .09 -|- .25 = .58 

a 2.2 

and for - = — = 3.8, 1 — P = .01 
r .58 
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Therefore the chance that the difference in the two deter- 
minations, due to accidental errors, will be as large as 2.2 is 
only one in a hundred. It is to be concluded then that the 
constant errors of observations at the two culminations differ 
by about 2.2, and as there is nothing to show that one meas- 
urement is more accurate than the other we will give them 
both the same weight and take the mean. Then the best 
value for the zenith distance is 

Mo — 14° 53' 13".20 ± 0".74 

12".10 4- 14".30 
For Jfo = 14° 53' + ^ 

= 14° 53' 13".20 



/ 2.42 
By (53) r„ = -S^^^y^^^ = .74 

For a more extended treatment of this subject see Johnson, 
" The Theory of Errors and Method of Least Squares," 
chap. vii. 

THE WEIGHTING OF OBSERVATIONS, 

100. In case the relative worth of observations is not 
settled by methods already discussed, the proper weight to 
assign to each quantity in the final adjustment can only be 
determined from a full knowledge of all the circumstances of 
the measurements. Even then considerable experience in 
the particular work in hand is required before the best values 
for these weights can be assigned. The weight given to a 
quantity should never be considered final, but always subject 
to revision whenever new information with regard to the 
quantity is obtained. Thus an observation which at first is 
supposed to deserve a high degree of confidence is often 
found on later investigation to possess very little value, and 
vice versa. 

See also Wright, page 118. 
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OTHER LAWS OF ERROR. 

101. Although in the great majority of cases the distribu- 
tion of errors follows the exponential law thus far considered, 
there are a few special cases in which some of the suppositions 
made in deriving that law do not hold, and hence for the 
adjustment of such observations the corresponding special 
laws of error must be determined. 

For instance, in applying the exponential law we assume a 
large number of observations, that each observation is sub- 
ject to the same law of error, that small errors are more 
likely to occur than large ones, and that positive and negative 
errors are equally probable. Now it is easy to conceive of 
cases where only positive errors can occur, or where the 
probability of the occurrence of a small error may not be 
greater than that of a larger one, etc. If we can determine 
the different sources of error in any case and the relative 
effect of each upon the quantity sought, we shall arrive at 
the law of error for that particular set of observations. The 
case of most common occurrence is the following. 

102. Suppose all errors between the limits a and — a 
are equally probable, and that there are no errors beyond 
these limits. Then \i y -^ <^(*) is the equation of the 
Curve of Error, and its area is represented as in para- 
graph 18, we have 

P<^(a!) dx = 1 (a) 

2</)(a!) C"dx = 1 (b) 

•/ 



or 



since by the supposition made <^ (x) must be a constant. 
Integrating and solving for <^(x), we have 



y = <t.(a5) = -L (-87) 

2 a 
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To find the Mean Error we have by definition as in para- 
graph 52 

ft" = f" a;' </> (a;) dx 

^ -a 



__ 1 /'o 

a Jo 



cb" dx 



a" 



K- = 



a 



The Probable Error is derived from the equation 

<f>(x) dx = 

r 2 

A Cdx = 1 

a Jo 2 



r = 



a 



Finally, for the Average Deviation we have 
a.d. := I X <t> (x) dx 

J ~a 

= _ I X dx 

a Jo 



a.d. = — 



And the Curve of Error has the form 



(88) 



(89) 



(90) 
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That the average deviation and probable error are in this 
case equal to one half of a may also be seen from the defini- 
tion of these quantities. 

Example. In taking a logarithm from a four place table, 
what is the probable error of the mantissa ? 

In this case the maximum error is .00005, and all 
errors between — .00005 and .00005 are equally probable. 
Therefore 

r = .000025 

103. The only other special case of common occurrence is 
that in which the error of a quantity is due to two sources, 
each of which can with the same probability assume all 
values between a and — a. Here it may be shown that 
the curve of error consists of two straight lines whose equa- 
tions are 

y = and y =z ■ (91) 

Also |i» = - a% r = (2-\/J)a (92) 

For a more extended discussion of special laws of error, see 
Wright, paragraphs 31 to 39. See also Example 151. 

104. Contradictory ObserTations. Suppose three observa- 
tions of a quantity give results 55, 56, 91. It is obvious that 
to take the arithmetical mean as the most probable value in 
this case would be contrary to the evidence furnished by the 
measurements, while in so small a number of observations it 
would not be allowable to reject the value 91 entirely. With 
such a series of observations no satisfactory solution can be 
obtained, but probably the best thing to do would be to take 
the value which has as many observations less than it as it 
has greater, or 56. 



CHAPTEK VI. 

GAUSS'S METHOD OP SUBSTITUTION. 

105. The most laborious part of the application of the 
Method of Least Squares to the adjustment of observations 
consists in the formation and solution of the "Normal Equa- 
tions," and this labor increases enormously with increase in 
the number of observations, of unknowns and of conditions. 
It is not at all unusual to find that the adjustment of a single 
set of observations takes several weeks, even with all the aid 
that can be obtained from tables of logarithms, of squares, of 
products, and of reciprocals, and also from the use of calcu- 
lating machines. In general the computations can be per- 
formed more rapidly and with less fatigue by using a machine 
or a table of products than by using logarithms, but a combi- 
nation of methods is often desirable. 

106. In any case however it is of the utmost importance 
that the formation and solution of the normal equations 
should be effected in a systematic way, and that as far as 
possible checks on the numerical work be carried along in 
the computations. By a slight amount of additional work 
checks upon the results at successive stages of the solution 
may be obtained by means indicated in the following dem- 
onstrations, while the most satisfactory form for the solution 
of the normal equations is given by the "Method of Substitu- 
tion" proposed by Gauss. 

This method will now be explained ; and it is to be observed 
that it is customary in this subject to enclose a quantity in 
brackets when the sum of a number of quantities of the same 
kind is to be denoted. Thus [a6] means the same as 2a6. 

For the sake of simplicity in demonstration it will be 
assumed that the observations are all reduced to weight unity. 
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107. Checks on the Formation of the Normal Equations. 
If, as in paragraph 76, we take for 

OBSERVATION EQUATIONS 
OiSi + Ji22 + . . . qiZq 4- Wli = Ui 



«22l + *2Ss + 



. . q^Zq + m^ = Vi 



(A) 



a,fii + K^i -\- • • ■ qrfiq -\- rrin = v, 
we shall have for 

NORMAL EQUATIONS 

[aa'] Zi -1- [a6] Zj -|- . . . [ag'] z, -|- [am'} 



[a&] z, + [W] z, + 



[a? J «! + [*2] 32 + • 



Ibq} Zq + [Swi] = 








(B) 



Let 






qi + mi = Si 
q^ + m^ = Sj 



(C) 



a„ +*»+... J'n + w„ = s„ 



... [a] + [J] + . . . [y] 4- [m] = [s] (D) 

Multiplying the first of (C) by wij, the second by w,, ... 
and adding, there results 

lam] + [&mj + • • . [<?m] + [mm] = [sm] (93) 

Next multiplying each of equations (C^ by its a and 
adding, and then each by its b and adding, etc., we have 
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[aa] + [CT&] + . . . [ag] + [awi] = [as] 
[a6J + [&&] + . . . [bq] + [bm] = [&s] 

(94) 

[aq] + [&(?] + . . . [gg] + [gm] = [gs] 

Equation (93) will be satisfied if the absolute terms in 
the normal equations are correct, and equations (94) when 
the coefiicients of the unknown quantities are correct. These 
check the formation of the normal equations. 

108. The Reduced Normal Equations and the Elimination 
Equations. 

The value of gj in terms of the remaining unknowns, 
derived from the first of equations (B), is 

Substituting this in the remaining n — 1 equations, they 
become 



= 






[5?] - [«?] ^1 .. + . . . j[?-] - [a?] ^1 = 

And letting 

[abl Vabl 

[5^ 1] = m - ^-^ 

rao") rami 
[gm, 1] = [qm2 - ^-^f^ 
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the above equations take the following form, which, being 
the same as that of the original normal equations, they are 
called the 

FIRST REDUCED NORMAL EQUATIONS. 

[&6,l]s, + [6c,l]23+...[ft?,l]sg + [5»«,ll = 
[be, 1] zj + \_ce, 1] Za + . . . {cq, 1] z, + [cwi, 1] = 

(G) 

\bq, 1] 2^ + Icq, 1] Zs + . . . ^qq, 1] z, + [g-m, 1] = 

An inspection of equations (F) will render it easy to form 
a rule for writing out any one of them. 

Now by means of the first of equations (G), eliminating 
Zj from each of the others in the same way that z^ was 
eliminated from the normal equations, there results the 

SECOND REDUCED NORMAL EQUATIONS 

[cc, 2] Zj + . . . \_cq, 2] Zg + [c»i, 2] = 



(H) 



Icq, 2] Za + . . . [qq, 2] z, + ^qm, 2] = 

In which 

[5c, 1] \bc, 1] 



[cc,2] = [cc,l] - 



\bb, 1] 



(I) 



[^m,2] = [^«^,l] - ^-^^^^^M 



Continuing this process we shall finally arrive at the single 
equation 

Iqq, q-1] z, + {qm, q-1] = (J) 

from which the value of Zg is determined. 
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The value of s^.j will then be obtained by substituting the 
numerical value of s^ in the first of the preceding set of 
equations, and so on, until finally Zj is obtained from the 
first of the original normal equations. The equations from 
which the unknowns are actually determined are then the 
following, called the 

ELIMINATION EQUATIONS. 

laa] «i + lab] Sj -f . . . [aq] », + [am] = 
Ibb, !]«,+... Ibq, 1] «g + [6m, 1] = 

(95) 

lqq,q-l]iSg + [qm, q-\] = 

It may be seen from the rule in paragraph 81 that 
[??) •1~1] i^ *^^ weight of Zq, and the weight of any 
unknown might be found at the same time as its value by 
making it the last in the order of elimination, but except in 
special cases the weights had best be obtained by the general 
process of paragraph 115. 

109. Check on the Solution of the Normal Equations. 
Multiplying the first of the Observation Equations (A) by 
nil, the second by m^, . . . and adding the results, we have 

[wjw] = [am] Zi -J- [ftwi] 32 + • • • [^^] 2g + ["inm] 

But in equation (T), paragraph 82, it was shown that 
[mv] = [w]- Therefore 

[vv] = [am] Zi -j- [bm] z^ -\- . . . [qm] z, + [mm] 

Substituting in this the value of Zj from the first of (95), 
we get the result 

[vv] = [bm, 1] Z2 + [cm, 1] Zg -|- . . . [qm, 1] Zg -f [mm, 1] 

in which 
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[aJ] [aw] 



[5m, 1] = [5m] — 



£mm, 1] = [mm] — 



[aa] 



[aa] 



(K) 



being similar in form to equations (F). 

Next eliminating 2j in a like manner, we get 

[yv'] = \_cm, 2] Zj -|- . . . [g-m, 2] 8, -j- [mm, 2] 

and continuing this process it finally appears that 

[vv'] = [rnm, q] [96] 

110. Arrangement of the Compntations. In computing 
the coefficients that appear in the "Auxiliary" or "Reduced 
Normal Equations" it is most convenient to arrange the work 
in tabular form. The arrangement of the solution will be 
illustrated for an example containing four unknowns but it 
will be evident that the process can be extended to cover any 
case. 

Let the Observation and Normal Equations be represented 
by equations (A) and (B), there being only four unknowns 
s,, Sj, 2a, S4, and arrange a table as on the next page. In this 
scheme of solution the upper lines of the rows in the first 
compartment contain all the quantities that appear in the 
Normal Equations, together with [mm] and the quantities 
in the column headed s which are used in checking the 
results in accordance with equations (93) and (94). The 
other compartments contain the corresponding quantities 
for the Reduced Normal Equations, and the first line in 
each compartment gives the coefficients in the Elimination 
Equations. 
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Scheme A. — Solution of the Normal Equations. 





a 


i 


t 


rf 


m 


J 


a 


[oa] 
Iog[<«.] 


[oJ] 
log [aJ] 


log[«] 


[«rf] 
log [nrf] 


log [om] 


[«] 
log [«] ■ 


b 


log ^6 
log^. 

log^„ 


[M] 

^»[<-i] 

Iog^»[ai] 


[J.] 

l0g^.[ac] 


[M] 
log.<.[<.rf] 


[*m] 

Ai [ami 

log Jj [am] 


[4,J . 
log-4j[a.] 


e 




log^cM 


[erf] 
logA[arf] 


[em] 

J, [am] 

log ^g [am] 


[«] 
log ^,, [as] 


d 




[rfrf] 

^„[<.rf] 

log^j[orf] 


[rfm] 

Ai [am] 

log ^rf [am] 


[A] 
log ^,, [as] 


m 




Am [am] 
log J„ [ami 


[ms] 
log -4„ [as] 


h 


logB. 
log5„ 
\ogB. 


[ii. 1] 
log[iJ.l] 


log [4c, 1] 


[M,l] 
log[M.l] 


[*m,l] 
log [Am, 1 ] 


[is, 1] 
log [4s, I] 


e 




[«, 1] 

B.[ie.l] 

log B, [id] 


[crf,l] 

B.[M,1] 

logS,[irf,l] 


[em, 1] 

B.[4m,l] 

log5.[im, 1] 


[e., 1] 

5.[4«. 1] 

logBj[4s, 1] 


d 




[rfrf, I] 

£d[M,I] 

logfl„[4rf,l] 


,[rfm,l] , 
£„[47n,l] 
log£i[Jm,l] 


[rfs, 1] 

Bi [4., 1] 

log Bi [4s, I] 


m 




[mm, 1] 

£«[*m, 1] 

logB„[Jm,I] 


[m., 1] 

-»™ [4s. I] 

logB„[4s.l] 


e 


logCi 
logC. 


log[«,2] 


[erf. 2] 
log [erf, 2] 


[™,2] 
log [em, 2] 


[«, 2] 
log [es, 2] 


d 




[rfrf, 2] 

C,,[crf.2] 

log Cj [erf, 2] 


[rfm,2] 

Ci[em,2] 

log Ci [cm, 2] 


[rfs, 2] 

Ci \c,, 2] 

logC„[es,2] 


m 




\mm, 2] 

C„[em,2] 

IogC„[em,2] 


[m,,2] 

C„[«,2] 

log C„[cs. 2] 


d 


logiJ. = 


log - «4 




[rfrf, 3] 
log [rfrf, 3] 


[rfm,3] 
log [rfm, 3] 


[rf«, 3] 
log [rfs, 3] 


ffl 




\mm, 3] 

.0„[rfm,3] 

log/)„[rfm,3] 


[ms, 3] 

-Om [rfs, 3] 

logZI„[rfs,3] 


m 




'i - 




M = 


[mm, 4] 


im>, 4] 
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The logarithms of the quantities in the first row of each 
compartment are also written in, and from these by proper 
subtractions are obtained the logarithms in the margin, where 

A _ [«^] A _ [«c] . _ jam-] 

* ~ [aa]' ' — [aaj " ' • ^"^ — [««] 



_ Ibc, 1] _ [bm, 1] 

" ~ Ibb, 1]' • • • -""• — [65, 1] 



(L) 

_ s - i> - '^'^^'^^ 
* ■" *" ~" [dd, 3] 

Now in each compartment adding the logarithms at the 
margin to each of the logarithms in the first row of that com- 
partment we obtain the corresponding logarithms written in 
the other rows. The numbers represented by these logarithms 
are next written above them, and if each of these quantities 
is then subtracted from the one above it the result will be 
the corresponding quantity in the compartment below. Some 
of the squares in each compartment are left vacant as the 
quantities belonging to them have already appeared above. 

111. Application of Checks. Also, by (93) and (94), in 
the first compartment the quantities in the first lines of the 
last column should be equal to the sum of all the quantities 
in the first lines of the corresponding rows plus the quantities 
similarly situated above the first terms of the rows. Similar 
checks will apply in each compartment; for if from the 
second of equations (94) we subtract the product of the first 
equation multiplied by Ab we have 

Ibb, 1] + Ibc, 1] + . . . [bm, 1] = lbs, 1] (97) 

In the same manner we may show that a corresponding 
check holds throughout, so that finally we shall have 

[mm, 4] = lm,s, 4] (98) 
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The last compartment of the table is added to give this 
final check and the value of [uy] in accordance with equa- 
tion (96). 

If the multiplications and divisions are simple or if a table 
of squares or products or a computing machine is used the 
logarithms will of course be omitted from the scheme of 
solution. 

112. Example. In order to illustrate the systematic for- 
mation of the coefficients that appear in the Normal Equations 
as well as the solution of the latter by the above method we 
will take the 



OBSERVATION EQUATIONS 

3l - Zj 



«i 4- 22 + Zs + a* = 0.1 
Sl + 22 - Ss — Z4 = 0.6 

Si — 2zj — ^^-\-z^ = 0.1 2; 



- 234 = 

- z, 4- z^ = 

— 3» = 



0.3 

— 0.1 

0.4 



First form a table containing the coefficients in these equa- 
tions and also the sums s. As a first check the sum of the 
quantities in column s should be equal to the sum of all the 
quantities in all the other columns. 



I. Coefficients ix the Observation Equations. 



No. 


a 


b 


c 


d 


m 


s 


1 


- 1 


1 


1 


1 


- .1 


1.9 


2 


1 


1 


- 1 


- 1 


- .6 


— .6 


3 


1 


— 2 


- 1 


1 


— .1 


- 1.1 


4 


1 


— 1 





- 2 


- .3 


- 2.3 


5 





1 


- 1 


1 


.1 


1.1 


6 


1 





- 1 





- .4 


- .4 


Sum. 


3 





- 3 





- 1.4 


— 1.4 
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From these we now compute the coeflScients in the Normal 
Equations (B), and also the necessary quantities for the 
check equations (93) and (94). 

II. Coefficients in the Normal Equations. 



N'o. 


aa 


ab 


ac 


ad 


am 


as 


bb 


be 


bd 


bm 


1 


1 


-1 


-1 


-1 


.1 


-1.9 


1 


1 


1 


-.1 


2 


1 


1 


-1 


-1 


-.6 


- .6 


1 


-1 


-1 


-.6 


3 


1 


-2 


-1 


1 


-.1 


-1.1 


4 


2 


-2 


.2 


4 


1 


-1 





-2 


-.3 


-2.3 


1 





2 


.3 


5 














.0 


.0 


1 


-1 


1 


.1 


6 


1 





-1 





-.4 


- .4 









1 


.0 


Sum 


5 


-3 


-4 


-3 


-1.3 


-6.3 


8 


1 


-.1 



N'o. 


bs 


cc 


cd 


cm 


cs 


dd 


dm 


ds 


771771 

.01 


ms 


1 


1.9 


1 


1 


-.1 


1.9 




-.1 


1.9 


-.19 


2 


- .6 


1 


1 


.6 


.6 




.6 


.6 


.36 


.36 


3 


2.2 


1 


-1 


.1 


1.1 




-.1 


-1.1 


.01 


.11 


4 


2.3 








.0 


.0 




.6 


4.6 


.09 


.69 


5 


1.1 


1 


-1 


-.1 


-1.1 




.1 


1.1 


.01 


.11 


6 


.0 


1 





.4 


.4 





.0 


.0 


.16 
.64 


.16 


Sum 


6.9 


5 





.9 


2.9 


8 


1.1 


7.1 


1.24 



If the coefficients in the Observation Equations are large, 
so that it becomes convenient to use logarithms, other tables 
corresponding to I and II would be formed to contain these 
logarithms. 
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Substituting these quantities now in the general tabular 
scheme of paragraph 110 we have the results on the following 
page. The work of the first compartment is performed 
without the use of logarithms, as the numbers are simple. 
The quantities in the last column should all be zero accord- 
ing to the check equations, what small differences there are 
being due to the rejection of figures beyond the third place 
in the decimals. The decimal points in logarithms to which 
correspond negative numbers have been replaced by the 
letter n. 

The demonstrations that have been made now enable us to 
see at once from an inspection of the results in this table 
that 

S4 = 0.238 

Pz. = 1-6 
[to] = .007 

Therefore substituting in equations (73) and (74) we 
have 

;.,. = .047 r,. = .032 

If the two values of [«v] obtained in the solution had 
differed at all we should have taken the mean of the two. 

113. If the Elimination Equations (95) are divided by 
l^aa^, [W, 1], [cc, 2], [<?t?, 3], respectively, they become 



»i -\- Ab»2 -\- Ac «3 -\- Ad»i -\- Am = 

»z + CdSi + Cm, = 



(99) 



And the solution for the unknowns can be effected most 
conveniently by arranging the computations in the manner 
illustrated on page 109. 
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Scheme A. Solution op the Normal Equation's. 
a b c d m 8 S 



A, 



d 

Aa 
m 



B. 



m 



B„ 



d 
m 



t'm 



d 



m 



B„ 



5 


-3 


-4 


-3 


-1.3 


-6.3 






-.6 


8 
1.8 


1 
2.4 


1 

1.8 


- .1 

.78 


6.9 

3.78 


-.8 




5 
3.2 




2.4 


.9 
1.04 


2.9 
5.04 





-.6 




8 
1.8 


1.1 

.78 


7.1 

3.78 





-.26 




.64 
.34 


1.24 
1.64 







6.2 
0.7924 


-1.4 
0„1461 


- .8 
9„9031 


— 88 
9„9445 


3.12 
0.4924 





9„3537 




1.8 
.316 

9.4998 


-2.4 

.181 
9.2568 


- .14 
.199 

9.2982 


-2.14 
- .705 
9„8479 





9„1107 




6.2 
.103 
9.0138 


.32 
.114 

9.0552 


3.32 
- .403 
9„6049 





9„1521 




.30 
.125 
9.0966 


- .40 

- .443 
9„6463 







1.484 
0.1715 


-2.581 
0„4118 


- .339 
9„5302 


- 1.435 
0„1569 


1 


0„2403 




6.097 

4.489 

0.6521 


.206 

.589 

9.7705 


3.373 

2.496 

0.3972 


1 


9„3587 




.175 
.077 

8.8889 


.043 

.328 

9.5156 


1 


9„3769 = log-Zi 


1.608 
0.2063 


- .383 
9„5832 


1.227 

0.0888 


2 




.098 

.091 

8.9601 


- .285 

- .292 
9„4657 





0.238 = Z4 \_vv'\ 


= .007 


.007 


.007 
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Scheme B. — Solution op the Eliminatiox Equations. 



- 1>^ 


-Caz, 


- ^Baz, 

- JBcZ, 


-Am 

- Aa 84 

- -^c 23 

- ^b 22 


S4 

log 84 


log 23 


2s 
log 22 


2i 


logC, 
log-Bd 
log^d 


log.B« 

log^o 


log ^6 




log ^6 22 


log -Sc 23 

log Ac Z3 


log C<jS4 

log^d 24 

log Aa Si 







114. Filling out this table for the example just solved, we 
have 



.238 


.228 
.414 


.142 
.031 
.145 


.260 
.143 
.514 
.191 


.238 
9.3769 


.642 

9.8075 


.318 
9.5024 


1.108 


0„2403 
9„1107 
9„7782 


9„3537 
9„9031 


9„7782 




9„2806 


9„1612 
9„7106 


9„6172 
8„4876 
9„1551 
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115. The Weights of the XJnknowns. In order to deter- 
mine the precision measures of gj, gji and z^, it would 
next be necessary to compute the weights of the latter quan- 
tities. The demonstration of the processes by which these 
weights may be found will not be taken up here, as the best 
method to adopt varies a good deal with the character of the 
example, but a statement of the results in the general form 
of solution will be given. 

By treating the Elimination Equations in a way similar to 
that used in deriving equation (E) of paragraph 76 from 
equations (B) of the same paragraph, we may show that 

^1 + -^m + -Sm cii + C^m "2 + -Z^m a, = 

z, + B„ + C„. A + D„y3a = (100) 

3. + t\ + B^y, = 

where the a's, /8's, y's, are determined from the equations 

Aa + Baa, + Caa, + a, = -B^ + Q yS, + /Ja = 
^, +^,ai+a, = ^„+^, = (101) 

^j -f- ai = Q + y, = 

Then by an application of the principles of Rule 1, para- 
graph (77), it may be shown that 



1 


= 


1 

\_aa'] 


+ 


ai' 




+ 


Pz, 


[66, 


1] 


1 


= 






1 




+ 


P', 


[66, 


1] 


I 

Ph 


z= 












1 
P'. 


= 













a.» , aa» 



ice, 2] + Idd, 3] 

(102) 



Ice, 2] + [<M73] 



1 , -H.' 
Ice, 2] "^ Idd, 3] 

1 
Idd, 3] 
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These equations can of course be extended to cover any 
number of unknown quantities, and tabular schemes for the 
computations of the a's, ^'s, y's, . . . and of the weights 
Pz^i Pz^^ Pz^t • ■ ■ can readily be arranged. 

For a general demonstration of these results, and also for 
a discussion of special methods of solution, consult — 

Johnson, "The Theory of Errors and Method of Least Squares," 

chap. ix. 
Wright, " Treatise on the Adjustment of Observations," chap. iv. 
Chawienet, " Spherical and Practical Astronomy," pp. 530-649. 



THE METHOD OF CORRELATIVES. 

116. The method of adjusting "Conditioned Observa- 
tions " explained in paragraph 33 is perfectly general, but 
where there are many conditions to be satisfied the solution 
is apt to be very laborious. For the case that occurs most fre- 
quently in practice, in which the observations are direct and 
equal in number to the number of unknown quantities, the 
process of solution devised by Gauss and called the " Method 
of Correlatives " is the most convenient. This method is 
derived as follows : — 

Let q observations, M^, M^, . . . Mg, of the respective 
weights pi, P21 • • • Pqi I'e made directly upon the values 
of q unknown quantities, and let the most probable values 
of the unknowns be 

2l = J/i 4" ''d 32 = ^2 + «3> . . . S5 = -S/g -f- Vq. 

"Where v^, v^, . . . Vg, are the most probable corrections to 
apply to the observed values as well as in this case the 
residuals of the observations. 

If the n' condition equations are not linear they may be 
reduced to that form by the method of paragraph 44, so that 
we may assume for our 
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CONDITION EQUATIONS 

«iWi + «2«a + • • • <^g'"q -\- mi = 
biVi -\- b^Vi + • • • *9^s -f Wj, = 



(A) 



ZjWl + l^Vi -\- . . . IqVq -f- »n„' = 

In which the quantities mj, wij, . . . ?w„', would all he zero 
if the observations were exact. It is to be observed that the 
coefficients a, b, ... I are not arranged in the same order in 
these equations as they are in the observation equations of 
paragraph 107. 

The values of Wj, v^, ... Vq, must be determined so as to 
satisfy the above equations and also by the principle of Least 
Squares, so as to make 

PiVi -\- p^v^ -|- . . . PqOq^ = a minimum. 
Corresponding with a minimum value of this last we have 

PlVidVi -\- p^V^dVi + . . . PqliqdVq = (B) 

for all possible simultaneous values of dv-i, dv^, . . . dvg ; that 
is, for all values which satisfy the equations, 



a^dvi -\- a^dVi -|- . . . aqdvq = 
bidvi -j- b^dv^ -f- • • • bqdvq = 



lidvi -\- Ijdvj -{-... Iqdvq = 



(C) 



obtained by differentiating equations (A). 

Therefore, denoting the first member of (B) by R and the 
first members of (C) by /S^, S^, . . . S^', it will be nec- 
essary that 

M - kiSi - k^S^ - ... kn.Sn, = (D) 

where ^i, k^, ... A„., are undetermined coefficients. 
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This last equation will be satisfied if the coefficient of each 
differential in it is made equal to zero, that is, if 






PqOq = ^ittg + ^s^g + ... *„'/, 



(103) 



All that remains therefore is to find values of v^, v^, . . . Vq 
and ^1, Aj, ... k^,, which will satisfy simultaneously equa- 
tions (A) and (103), and that this may be done is easily 
seen from the fact that we have the same number of equations 
as unknowns. 

Substituting the values of Vi, v^, ... v„ from equations 
(103) inequations (A) we have the following : 



, _ aa , , ^ ab , , „ al , 

A;x S — + fcj S— + . . . A!„. S — + wi = 

, „ ab , , _ &6 , ■, ^ bl , 

fcx S — - + A;^ S — +...*„' S — + m. = 



p ' ' ' p ' "«■ - p 



(104) 



fci S — + &, S — + . . . A!„. S - + m„' = 



The solution of these equations gives at once the values of 
Ai, ^2, ... kn', which are called the "Correlatives" of the 
Condition Equations. The values of Vi, v^, . . . Vg are then 
found by substituting the values of the k's in equation (103). 

117. As equations (104) are of the same general form as 
a set of Normal Equations, Gausses Method of Substitution 
can be advantageously employed in the solution. 
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When there is but a single equation of condition the second 
members of equations (103 J reduce to their first terms, and 
equations (104) reduce to the single equation 

aa 

k^% [- wii = (105) 

and the values of Vi, v^, ... Vq in (103) become 



Vi = — nil, Vi ^ ^-— niu 



CTi 


an 


aa 
^ P 


Pi 

Vo ^ 

aa 
^ P 




Clq 


Vg = 


aa 
^ P 



(106) 



It is from these results that the rules in paragraph 35 are 
derived. 

118. Example. Suppose we have given the observations 



(a) 



and let the most probable values of the unknowns be repre- 
sented by 

Zi = JIfi + «i, Zj = Jfa + «j, . . . 25 = Jfj + W5 (b) 

Also suppose that the unknowns are subject to the conditions 

Si + 2j + Z3 + Z4 -f Z5 = 14.0 

Sa - s^ = 1.5 (•"•' 



M^ = 2.02, 


weight 3 


M^ = 4.13, 


" 2 


M, = 2.52, 


" 5 


M, = 2.67, 


" 7 


M, = 2.84, 


" 4 
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Then expressing these conditions in terms of the corrections 
by means of (a) and (b), we have the 



(d) 



CONDITION EQUATIONS 

•Oi -\- Vi + v^ + v^ ^ v^ -\- .18 = 
■yj — Vi — .04 = 

Referring to paragraph 116, we see that in this example 

n' = 2, »ii = 0.18, m^ = - 0.04 



For the purpose of computing the coefficients in equations 
(104) we next arrange the following table. 



p 


a 


b 


a a 
~P 


ab 


bb 
P 


3 
2 
5 

7 
4 


1 
1 

1 
1 

1 




1 


- 1 




1 

1 
1 

2 
1 

1 

7 
1 

4 




1 
2" 



1 

7 






1 
2 


1 

Y 









599 

420 


5 
14 


9 

14 
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Substituting these results in equations (104), we have 



599 5 

— Ai -I ;fcj 4- .18 = 

420 ^14 ^ 



— Ai 4- — ^j - .04 = 
14 ' ^ 14 

Solving, Ai = - ,1647 

A, = .1537 

Then from equations (103) we get at once 



(f) 



Vi = — .0549, 1)4 = — .0455, 

Wj = — .0055, Vi = — .0412. (g) 

«s = —.0329, 

And by substituting these results in (b) we can obtain the 
values of Zj, z^, z,, z,, Zj. 

The above is the solution of Example 70, page 127. 



EXAMPLES. 



1. An urn contains five black balls, three red balls and 
two white balls. If three balls are drawn from the urn what 
different combinations may result, and what is the probability 
of each ? 

2. In a single throw with a pair of dice what is the 

5 
probability that neither ace nor doublets will appear ? 

y 

3. Four cards are drawn from a pack. What is the 
probability of getting four aces? Of getting one of each 
suit? 

4. From a lottery of thirty tickets, marked 1, 2, . . . 30, 
four tickets are drawn. What is the probability that num- 

2 
bers 1 and 15 will be among them ? - — 

^ 145 

5. Find the odds against the appearance of 7 or 11 in a 
single throw with a pair of dice. 7 : 2 

6. I toss up n coins. What is my chance of getting just 
one head ? 

7. In a single throw what are the relative chances of 
throwing 9 with two dice and with three dice ? 24 : 25 

8. From 2 n counters marked with consecutive numbers 
two are drawn. What are the odds against having an even 
sum ? n : n — 1 

9. In two trials with a single die what is the probability 
of throwing (a) an ace the first time only? (J) at least one 
ace? 

10. Find the probability of throwing doublets one or 

91 
more times in three trials with a pair of dice. 

216 
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11. Find the probability of throwing exactly three aces 

125 
in five trials with a single die. 

3888 

12. A certain stake is to be won by the first person who 
throws 5 with a die of twelve faces. What is the chance of 
the sixth person ? 

13. A and B play chess. A wins on the average two 

games out of three. What is A's chance of winning just 

80 

four games out of the first six ? 

^ 243 

14. A and B shoot alternately at a mark. A hits once in 
n times and B once in w — 1 times. Find their chances of 
first hit, and the odds in favor of B if A misses on his first 
shot. Even, n : n — 1 

15. In how many trials will it be a wager of 4 to 3 that 
double five will be thrown with a pair of dice ? 30 

16. Find the probability of throwing one and only one 

5 

ace in two trials with a single die. — 

18 

17. If I have three tickets in a lottery of four prizes and 

41 

eight blanks, what is my chance of drawing a prize ? — 

55 

18. Find the probability of throwing at least four aces in 

six trials with a single die. 

23328 

19. On an average seven ships out of eight return to port. 

Find the chance that out of five ships expected at least three 

.„ ^ 16121 
will return. 

16384 

20. In a lottery containing a large number of tickets, 
where the prizes are to the blanks as 1 : 9, find the chance 

of drawing at least two prizes in five trials. 

100000 
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21. In a purse are ten coins, all nickels except one which 
is a five-dollar gold piece ; in another are ten coins, all 
nickels. Nine coins are taken from the first purse and 
placed in the second, and then nine coins are taken from the 
latter and placed in the former. If you now had your choice 
which purse would you take ? 

22. A and B engage in a game in which A's skill is to 
B's as 2:3. What is A's chance of winning at least two 
games out of five ? 

23. If A's skill at a certain game is double that of B, 
what are the odds against A's winning four games before B 
wins two ? 131 : 112 

24. A party of twenty-five take seats at a round table. 
What are the odds against any two specified persons sitting 
next to each other ? 

25. A has three shares in a lottery in which there are three 
prizes and six blanks. B has one share in another where 
there is but one prize and two blanks. What are their 
relative chances of getting a prize? A : B ^ 16 : 7 

26. Expand through the terms involving A' and k^, the 
expression 

— ^ + (2/ + *)" 

When a is 1 and y is -, does 1- y' increase or 

2 X 

diminish when x and y begin to increase at the same rate ? 

27. Given f {x, y) = x^{a -f y)\ expand the expres- 
sion (X + A)=(a + y -f ky. 

28. Find the value of logio a -\- cos b, when a == 1001 
and b = 0.1°. Give the result first to five places and then 
to seven places of significant figures, in each case without the 
aid of tables. 4.000433 

29. Transform to the new origin, (2, — 3, 1), the equa- 
tion, x» + y' + z* - 4 a; -f 6 2/ - 2 z - 11 = 0, the 
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axes remaining parallel to the original ones. The eqnations 
of transformation are x = 1 -\- x' , y = — 3 -\- y' , 
2 = 1 + z'. SB" + y" + 2= = 25 

30. Find the minunmn value of x^ -\- xy -\- y^ — ax — by. 

(ab — a' — b^) 
3 

31. Find the values of x, y and z that render a maximum 
or a minimum the function x^ -\- y^ -\- z" -\- x — 2 z — xy. 

32. Find the values of x and y that render a maximum 
or a minimum the expression sin x -\- sin y -\- cos (x -\- y). 

33. Find the co-ordinates of a point the sum of the squares 

of whose distances from three given points, (Xi, y,), (x^, 2/2)) 

s • • • (a^ + aij -I- Xs) 
(ajj, 2/3), IS a minimum. x = -^— = — ■ = — ■ — -: 

o 

34. Given the volume, a\ of a rectangular parallelopiped, 
find its shape when its surface is a minimum. 

35. Find the volume of the greatest rectangular parallelo- 
piped that can be inscribed in the ellipsoid 

x^ y^ z'^ Sabc 

«' b' c^ 3^ 

36. In the Physical Laboratory apparatus for illustrating 
the Estimation of Tenths, the reading of the micrometer head 
for a certain setting is computed to be 2.3038, which may be 
taken as exact. Setting the apparatus by the eye the follow- 
ing readings are obtained : — 

2.314 2.324 2.310 2.519 2.326 
2.320 2.302 2.313 2.305 

What are the accidental and real errors and the residuals ? 
Are there any constant errors or mistakes ? If there are 
constant errors are they of the first, second, or third class? 
How do you tell ? 
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37. Are the following observations such as to call for an 
application of the Method of Least Squares in their adjust- 
ment? 

x-\- y — z-\-u ■= 5 3x — y— 2— M=l 

X — 2y -{- 1z -\- u z= —1 SB — y + 4z-f6M = 9 

38. In the case of direct observations, what other quanti- 
ties besides the arithmetical mean might reasonably be 
assumed to give plausible values of the unknowns ? Why is 
the arithmetical mean preferred to these ? 

39. Find the most probable value of a quantity M from 
the observations 

216.58 



216.27 


216.16 


216.04 


216.19 


216.44 


.29 


.43 


215.99 


.39 


.51 


.33 


.09 


216.23 


.14 


215.94 



Also test the result by finding the sum of the residuals. 

M = 216.251 

40. In the determination of a certain wave length, Row- 
land made the following observations. Find the most prob- 
able wave length. 

4.524 4.515 4.513 4.507 4.501 4.485 4.517 4.493 4.505 
.500 .508 .511 .497 .502 .519 .504 .492 

4.5055 ± 0.0017 

41. Ten measurements of the density of a body gave 
results as follows. Find the most probable density. 



9.662 


9.664 


9.677 


9.663 9.645 


.673 


.659 


.662 


.680 .654 

9.6639 ± 0.0022 



42. In a triangulation of the U. S. Coast Survey an angle 
was measured twenty-four times with results 
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116° 43' 



44" 


.45 


48".90 


47".40 


47".85 


51".75 


50 


.55 


49 .20 


47 .75 


50 .60 


49 .00 


50 


.95 


48 .85 


51 .05 


48 .45 


52 .35 






51".30 


49".05 


46".75 








51 .05 


50 .55 


49 .25 








51 .70 


49 .25 


53 .40 





Find the best value for the angle. 116° 43' 49".64 ± 0".28 

43. The following observations were made with a sextant 
in order to determine the latitude of a place : — 

43° 4' 

What is the most probable latitude? 43° 4' 27" ± 5" 

44. Determine the quantity M from the observations 



46" 


7" 


59" 


52" 


47" 


36" 


24 


28 


39 


52 


15 


40 



M = 


81.55 


.41 


.68 


.25 


.27 


.77 


.13 


.86 


.10 


.03 


Wts. 


25 


25 


16 


9 


9 


9 


4 


1 


1 


1 



45. An angle M is measured with the following results : 

M p M p M p M p M p 

45".00 5 42".50 5 27".50 3 36".25 2 45".00 2 

31 .25 4 37 .50 3 43 .33 3 42 .50 3 40 .83 3 

45 .00 3 38 .33 3 40 .63 4 39 .17 3 

Find the most probable value of the angle and test the 
result by computing 2/>«. 39".78 ± 0".94 

46. In one hundred measurements of angles made in the 
primarry triangulation in Massachusetts by the Coast Survey 
there were found between 



6".0 


and 


5'' 


'.0 


1 


error. 


1' 


.0 


and 0".0 


26 errors. 


5 .0 


(I 


4 


.0 


2 


errors. 





.0 


" - 1 .0 


26 " 


4 .0 


(I 


3 


.0 


2 


(t 


-1 


.0 


" - 2 .0 


17 " 


3 .0 


(C 


2 


.0 


3 


(( 


-2 


.0 


" - 3 .0 


8 « 


2 .0 


<1 


1 


.0 


13 


t( 


-3 


.0 


« — 4 .0 


2 " 



Plot these results as in Figure 1, page 8. 
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47. In sixty-six determinations of the velocity of light 
made at Washington the percentage of errors of different 
magnitudes was found to be as follows : — 



Over - 7.0 


.8% 


Over 


7.0 


0.0% 


Equal to - 7.0 


.8 


Equal to 


7.0 


1.6 


-6.0 


1.6 




6.0 


3.5 


-5.0 


2.0 




5.0 


4.2 


-4.0 


3.9 




4.0 


5.8 


-3.0 


7.8 




3.0 


8.6 


-2.0 


12.8 




2.0 


12.4 


— 1.0 


17.1 




1.0 


17.1 



Plot the results and draw a smooth curve. 

48. Given the observations 

X — y + 2s = 3 4a: + y -f 4z = 21 

3a; + 22/ - 5a = 5 — a; -}- Sy -|- Sz = 14 

find the most probable values of x, y, and z. 

X = 2.470 ± .038 

49. being the level of the sea and Pi, Pj, Pg, three 
points whose altitudes are to be determined, the following 
observations are made : — 



Pi above = 10 ft. 


Pj above P, = 9 ft. 


P, " Pi = 7 


P, « Pg = 2 


Pj " = 18 





Find the most probable altitudes. P, = 8.50 ± .29 

50. The altitudes of A above 0, B above A, and B 
above are found by measurements to be respectively 
12.3, 14.1, and 27.0 feet. What is the most probable value 
of each of these differences in level? A = 12.50 ± .17 

51. Measurements of the ordinates of points on a straight 
line corresponding to abscissas 4, 6, 8, 9, are made with 
results 5, 8, 10, 12. What is the most probable equation of 
the line in the form y = mx -\- bf b = — 0.29 
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52. Find the altitudes in Example 49 if the observations 
have weights 5, 3, 6, 2, 4, respectively. 

53. Solve the example in paragraph 31, giving the obser- 
vations the weights 25, 25, 4, 4, 4, 4, 4, 4, 1, respectively. 

Elevation of A = 320.25 

54. Find the most probable values of s,, z^, and a;, from 
the observations 

2i = 552.10 wt. 16 Si 

— 22 -f Zj ^ .15 " 9 Zi + Zj 

Zj = 551.28 « 4 Zi 
z„ = 551.30 " 4 



55. In the triangalation of Lake Superior there were 
measured at station O the angles 

F P = 62° 59' 40".33 wt. 5 



S2 


= 


.75 


wt. 


1 


S3 


= 


552.05 


C( 


1 


Sa 


= 


.70 


C( 


1 




Zj 


= 551.2345 



F B = 64 


11 


34 


.92 


u 


7 


FOB = 100 


20 


29 


.12 


u 


4 


P B = 37 


20 


49 


.55 


(C 


7 


B, B = 36 


8 


55 


.86 


(( 


4 



Required the adjusted values of the angles. 

F P = 40".28 ± 0".34 
56. In the U. S. Lake Survey the following angles were 
measured at station North Base : — 

(1) Crebassa — Middle 55° 57' 58".68 wt. 3 

(2) Middle — Quaquaming 48 49 13 .64 " 19 

(3) Crebassa — Quaquaming 104 47 12 .66 " 17 

(4) Quaquaming — South Base 54 38 15.53 " 13 

(5) Middle — South Base 103 27 28.99 " 6 

Find the adjusted values of the angles. 

(1) = 58".965; r = 0".28 
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57. Adjust the following obBervations of difEerences in 
level : — 

Altitude of A 401.3 wt. 16 C above B 12.b wt. 9 

A above B 220.8 " 16 A " B 222.0 " 1 

^ " C 150.2 " 4 Altitude of B 180.7 " 1 

58. In "Conditioned Observations" can the number of 
observations required be less than the number of unknown 
quantities? Why must the number of conditions be less 
than the number of unknowns? 

59. From the following measurements of the angles formed 
at the centre of a disk by four radial lines, find the most 
probable values of the angles. 



A = 104° 25' 13" 


C = 86° 33' 20" 


B = 98 13 47 


2) = 70 48 23 




A = 104° 25' 2".25 



Also solve giving the observations the weights 5, 2, 1, 4, 
respectively. 

60. Four observations on the angle -4 of a triangle gave 
a mean of 36° 25' 47", two observations on B gave a 
mean of 90° 36' 28", and three on C gave 52° 67' 57". 
Adjust the triangle. A = 36° 25' 44".2; r = 7".7 

61. Five angles at a station are measured, and also their 
sum. The observed sum differs from the sum of the five 
observed parts by the amount d. What are the adjusted 
values of the angles ? 

62. The three angles of a spherical triangle are measured 
with results 

A = 46° 17' 38".32 B = 73° 35' 16".15 
C = 60° 7' 5".16. 

Adjust the triangle, knowing that the spherical excess is 
2".475. A = 39".3; ^ = 1".6 
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63. At the station Pine Mountain the following angles 
were observed between surrounding stations : — 

Jocelyne — Deepwater 65° 11' 52".500 wt. 3 

Deep water — Deakyne 66 24 15 .553 " 3 

Deakyne — Burden 87 2 24 .703 " 3 

Burden — Jocelyne 141 21 21 .757 " 1 

Find the most probable values of the angles. 

64. Solve Examples 55 and 56 by the method of "Con- 
ditioned Observations." 

65. A la a. station whose altitude is known to be 5240.1 
feet. -B and are floats on a lake, and Z> is a signal point. 
From the following observations determine the most prob- 
able altitudes of B, C and D. 



C below A 


720.1 


wt. 3 


B below A 719.7 


wt. 3 


D ^' A 


200.3 


" 5 


B '' J) 520.9 


« 2 


C <■' I) 


520.4 


" 2 







66. Given the following observations, subject to the con- 
dition Si -|- Sj = Sj, find the most probable values of 
Zxi ^2, and S|. 

2zi — s, + 2, = 3.0 2z2 — 2, = 1.0 

2zi — 322 = — 4.5 2, -f 2zj = 5.1 

Zi 4- z, = 3.8 

67. The chemical composition of a specimen was found 
by several observers to be as follows : — 

Pb = .52 Other substances = .09 Au and Ag = .39 
Ag = .27 Fb and Ag = .78 Impurities = .10 

Au = .11 P6 and impurities = .62 Au =: .12 

From these observations find the most probable composition 
of the specimen. 



- Washington 


23» 


41^041 


wt. 30 


- Cleveland 


42 


14.875 


" 7 


- Columbus 


47 


27.713 


" 8 


— Columbus 


23 


46.816 


" 7 


Columbus 


5 


12 .929 


« 5 
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68. From the following observations what are the best 
values of the unknowns, supposing that y and z must be 
equal ? 

a: + y = 5.2 wt. 4 y + s = 4.2 wt. I 
a: = 3.0 " 9 3 = 2.0 « 4 

X — z = 1.1 "1 

69. In determining the difference in longitude between 
various cities the results obtained were 

(1) Cambridge — 

(2) Cambridge — 

(3) Cambridge — 

(4) Washington ■ 

(5) Cleveland — 

Adjust these observations. 

70. The capacity of a condenser is known to be 14.0 m. f. 
It is divided into five sections, a, b, c, d, e, and it is known 
that the difference between b and d is 1.5 m. f. Find the 
most probable capacities of the sections from the observa- 
tions 

a = 2.02 wt. 3 d = 2.67 wt. 7 

b = 4.13 "2 e z= 2.84 " 4 

'^ = 2.52 "5 „ _ 19651 

71. If the unknowns in the following observations are 
subject to the condition x -\- 2y -\- Sz = 36, what are 
their adjusted values ? 

X = 4.3 wt. 1, y = 5.7 wt. 4, z = 7.3 wt. 9 

X = 3.77 

72. A cannon is discharged horizontally from the top of 
a bluff. Observations on the time, and distance of fall of the 
ball gave the results 

t = 0.5 1.0 1.5 2.0 seconds 
s = 1.2 4.0 9.1 15.0 metres 
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What curve, passing through the point of departure of the 
ball, will represent the above observations ? 

73. An Argand burner shows the following efficiencies 
with varying rates of gas consumption : 

g = 2.0 2.3 2.8 8.3 4.0 4.5 5.0 feet 
B — 2.1 2.4 2.5 3.0 3.2 3.8 4.1 

Find the equation of the straight line which best rep- 
resents the relation between ff and £J. The measurements 
on ff are without appreciable error. 

74. Observations are made upon the expansion of Amyl 
alcohol with change in temperature as follows : — 



V = 1.04 


1.12 


1.19 


1.24 


1.27 


cu. cm. 


t = 13.9 


43.0 


67.8 


89.0 


99.2 


C. degrees 



If V ^ 1 -\- J3 t -{- C t^ expresses the law connect- 
ing the volume and temperature, find the most probable 
values of £ and C. 

75. In a Hooke's joint where the angle between the axes 
is 45°, X being the angular rotation of the driver, and y 
that of the follower, from the following measurements find 
the equation of a curve that will represent the relation 
between x and y — x. 



X 


y — X 


X 


y—x 


x 


y—x 


0° 


0°.0 


80° 


- 5°.8 


140° 


8°.8 


20O 


- 5°.7 


90° 


- 2°.0 


160° 


5°.3 


40° 


— 9°.9 


100° 


2°.3 


180° 


0°.0 


60° 


- 10°.4 


120° 


8°.0 







y - X = — 0.85 — 9.82 sin 2a! + 0.92 cos 2x 

76. A series of observations extending over a period of 
thirty years was made by Quetelet to determine the daily 
variation in temperature at Brussels. The mean results of 
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the measurements are given below. From them derive an 
equation to express the temperature at any time of the year. 



Jan. 


4°.66 


May 9°. 83 


Sept. 8°.16 


Feb. 


5°.42 


June 10°.09 


Oct. 6°.55 


Mar. 


6°.77 


July 9°.71 


Nov. 5°.10 


Apr. 


8°.59 


Aug. 9°.14 


Dec. 4°.41 




y = 


7.369 + 0.9854 sin 


30a! - 2.7084 cos30£b 






4- 0.0100 sin 


60a;- 0.1950 cos 60a; 






- 0.0133 sin 


90a; 4- 0.1788 cos 90a; 



In this answer the values of x begin at the 15th of Janu- 
ary, and represent the time in months. 

77. The law connecting the time of vibration of a pendu- 
lum with its length is assumed to be of the form, T ■=^ m i". 
From the following observations find the most probable values 
of m and n. 



T = 12.9 


11.6 


10.4 


9.7 


5.3 


4.6 


L — 164.4 


132.9 


107.6 


93.5 


28.4 


20.6 



L is in centimetres, T in tenths seconds. n = 0.5000 

m — 1.0044 
78. Determine the equation of a curve which will repre- 
sent the following observations : 



X 


y 


X 


y 


X 


y 


0.0 


0.00 


1.5 


1.09 


3.0 


8.65 


0.5 


0.04 


2.0 


2.56 


3.5 


13.72 


1.0 


0.31 


2.5 


4.99 


4.0 


20.47 



79. Determine the equation of a curve which will repre- 
sent the relation between x and y in the observations. 



X 


y 


X 


y 


X 


y 


X 


y 


0.0 


4.51 


0.3 


4.09 


0.6 


3.03 


1.2 


0.92 


0.1 


4.44 


0.4 


3.76 


0.8 


2.24 


1.5 


0.38 


0.2 


4.31 


0.5 


3.42 


1.0 


1.49 


2.0 


0.05 



130 METHOD OF LEAST SQUARES. 

80. At a station P the angles between a straight line 
passing through P parallel to the axis of ^ and the direc- 
tions from P of four points Pi, Pj, Pj, P4, are measured. 
Having given the coordinates, (a, b), of the four points, 
find the coordinates of P. 



^oint. 


Coordinates. 


Angle. 




a 


6 




A 


4.21 


3.24 


39° 18' 


P, 


1.21 


2.10 


147° 54' 


^. 


— 0.51 


0.22 


205° 24' 


A 


2.50 


- 1.10 


277° 15' 



If the coordinates of the point P are (a:, y), and the 
angle is denoted by A, we have 

y — b 

t&uA = 

X — a 

81. If a sin 5a; = M, and values of M are observed 
for known values of a and b, determine the most probable 
value of X. If x' is an approximate value of x found by- 
trial, and m = a sin bx' — M, we shall have 

, 'S.a b m cos bx' 
^ — ^ — -^i^ah cos bxf 

82. If in one series of observations the value of h is 
twice what it is in another, what is the relative probability of 
the occurrence of an error of given magnitude a in the two 
series ? Show what the curves of error will be in the two 
cases. What error has the same probability for its occurrence 
in each series ? What is the relative probability of the occur- 
rence of an error not greater than a in the first case and not 
greater than 2a in the second case ? 

83. From 64 observations the latitude of a station was 
found to be 49° 10' 9".110 ± 0".051. What was the prob- 
able error of a single observation ? 0".41 
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84. If twenty measurements of an angle give a result with 
an A.D. of 0".38, and it is required to find the angle so 
that the A.I), shall be only 0".25, how many more observa- 
tions must be made ? 27 

85. From the following determinations of the area of a 
field find the most probable area and its probable error. 

5674 ± 12, 5680 ± 4, 5685 ± 3, 5682 ± 1, 5678 ± 2 

A = 5681.41 ± 0.84 

86. From the following measurements by Fizeau and 
others, find the most probable value for the velocity of light 
together with its probable error. Measurements are in kilo- 
meters. 

298000 ± 1000 299990 ± 200 299930 ± 100 
298500 ± 1000 300100 ± 1000 

V = 299917 ± 88 

87. Two different instruments give for the value of an 
angle, 

34° 55' 33".0 ± 4".l, 34° 55' 36".0 ± 6".3 

What is the best value to take for the angle ? 

34° 55' 33".9 ± 3".4 

88. Determinations of the difference in longitude between 
Washington and Key West made on seven different days 
gave the results 

19"* 1'.42 ± 0*.044 19*" 1'.60 ± 0'.046 

1 .37 ± .037 1 .55 ± .045 

1 .38 ± .036 1 .57 ± .047 
1 .45 ± .036 

What is the best value and its probable error ? 

i'.46o ± o'.oie 
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89. In the triangulation of Lake Ontario two different 
instruments gave for an angle, 74° 25' 5".429 ± 0".29 from 
sixteen readings, and 74° 25' 4".611 ± 0".22 from twenty- 
four readings. Find the most probable value of the angle 
and its probable error. 

90. In each of Examples 39-45 find the mean and prob- 
able errors and average deviation of each observation and of 
the most probable value, using formulas from (50) to (63) 
according as they apply. 

91. In Example 42 divide the observations in their order 
into six groups of four observations each and compute the 
mean of each group. Then determine the probable error of 
tlie first of these means: (1) considered as a single measure 
of four times the weight of those in Example 42 ; (2) directly 
as one of six observations of equal weight; (3) as a deter- 
mination from its four constituents. 0".67 ; 0".72 ; 1".00 

92. The following twenty-nine measurements on the den- 
sity of the earth, made by Cavendish, give as a, mean result 
5.48. What is the probable error of an observation ? Solve 
by the usual method and also by taking the residual that 
occupies the middle position. 0.14 



5.50 


5.55 


5.57 


5.34 


5.42 


5.30 


.61 


.36 


.53 


.79 


.47 


.75 


.88 


.29 


.62 


.10 


.63 


.68 


.07 


.58 


.29 


.27 


.34 


.85 


.26 


.65 


.44 


.39 


.46 





93. What is the probable error of the mean of two obser- 
vations which differ by the amount a ? 

94. A base-line is measured five times with a steel tape 
reading to hundredths of a foot, and five times with a chain 
reading to tenths of a foot, with results 



By tape, 


741.17 


741.09 


741.22 


741.12 


741.10 


By chain, 


741.2 


741.4 


741.0 


741.3 


741.1 
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Find the probable errors and weights for a single observa- 
tion in each case, and also the adjusted length of the line 
and its probable error. 741.146 ± 0.015 

95. Twenty-one determinations of a chronometer correc- 
tion gave results 

- 8.78 - 8.78 - 8.68 - 8.80 - 8.96 - 8.83 - 8.79 
.76 .51 .63 .75 .64 .70 .90 

.85 .64 .58 .78 .65 .64 .93 

Find the probable error of the mean by using both formulas 
(53) and (57), and also determine the probable error of a 
single observation by taking the middle residual. 

0.017; 0.018; 0.09 

96. In the following observations show that Mg := 49.64, 
^ — 1.95, r = 1.31, ^ = 0.40, ^ = 0.27, fi^ = 0.87, 
rg = 0.59. 

M = 48.81 48.76 49.53 51.56 50.38 49.84 
p = 5 4 5 3 2 5 

97. Observations on the time of ending of a transit of 
Mercury are made by different observers with a variety of 
instruments and under more or less favorable circumstances. 
If the weights assigned by the computer are as indicated, find 
the best value for the time and its probable error. 



5* 38"* 23' 


wt. 1 


38™ 26* 


wt. 3 


38" 19' wt. 3 


37 55 


« 


38 21 


" 2 


38 21 « 2 


38 10 


« 1 


38 18 


« 2 


38 15 « 2 
to = 5* 38" 19'.9 



98. An angle is measured five times with a theodolite, and 
seven times with a transit, giving results 

Theodolite, 31".7, 39".8, 40".7, 28".6, 32".3 

Transit, 32 .8, 36 .7, 38 .2, 29 .3, 41 .6 35".3, 36".2 
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If the relative values of readings by the two instruments 
are as 3 to 2, what is the most probable value of the angle ? 
What is the mean error of the result ? 

99. Given 3I-, = 65.58 ± .59, M^ = 35.15 ± .93, 
Ms = 49.64 ± .27, find the probable errors of 4 Jfj — 

3 J!f2 + 2 Jfs and of ^ + ^ ■- ^. 3.69 ; 0.43 

ii O *i 

100. The three angles of a triangle are measured, and the 
probable error of each observation is r . What is the prob- 
able error of the triangle error ? r f/W 

101. The zenith distance of a star on the meridian is 
observed to be z = 21° 17' 20".3 ± 2".3. The declina- 
tion of the star is given as <? = 19° 30' 14".8 ± 0".8. 
What is the latitude of the place and its probable error ? 

i = g + (f = 40° 47' 35".l ± 2"4. 

102. The zenith distance z of a star at upper culmina- 
tion is observed n times, and its zenith distance z' at 
lower culmination ?i' times. If the latitude is given by 
X, = 90° — i (3 + s')j ^^^ *li6 probable error of an 
observation is r, what is the probable error of the latitude ? 

103. The horizontal force necessary to start a 100-pound 
•weight sliding along a table is observed to be 15.5 ± 0.2 
pounds. Find the probable error of the coefficient of friction. 

104. If a line is measured by the continued application of 
a unit of measure, and r is the probable error of the placing 
and reading of this measure, what is the probable error of 
the length Z ? ^ yT^" 

105. If the average deviations of Zj, Zj, z,, are a, b, c, 
respectively, what is the average deviation of Sj^ -|- z^' -f- ^s' ? 

106. If the radius of a circle is measured with result 
1000.0 ± 2.0, how should the circumference and area be 
expressed ? 

107. Two sides, a and b, and the included angle C of 
a triangle are measured with results a ^ 252.52 ± .06 
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feet, b = 300.01 ± .06 feet, C = 42° 13' 00" ± 30". 
What is the area and its probable error ? 25452 ± 9 

108. Measurements of adjacent sides of a rectangle gave 
a ± ri, and b ± r^. What is the probable error of the 
area, and for what kind of a rectangle will this probable 
error be the least ? 

109. If the measured sides of a rectangle have the same 
a.d., what is the a.d. of the diagonal determined from 
them ? Same 

110. If the sides of a rectangle are measured in the 
manner indicated in Example 104 and found to be a and b, 
what is the probable error of the area ? 

111. The correction to be applied to a chronometer is 
found to be + 12"* 13'.2 ± 0'.3. Ten days later the cor- 
rection is again determined and found to be 12™ 21'.4 ± 0'.3. 
What is the mean daily rate and its probable error ? 

0'.820 ± 0'.042 

112. Measurements of the compression of the earth's 
meridian have resulted in 

± .000046 

294 

What is the probable error of the denominator 294 ? 3.98 

113. The current flowing in a circuit is due to two 
sources whose electromotive forces are determined to be 
Ci = 200 ±2, ej = 400 ± 3. The resistance of the 
circuit is 30 ± 1. Find the current and its probable error. 

20 ± 0.68 

114. The side b and angles -B and C of a triangle are 
measured with results b = 106 ± .06 metres, J3 = 
29° 39' ±1', C = 120° 7' ± 2'. What is the most 
probable value of the angle A and of the side c ? 

A = 30° 14' ± 2'.2; c = 185.35 ± 0.15 

115. The distance between two divisions on a graduated 
scale is measured by a micrometer. Show that the average 
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deviation of the mean of two results is the same as the aver- 
age deviation of a single reading. 

116. If the weights of the determinations of three angles 
A, B, C, are 3, 3, 1, respectively, what is the weight of 
the sum of the three angles ? 0.6 

117. If the weight of x is />, what is the weight of 
logo « ? 

118. If a; = — and the weight of y is p, what is the 
weight of X ? c^p 

119. In Example 107, how closely must the parts be 
measured in order to obtain the area within 0.5 per cent? 

120. From observations on I and t the value of g is 
to be computed by the pendulum formula 



v/4 



9 

What changes in g will be produced by changes in I and 
t of 8i and 82 units, respectively, and what are the allow- 
able errors in I and { if ^ is to be determined within 
1 per cent ? 

121. The moment of inertia of a cylindrical bar is to be 
obtained from measurements on its mass m, its length h, 
and its diameter d. The error in the determination of m 
is negligible, the precision of the determination of d is four 
times that of h. If the measurements give m = 48, 
h = 8.000, d = 1.200 ± 0.10, and 



I = ml U — 

' 12 ^ 16 



what is the probable error of I, and what should be the 
ratio of d to h to determine J most accurately ? 

d : h = 256 : 9 
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122. If observations give for a certain quantitj' x the 
value 303, with a mean error of 2, wh:it is the mean error 
of the expression 'dx -\- log^ 2 cc ? 

123. The probable error of the determination of the 
angle A is 20". "What is the maximum probable error of 
sin A -\- cos A ? 

124. If the probable error of an observation on an angle 
is 10", is there any difference between the probable error 
of the function sin A -\- cos A -\- sin C and of tlie func- 
tion sin A -f- cos B -\- sin C, supposing A and B are 
of the same magnitude ? 

125. Given the observations, 

2j - 222 + Ss - 3 = 331 + Sj + 2z3 - 17 = 
23i -f 32j - 4S3 - 2 = - si + 4s, + 323-10 = 

Find the most probable values of z^, z,, 23, and also their 
weights and precision measures. 

z, = 3.541; p,, = 29; »•„ = .024 

126. Find the weights and precision measures of the 
unknowns in Examples 48 to 57. 

127. Determine the probable errors of the constants in 
Examples 72 to 79, inclusive. 

128. The length of a pendulum which beats seconds is 
given by 

I =z I' J^ ( ^ g — s] r sin^i 



where I' is the length at the equator, a the ratio of 

^ 1 ' i 289 

the centrifugal force at the equator to the weight, and s the 
compression of the meridian regarded as unknown. Putting 

I' = 991'"-'"- + x, f^q - s\ r = y, 

observations in different latitudes gave in millimetres the 
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following equations, from which we are to determine I and 
« together with their probable errors : — 

X -\- 0.969y = 5.13 x + 0.152y = 0.77 

X -I- 0.749y = 3.97 x + 0.327y = 1.70 

X 4- 0.426y = 2.24 x -j- O.GSSy = 3.62 

X -\- 0.095y = 0.56 x + 0.793y = 4.23 

X = 0.19 

V = 991.069 ± .026; s = — ± 0.00046 

294 

129. Find the weights and precision measures of the 
unknowns in Examples 64 to 70, and also in Examples 59 to 
63, and in 71. 

130. In Example 95 the probable error of a single obser- 
vation is 0.08 seconds. Find the number of errors which 
should fall between 0.00 seconds and 0.10 seconds, between 
0.10 seconds and 0.20 seconds, and also the number that 
should be over 0.20 seconds. Compare the results with the 
number actually found. 

131. In 470 determinations of the right ascensions of 
Sirius and Altair made by Bradley, the probable error of a 
single observation was 0".2637. The number of errors falling 
between specified limits was as shown below. Compare this 
result with the distribution of errors called for by the theory. 

Limits. Errors. Limits. Errors. 



0".0 to 0".l 


94 


0".6 to 0".7 


26 


.1 to .2 


88 


.7 to .8 


14 


.2 to .3 


78 


.8 to .9 


10 


.3 to .4 


58 


.9 to 1 .0 


7 


.4 to .5 


51 


Over 1 .0 


8 


.5 to .6 


36 







132. What is the probability that the error of a single 
observation will be as large as twice the probable error ? 
As large as five times the probable error ? 
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133. On the average how many observations must be made 
before an error as large as three times the mean error will 
occur ? 

134. In Example 46, assuming that all errors between any 
two limits fall half way between those limits, compute the 
average deviation and mean error of an observation and com- 
pare their ratio with the theoretical value given in the table 
in paragraph 55. 

135. A line is measured 500 times and the probable error 
of each observation is 0.6 cm. How many errors should 
occur between 0.4 cm. and 0.8 cm. ? 

136. Show how the value of n- could be determined 
experimentally from observations such as those in Example 
181. 

137. In a system of observations all equally good, r being 
the probable error of a single observation, if two observations 
are taken at random, what quantity is their difference as 
likely as not to exceed, and what is the probability that the 
difference will be less than r? '>'^~2'; 0.367 

138. In the following measurements of an angle, ought 
any -of the observations to be rejected ? 

12' 51".75 47".85 47".40 48".90 44".45 
48 .45 51 .05 48 .85 50 .95 
50 .60 47 .75 49 .20 50 .55 

139. Determine whether any of the observations in 
Example 44 should be rejected. 

140. A quantity M is measured with the results given 
below. Ought all the observations to be retained ? 

M = 236, 251, 249, 252, 248, 254, 246, 257, 243, 274 

141. A certain angle has been laid out with such accu- 
racy that its true value may be taken as exactly 90°. Twenty- 
five observations are made upon it with a transit that it is 
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desired to test, and the result obtained is 89° 59' 57" ± 0".8. 
What are the odds in favor of a constant error in the instru- 
ment between - 1" and — 5"? Between 0" and — 6"? 

908 : 92 ; 86 : 1 

142. Repeated measurements of a standard metre bar 
with a decimetre scale gave a result 10.032 ± 0.010. What 
are the odds in favor of a constant error in the scale between 

^-I^and^:^? 43:7 

y/TO \JlQ 

143. Two determinations of the length of a line gave 
683.4 ± 0.3 and 684.9 ± 0.3, respectively. Show that the 
best value for the length is 684.15 ± 0.51, and that the 
probable systematic error of each determination is 0.65. 

144. Two men A and -B observe an angle repeatedly 
with the same instrument with results 

A. B. 

47° 23' 40" 28' 35" 47° 23' 30" 24' 00" 

23 45 23 40 23 40 23 20 

23 30 23 50 

Is there any relative personal error, and what is the best 
final value? 47° 23' 38".2 ± 1".6 

145. Three independent determinations of the capacity 
of a condenser made with three different instruments gave 
results 42.22 ± 0.21, 43.40 ± .15 and 44.20 ± 0.18. 
What is the most probable value of the capacity ? 

For extended treatment of the subject illustrated in Exam- 
ples 143 to 145 see Johnson, "The Theory of Errors and 
Method of Least Squares," chap. vii. 

146. In an estimation of tenths what is the probable error 
of an observation ? What is the average deviation ? 0.025 

147. In obtaining the angle of deflection of the needle of 
a tangent galvanometer by the usual method what is the 
probable error of the result ? 
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148. If all the errors of a series of observations must fall 
between and a, and the frequency of any error is pro- 
portional to its magnitude, what is the Curve of Error? 

What are the values of r, a.d., and u, ? r ^ 

sfT 

149. In Example 148 what is the probability that the error 
of a single observation will be as large as 0.5 a. 

150. If all values of x between and a are possible, 
and their probabilities are proportional to their squares, find 
the mean value of x and the probability that x will be as 
large as 0.5 a. Also draw the Curve of Error. 

151. What is the greatest probable error of a logarithm 
found by interpolation in a seven place table ? .000000015 

152. Given the following set of Normal Equations, together 
with [»iwij = 1.3409, find the most probable values of 
the unknowns and their weights and probable errors. There 
were sixteen observations. 

3.1217 Zi + .5756 s, - .1565 Zj - .0067 24 

- 1.5710 = 

.5756 Zi + 2.9375 s^ + .1103 z, — .0015 S4 

4- .9275 = 

- .1565 Si + .1103 Z2 + 4.1273 z^ -f -2051 z^ 

-j- .0652 = 
_ .0067 Zi - .0015 Zj + .2051 z^ -j- 4.1328 z^ 

-I- .0178 = 
Zi = 0.583 ± 0.018 
Z4 = - 0.004 ± 0.015 
p., = 4.12 

153. From ten observation equations, for which was found 
[mm] = 2.6322, there resulted the normal equations 

5.2485 Zi - 1.7472 z^ - 2.1954 z, + 0.5399 = 

- 1.7472 Zi -t- 1.8859 Zj + 0.8041 z, - 1.4493 = 

- 2.1954 Zi + 0.8041 Zj, -\- 4.0440 23 - 1.8681 = 
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Find the most probable values of Zj, Zj and Zj together with 
their probable errors. Zj = 0.42 ± 0.11 

154. Find the most probable values of the unknowns in 
the normal equations 

459 Zi - 308 Zj - 389 z, -f 244 z^ _ 507 = 

- 308 Zi + 464 Zj + 408 z, — 269 z^ + 695 = 

— 389 2i -j- 408 Zj -i- 676 z, - 331 z^ -f- 653 = 
244 Zx - 269 Zj - 331 z, + 469 z^ — 283 = 

[wm] = 1129 
Z4 = - 0.488 ; JO,, = 281 

155. If thirteen observation equations give rise to the 
result [mwi] = 100.34 and to the normal equations 

17.50 Zi - 6.50 z, - 6.50 z, - 2.14 = 

- 6.50 Zi -)- 17.50 Zj - 6.50 z, — 13.96 = 

- 6.50 zi - 6.50 Zj + 20.50 z, + 5.40 = 

show that the most probable values of the unknowns are 

Zi = 0.67 ± 0.60, Zj = 1.17 ± 0.60, z, = 0.32 ± 0.55 
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ELEMENTS OP THE THEORT OF PROBABILITY. 

200. Definition. If an event can happen in a ways, and 
fail in b ways, and all these ways are equally likely to occur, 

the probability of the happening of the event is i~x> 

and the probability of its failure is . 

Since the event must either happen or fail, the sum of the 
above probabilities must represent a certainty. But 

a b 

a -\- b + a -\- b ~ ^ 

That is, the probability of a certainty is expressed by unity. 
Also, if the probability, P, of the happening of an event is 
known, the probability of its failure is given at once by 
1 - P- 

201. Example A. A single throw is made with a pair of 
dice. What is the probability that the sum of the spots 
turned up will be 5 ? 

Number of ways of throwing the dice is 6 X 6 = 36 
Number of ways of throwing five is 4 

.,. ^ ^ • ^. • 4 1 

.•. Probability of throwing five is — = - 

Example B. A coin is tossed up six times. Find the 
chance that three heads and three tails will be the result. 
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Number of ways of throwing the coin is 2' = 64 

6x5x4 

Number of ways of throwing three heads is = 20 

•' ^ 1X2X3 

20 5 

Probability of throwing three heads is — = — 

•^ ^ 64 16 

1202. Compound Events. A certain eyent can happen in 
a ways and fail in b ways : a second independent event can 
happen in a' ways, and fail in b' ways, all of these ways 
being equally likely to occur. 

To find the probability of the simultaneous occurrence of 
the two events. 

The total number of ways in which the events can take 
place together is (a -\- b) (a' -{- b') 

( 1 ) Both events can happen in a a' ways. 

(2) Both events can fail in b b' ways. 

(3) First event can happen and second fail in a b' ways. 

(4) First event can fail and second happen in a' b ways. 



The probability of (1) is 

The probability of (2) is 

The probability of (3) is 

The probability of (4) is 



a a 




(a + *) (a' 


+ i') 


bb' 




(« + b) (a' 


+ *') 


ab' 




(a + b) (a' 


+ *') 


a' b 





(« _|_ b) {a' + b') 



But the probability of the happening of the first event is 
and of the second event is — ; — | — ^„ etc. Hence it 



a -f 6' " """ """ a' + b 

will at once be seen that the probability of the simultaneous 
occurrence of two independent events is equal to the product 
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of the probabilities of the occurrence of the component events. 

Or, in general, if Pj, Pj, . . . P„ are the probabilities of 

the occurrence of any number, n, of independent events, the 

probability of the simultaneous occurrence of all the events is 

Pi X A X ... ^„ (A) 

By independent events is meant those such that the manner 
of occurrence of one has no influence upon the manner of 
occurrence of the others. 

203. Example G. The chance that A can solve a cer- 
2 
tain problem is -, and the chance that B can solve it is 

—. Find, 
12 

(a) The probability that both will solve it. 

(b) The probability that the problem will be solved. 

For (a). This is a question as to the probability of the 
concurrence of two independent events. Therefore by an 
application of (A), the probability that both will solve the 
problem is 

2 A _ A 

3 ^ 12 ~ 18 

For (b). The problem will be solved unless both fail. 

17 7 

The probability that both will fail is - x — = — 

o \2t oD 

7 29 

The probability of getting a solution is 1 =: — 

36 36 

Example 2>. A pack of cards is cut, and those taken off 
then replaced. In how many trials will it be an even wager 
that an ace will be cut ? 
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Let n be the number of trials. Then n is to be found 
from 

\52y 2 

where the first member of the equation represents the proba- 
bility that we shall not fail n times in succession. 

Solving for n, 

log 2 

n = 5 

log 52 - log 48 

= 8.7 

In nine trials then there is a little more than an even chance 
of cutting an ace. 

204. Dependent Events. If we have a number of events 
whose modes of occurrence are aepenaent one upon another, 
the probability of their concurrence will be found by the same 
method as in pai-agraph 202 ; a now denoting the number 
of ways in which after the first event has happened the 
second will follow, and b' the number of ways in which after 
the first has happened the second will not follow, etc. Accord- 
ingly, the general formula (A) of paragraph 202 applies to 
dependent events as well as to independent ones. 

Also, if an event can take place in a variety of ways, the 
total probability of its occurrence will be the sum of the 
probabilities of its occurrence in each of the different ways. 

205. Example E. Suppose two purses contain respect- 
ively five dimes and a copper, and six dimes. A coin is taken 
at random from the first purse and placed in the second, and 
then a coin is transferred from the second to the first. What 
is the probability that the copper will remain in the first purse ? 

The probability that the copper will be taken from the first 
purse and placed in the second, and then returned to the first 
purse is 

1 1 _ _L 

6 '*^ 7 ~ 42 
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and the probability that the copper will not be taken from 
the first purse at all is 

5 
6" 

Therefore the probability that the copper will finally remain 
in the first purse is 

J_ ^ _ 36 _ 6^ 

42 ~^ "e" ~ 42 ~" y 



FUNCTIONS OF SEVERAL VARIABLES. 

206. For the application of Taylor's Theorem to the 
expansion of a function of several independent variables, see 
Osborne's " Differential and Integral Calculus," page 145. 
And for the conditions that lead to maxima and minima 
values of such functions, see page 155 of the same work. 
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TABLE I. 

2 /•< t a 

Values of the Integral — ^ I e-i'dt for Areument or - 

y/^ Jo 0.4769 r 



a 


12 3 4 


5 6 7 8 9 


Diff. 


0.0 


0.0000 0.0054 0.0108 0.0161 0.0215 


0.0269 0.0323 0.0377 0.0430 0.0484 


54 


0.1 


0638 Oo'Jl 0G45 0699 0752 


0806 0859 0913 0966 1020 


54 


0.2 


1073 1126 1180 1233 1286 


1339 1392 1445 1498 1551 


53 


0.3 


1603 1056 1709 1761 1814 


1866 1918 1971 2023 2075 


52 


0.4 


2127 2179 2230 2282 2334 


2385 2436 2488 2539 2690 


51 


0.5 


0.2641 0.2691 0.2742 0.2793 0.2843 


0.2893 0.2944 0.2994 0.3043 0.3093 


50 


0.6 


3143 3192 3242 3291 3340 


3389 3438 3487 3535 3583 


49 


0.7 


3632 3680 3728 3775 3823 


3870 3918 3965 4012 4059 


46 


0.8 


4105 4152 4198 4244 4290 


4336 4381 4427 4472 4517 


45 


0.9 


4562 4606 4651 4695 4739 


4783 4827 4860 4914 4967 


43 


1.0 


0.5000 0.5043 0.5085 0.5128 0.5170 


0.5212 0.5254 0.5295 0.5337 0.5378 


41 


1.1 


5419 5460 5500 5540 5581 


5620 5660 5700 5739 5778 


39 


1.2 


5817 5856 6894 5932 5970 


6008 6046 6083 6120 6157 


37 


1.3 


6194 6231 6267 6303 6339 


6375 6410 6445 6480 6515 


35 


1.4 


6560 6584 6618 6652 6686 


6719 6753 6786 6818 6851 


32 


1.5 


0.6883 0.6915 0.6947 0.6979 0.7011 


0.7042 0.7073 0.7104 0.7134 0.7165 


30 


1.6 


7195 7225 7255 7284 7313 


7342 7371 7400 7428 7457 


28 


1.7 


7485 7512 7540 7567 7594 


7621 7648 7675 7701 7727 


26 


1.8 


7753 7778 7804 7829 7854 


7879 7904 7928 7952 7976 


24 


1.9 


8000 8023 8047 8070 8093 


8116 8138 8161 8183 8205 


22 


2.0 


0.8227 0.8248 0.8270 0.8291 0.8312 


0.8332 0.8353 0.8373 0.8394 0.8414 


19 


2.1 


8433 8453 8473 8492 8511 


8530 8549 8567 8585 8604 


18 


2.2 


8622 8639 8657 8674 86D2 


8709 8726 8742 8759 8775 


17 


2.3 


8792 8808 8824 8840 8855 


8870 8886 8901 8916 8930 


15 


2.4 


8945 8960 8974 8988 9002 


9016 9029 9043 9056 9069 


13 


2.5 


0.9082 0.9095 0.9108 0.9121 0.9133 


0.9146 0.9158 0.9170 0.9182 0.9193 


12 


2.6 


8205 9217 9228 9239 9250 


9261 9272 9283 9293 9304 


10 


2.7 


9314 9324 9334 9344 9354 


9364 9373 9383 9392 9401 


9 


2.8 


9410 9419 9428 9437 9446 


9454 9463 9471 9479 9487 


8 


2.9 


9495 9503 9511 9519 9526 


9534 9541 9548 9556 9563 


T 


3.0 


0.9570 0.9577 0.9583 0.9590 0.9597 


0.9603 0.9610 0.9616 0.9622 0.9629 


6 


3.1 


9635 9641 9647 9652 9658 


9664 9669 9675 %>0 9686 


6 


3.2 


9691 9696 9701 9706 9711 


9716 9721 9726 9731 9736 


6 


3.3 


9740 9744 9749 9753 9757 


9761 9766 9770 9774 9778 


4 


3.4 


9782 9786 9789 9793 9797 


9800 9804 9807 9811 9814 


4 


3. 


0.9570 0.9635 0.9691 0.9740 0.9782 


0.9818 0.9848 0.9874 0.9896 0.9915 




4. 


9930 9943 9954 9963 9970 


9976 9981 9985 9988 9990 




6. 


9993 9994 9996 9997 9997 


9998 9998 9999 9999 9999 




00 


1.0000 






a 
r 


12 3 4 


5 6 7 8 9 


Diff. 
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TABLE II. — Common LogarithmB. 



n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


DM. 


10 


0000 


0043 


0086 


0128 


0170 


0212 


0253 


0294 


0334 


0374 


42 


11 


0414 


0463 


0492 


0631 


0669 


0607 


0646 


0682 


0719 


0765 


38 


12 


0792 


0828 


0864 


0899 


0934 


0969 


1004 


1038 


1072 


1106 


36 


13 


1139 


1173 


1206 


1239 


1271 


1303 


1336 


1367 


1399 


1430 


32 


14 


1461 


1492 


1523 


1553 


1584 


1614 


1644 


1673 


1703 


1732 


30 


15 


1761 


1790 


1818 


1847 


1875 


1903 


1931 


1959 


1987 


2014 


28 


16 


2041 


2068 


2095 


2122 


2148 


2175 


2201 


2227 


2263 


2279 


27 


17 


2304 


2330 


2355 


2380 


2405 


2430 


2455 


2480 


2604 


2529 


25 


18 


2553 


2677 


2601 


2625 


2648 


2672 


2695 


2718 


2742 


2766 


24 


19 


2788 


2810 


2833 


2856 


2878 


2900 


2923 


2945 


2967 


2989 


22 


20 


3010 


3032 


3054 


3075 


3096 


3118 


3139 


3160 


3181 


3201 


21 


21 


3222 


3243 


3263 


3284 


3304 


3324 


3345 


3305 


3385 


3404 


20 


22 


3424 


3444 


3464 


3483 


3602 


3522 


3541 


3560 


3579 


3598 


19 


23 


3617 


3C36 


3655 


3674 


3(>g2 


3711 


3729 


3747 


3766 


3784 


18 


24 


3002 


3820 


3838 


3856 


3874 


3892 


3909 


3927 


3945 


3962 


18 


25 


3979 


3997 


4014 


4031 


4048 


4065 


4082 


4099 


4116 


4133 


17 


26 


4100 


4166 


4183 


4200 


4216 


4232 


4249 


4265 


4281 


4298 


17 


27 


4314 


4330 


4346 


4362 


4378 


4393 


4409 


4425 


4440 


4466 


16 


28 


4472 


4487 


4502 


4518 


4633 


4648 


4564 


4579 


4694 


4609 


15 


29 


4624 


4639 


4654 


4669 


4683 


4698 


4713 


4728 


4742 


4757 


15 


30 


4771 


4786 


4800 


4814 


4829 


4843 


4857 


4871 


4886 


4900 


14 


31 


4'J14 


4928 


4942 


4965 


4969 


4983 


4997 


5011 


6024 


6038 


14 


32 


6051 


5065 


5079 


5092 


6105 


5119 


6132 


6145 


6159 


6172 


13 


33 


5185 


6198 


6211 


6224 


6237 


6250 


5263 


5276 


6289 


6302 


13 


34 


6315 


5328 


5340 


63,53 


63C6 


6378 


5391 


6403 


6416 


6428 


13 


35 


5441 


6453 


6465 


6478 


5490 


6502 


6514 


5527 


6539 


6561 


12 


36 


6563 


6575 


5587 


6699 


5011 


6623 


6635 


6647 


5658 


6670 


12 


37 


6682 


6694 


5705 


6717 


6729 


5740 


6752 


5763 


6775 


5786 


12 


38 


5798 


6809 


6821 


5832 


6843 


5865 


6866 


5877 


6888 


6899 


U 


39 


5911 


5922 


5933 


6944 


69D5 


5966 


5977 


6988 


6999 


6010 


11 


40 


6021 


6031 


6042 


6053 


6064 


6075 


6086 


6096 


6107 


6117 


11 


41 


6128 


6138 


6149 


6160 


6170 


6180 


6191 


6201 


6212 


6222 


11 


42 


6232 


6243 


6253 


6263 


6274 


6284 


6294 


6304 


6314 


6326 


10 


43 


6335 


6345 


6355 


6365 


6375 


6386 


6395 


6405 


6415 


6425 


10 


44 


6435 


6444 


6454 


6464 


6474 


6484 


6493 


6503 


6513 


6622 


10 


45 


6532 


6642 


6551 


6561 


6571 


6580 


6590 


6699 


6609 


6618 


10 


46 


6628 


6637 


6646 


6656 


6665 


6675 


6684 


6693 


6702 


6712 


9 


47 


6721 


6730 


6739 


6749 


6758 


6767 


6776 


6786 


6794 


6803 


9 


48 


6812 


6821 


6830 


6839 


6848 


6857 


6866 


6875 


6884 


6893 


9 


49 


6902 


6911 


6920 


6928 


6937 


6946 


6955 


6964 


6972 


6981 


9 


50 


6990 


6998 


7007 


7016 


7024 


7033 


7042 


7050 


7059 


7067 


9 


51 


7076 


7084 


7093 


7101 


7110 


7118 


7126 


7135 


7143 


7152 


8 


62 


7100 


7108 


7177 


7186 


7193 


7202 


7210 


7218 


7226 


7235 


8 


63 


7243 


7261 


7259 


7267 


7275 


7284 


7292 


7300 


7308 


7316 


8 


64 


7324 


7332 


7340 


7348 


7366 


7364 


7372 


7380 


7388 


1396 


8 


n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


DifE. 
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65 
56 
B7 
68 
69 

60 
61 
62 
63 
64 

65 
66 
67 
68 
69 

70 
71 
72 
73 

74 

76 
76 
77 
78 
79 

80 
81 
82 
83 
84 

85 
86 

87 
88 
89 

90 
91 
92 
93 
94 

96 
96 
97 
98 



7404 
7482 
7559 
7634 
7709 

7782 
7863 
7924 
7993 
8062 

8129 
8195 
8261 
8325 



8451 
8513 
8573 
8G33 
8692 

8751 
8808 
8865 
8921 
8976 

9031 
9085 
9138 
9191 



9294 
9345 
9395 
9445 
9494 

9542 
9590 
9638 
9685 
9731 

9777 
9823 
9868 
9912 
9956 



7412 
7490 
7666 
7642 
7716 

7789 
7860 
7931 
8000 
8069 

8136 
8202 
8267 
8331 
8395 

8457 
8519 
8679 
8639 
8698 

8756 
8814 
8871 
8927 



9036 
9090 
9143 
9196 
9248 

9299 

9350 
»i00 
9450 
9499 

9647 
9595 
9643 
9689 
9736 

9782 
9827 
9872 
9917 
9961 



7419 
7497 
7674 
7649 
7723 

7796 
7868 
7938 
8007 
8075 

8142 
8209 
8274 
8338 
8401 

8463 
8625 
8586 
8645 
8704 

8762 
8820 
8876 
8932 
8987 

9042 
9096 
9149 
9201 
9253 

9304 
9356 
9406 
9456 
9504 

9552 
9600 
9647 
9694 
9741 

9786 
9832 
9877 
9921 
9965 



7427 
7606 
7682 
7657 
7731 

7803 

7875 
7946 
8014 
8082 

8149 
8216 
8280 
8344 
8407 

8470 
8631 
8691 
8651 
8710 

8768 
8825 
8882 
8938 



9047 
9101 
9154 
9206 
9268 

9309 
9360 
9410 
9460 
9609 

9567 
9605 
9652 
9699 
9746 

9791 
9836 
9881 
9926 
9969 



7436 
7513 
7589 
7664 
7738 

7810 
7882 
7962 
8021 
8089 

8156 
8222 
8287 
8351 
8414 

8476 

8637 
8597 
8657 
8716 

8774 
8831 
8887 
8943 
8998 

9053 
9106 
9169 
9212 
9263 

9315 
9365 
9415 
9465 
9613 

9662 
9b09 
9657 
9703 
9750 

9796 
9841 
9886 
9930 
9974 



5 


6 


7 


8 


9 


7443 


7461 


7459 


7466 


7474 


7620 


7528 


7536 


7643 


7661 


7597 


7604 


7612 


7619 


7627 


7672 


7679 


7686 


7694 


7701 


7746 


7752 


7760 


7767 


7774 


7818 


7826 


7832 


7839 


7846 


7889 


7896 


7903 


7910 


7917 


7959 


7966 


7973 


7980 


7987 


8028 


8035 


8041 


8048 


8055 


8096 


8102 


8109 


8116 


8122 


8162 


8169 


8176 


8182 


8189 


8228 


8235 


8241 


8248 


8264 


8293 


8299 


8306 


8312 


8319 


8357 


8363 


8370 


8376 


8382 


8420 


8426 


8432 


8439 


8446 


8482 


8488 


8494 


8600 


8606 


8543 


8549 


8556 


8661 


8567 


8603 


8609 


8615 


8621 


8627 


8663 


8669 


8676 


8681 


8GS6 


8722 


8727 


8733 


8739 


8745 


8779 


8785 


8791 


8797 


8802 


8837 


8842 


8848 


8854 


8859 


8893 


8899 


8904 


8910 


8915 


8949 


8954 


8960 


8965 


8971 


9004 


9009 


9016 


9020 


9025 


9058 


9063 


9069 


9074 


9079 


9112 


9117 


9122 


9128 


9133 


9165 


9170 


9175 


9180 


9186 


9217 


9222 


9227 


9232 


9238 


9269 


9274 


9279 


9284 


9289 


9320 


9325 


9330 


9335 


9340 



9370 
9420 
9469 
9518 

9566 
9614 
9661 
9708 
9754 

9800 
9845 
9890 
9934 
9978 



9375 9380 

9425 9430 

9474 9479 

9523 9528 

9571 9576 

9619 9624 

9666 9671 

9713 9717 

9769 9763 

9805 9809 

9850 9854 

9894 9899 

9939 9943 

9983 9987 



9385 
9435 
9484 
9533 

9581 
9628 
9675 
9722 
9768 

9814 
9859 
9903 
9948 
9991 



9390 
9440 
9489 
9638 

9586 
9633 
9680 
9727 
9773 

9818 
9863 
9908 
9952 
9996 



Dift. 
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TABLE III. — Squares of Numbers. 





T* 





1 


2 


3 


4 


5 


6 


7 


8 


9 


DifE. 




1.0 


1.000 


1.020 


1.040 


1.061 


1.082 


1.103 


1.124 


1.145 


1.166 


1.188 


22 




1.1 


1.210 


1.232 


1.254 


1.277 


1.300 


1.323 


1.346 


1.369 


1.392 


1.416 


24 




1.2 


1.440 


1.464 


1.488 


1.513 


1.538 


1.563 


1.588 


1.613 


1.638 


1.664 


26 




1.3 


1.690 


1.716 


1.742 


1.769 


1.796 


1.823 


1.850 


1.877 


1.904 


1.932 


28 




1.4 


1.960 


1.988 


2.016 


2.016 


2.074 


2.103 


2.132 


2.161 


2.190 


2.220 


30 




1.5 


2.250 


2.280 


2.310 


2.341 


2.372 


2.403 


2.434 


2.465 


2.496 


2.528 


32 




1.6 


2.560 


2.692 


2.624 


2.657 


2.690 


:i.723 


2.756 


2.789 


2.822 


2.866 


34 




1.7 


2.890 


2.924 


2.968 


2.993 


3.028 


3.UU3 


3.098 


3.133 


3.168 


3.204 


36 




1.8 


3.240 


3.276 


3.312 


3.349 


3.386 


3.423 


3.460 


3.497 


3.534 


3.572 


38 




1.9 


3.610 


3.648 


3.686 


3.726 


3.764 


3.803 


3.842 


3.881 


3.920 


3.960 


40 




2.0 


4.000 


4.040 


4.080 


4.121 


4.162 


4.203 


4.244 


4.286 


4.326 


4.368 


42 




2.1 


4.410 


4.452 


4.494 


4.537 


4.580 


4.623 


4.666 


4.709 


4.752 


4.796 


44 




2.2 


4.840 


4.884 


4.928 


4.973 


5.018 


6.063 


5.108 


6.153 


5.198 


5.244 


46 




2.3 


6.290 


6.336 


5.382 


6.429 


6.476 


6.023 


5.570 


5.617 


5.664 


5.712 


48 




2.4 


6.760 


6.808 


5.866 


5.905 


6.964 


6.003 


6.052 


6.101 


6.150 


6.200 


60 




2.5 


6.250 


6.300 


6.350 


6.401 


6.452 


6.503 


6.554 


6.605 


6.656 


6.708 


52 




2.6 


6.760 


6.812 


6.864 


6.917 


6.970 


7.023 


7.076 


7.129 


7.1S2 


7.236 


64 




2.7 


7.290 


7.344 


7.398 


7.453 


7.608 


7.563 


7.618 


7.673 


7.728 


7.784 


66 




2.8 


7.840 


7.8% 


7.952 


8.009 


8.066 


8.123 


8.180 


8.237 


8.294 


8.352 


58 




2.9 


8.410 


8.468 


8.526 


8.686 


8.644 


8.703 


8.762 


8.821 


8.860 


8.940 


60 




3.0 


9.000 


9.060 


9.120 


9.181 


9.242 


9.303 


9.364 


9.425 


9.486 


9.548 


62 




3.1 


9.610 


9.672 


9.734 


9.797 


9.860 


9.923 


9.986 


10.05 


10.11 


10.18 


6 




3 2 


10.24 


10.50 


10.37 


10.43 


10.50 


10.56 


10.63 


10.69 


10.76 


10.82 


7 




3.3 


10.89 


10.96 


11.02 


11.09 


11.16 


11.22 


11.29 


11.36 


11.42 


11.49 


7 




3.4 


11.66 


11.63 


11.70 


11.76 


11.83 


11.90 


11.97 


12.04 


12.11 


12.18 


7 




3.5 


12.25 


12.32 


12.39 


12.46 


12.53 


12.60 


12.67 


12.74 


12.82 


12.89 


7 




3.6 


12.96 


13.03 


13.10 


13.18 


13.25 


13.32 


13.40 


13.47 


13.54 


13.62 


7 




3.T 


13.69 


13.76 


13.84 


13 91 


13.99 


14.06 


14.14 


14.21 


14.29 


14.36 


8 




3.8 


14.44 


14.52 


14.59 


14.67 


14.75 


14.82 


14.90 


14.98 


15.05 


15.13 


8 




3.9 


16.21 


15.29 


15.37 


15.44 


15.52 


15.60 


15.68 


15.76 


15.84 


15.92 


8 




4.0 


16.00 


16.08 


16.16 


16.24 


16.32 


16.40 


16.48 


16.56 


16.65 


16.73 


8 




4.1 


16,81 


16.89 


16.97 


17.06 


17.14 


17.22 


17.31 


17.39 


17.47 


17.56 


8 




4.2 


17.64 


17.72 


17.81 


17.89 


17.98 


18.06 


18.15 


18.23 


18.32 


18.40 


9 




4.3 


18.49 


18.68 


18.66 


18.76 


18.84 


18.92 


19.01 


19.10 


19.18 


19.27 


9 




4.4 


19.36 


19.45 


19.54 


19.62 


19.71 


19.80 


19.89 


19.98 


20.07 


20.16 


9 




4.5 


20.25 


20.34 


20.43 


20.52 


20.61 


20.70 


20.79 


20.88 


20.98 


21.07 


9 




4.6 


21.16 


21.25 


21.34 


21.44 


21.53 


21.62 


21.72 


21.81 


21.90 


22.00 


9 




4.7 


22.09 


22.18 


22.28 


22.37 


22.47 


22.56 


22.66 


22.75 


22.85 


22.94 


10 




4.8 


23.04 


23.14 


23.23 


23.33 


23.43 


23.52 


23.62 


23.72 


23.81 


23.91 


10 




4.9 


24.01 


24.11 


24.21 


24.30 


24.40 


24.50 


24.60 


24.70 


24.80 


24.90 


10 




5.0 


25.00 


25.10 


25.20 


25.30 


25.40 


25.50 


25.60 


25.70 


26.81 


25.91 


10 




6.1 


26.01 


26.11 


26.21 


26.32 


26.42 


26.52 


26.63 


26.73 


26.83 


26.94 


10 




5.2 


27.04 


27.14 


27.25 


27.35 


27.46 


27.66 


27.67 


27.77 


27.88 


27.98 


11 




6.3 


28.09 


28.20 


28.30 


28.41 


28.52 


23.62 


28.73 


28.84 


28.94 


29.05 


11 




6.4 


29.16 


29.27 


29.38 


29.48 


29.59 


29.70 


29.81 


29.92 


30.03 


30.14 


11 




H 





1 


2 


3 


4 


5 


6 


7 


8 


9 
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TABLE III. — Squares of Numbers. 



n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


Diff. 


6.5 


30.26 


30.36 


30.47 


30.58 


30.69 


30.80 


30.91 


31.02 


31.14 


31.25 


11 


6.6 


31.36 


31.47 


31.58 


31.70 


31.81 


31.92 


32.04 


32.15 


32.26 


32.38 


11 


8.7 


32.49 


32.60 


32.72 


32.83 


32.95 


33.06 


33.18 


33.29 


33.41 


33.52 


12 


6.8 


33.64 


33.76 


33.87 


33.99 


34.11 


34.22 


34.34 


34.46 


34.57 


34.69 


12 


6.9 


34.81 


34.93 


36.06 


35.16 


35.28 


35.40 


35.52 


.'i5.64 


35.76 


35.88 


12 


6.0 


36.00 


36.12 


36.24 


36.36 


36.48 


36.60 


36.72 


36.84 


36.97 


37.09 


12 


6.1 


37.21 


37.33 


37.45 


37.68 


37.70 


37.82 


37.95 


38.07 


38.19 


38.32 


12 


6.2 


38.44 


38.56 


38.69 


38.81 


38.94 


39.06 


39.19 


39.31 


39.44 


39.56 


13 


6.3 


39.69 


39.82 


39.94 


40.07 


40.20 


40.32 


40.46 


40.58 


40.70 


40.83 


13 


6.4 


40.96 


41.09 


41.22 


41.34 


41.47 


41.60 


41.73 


41.86 


41.99 


42.12 


13 


6.5 


42.25 


42.38 


42.51 


42.64 


42.77 


42.90 


43.03 


43.16 


43.30 


43.43 


13 


6.6 


43.56 


43.69 


43.82 


43.96 


44.09 


44.22 


44.36 


44,49 


44.62 


44.76 


13 


6.7 


44.89 


45.02 


45.16 


45.29 


45.43 


45.66 


45.70 


45.83 


45.97 


46.10 


14 


6.8 


46.24 


46.38 


46.51 


46.66 


46.79 


46.92 


47.06 


47.20 


47.33 


47.47 


14 


6.9 


47.61 


47.75 


47.89 


48.02 


48.16 


48.30 


48.44 


48.68 


48.72 


48.86 


14 


7.0 


49.00 


49.14 


49.28 


49.42 


49.56 


49.70 


49.84 


49.98 


60.13 


50.27 


14 


7.1 


60.41 


60.55 


50.69 


60.84 


50.98 


51.12 


51.27 


51.41 


51.55 


51.70 


14 


7.2 


61.84 


61.98 


62.13 


52.27 


52.42 


52.56 


52.71 


52.85 


53.00 


53.14 


15 


7.3 


63.29 


53.44 


53.68 


63.73 


53.88 


54.02 


64.17 


54.32 


54.46 


54.61 


15 


7.4 


64.76 


54.91 


55.06 


66.20 


65.35 


55.50 


65.65 


65.80 


65.95 


56.10 


15 


7.5 


56.25 


66.40 


66.55 


66.70 


66.85 


57.00 


57.15 


57.30 


57.46 


57.61 


15 


7.6 


57.76 


57.91 


58.06 


68.22 


58.37 


58.52 


58.68 


68.83 


68.98 


59.14 


15 


7.7 


69.29 


59.44 


59.60 


69.75 


59.91 


60.06 


60.22 


60.37 


60.53 


60.68 


16 


7.8 


60.84 


61.00 


61.15 


61.31 


61.47 


61.62 


61.78 


61.94 


62.09 


62.26 


16 


7.9 


62.41 


62.57 


62.73 


62.88 


63.04 


63.20 


63.36 


63.52 


63.68 


63.84 


IS 


8.0 


64.00 


64.16 


64.32 


64.48 


64.64 


64.80 


64.96 


65.12 


65.29 


65.45 


16 


8.1 


65.61 


65.77 


65.93 


66.10 


66.26 


66.42 


66.59 


66 75 


66.91 


67.08 


16 


8.2 


67.24 


67.40 


67.57 


67.73 


67.90 


68.06 


68.23 


68.39 


68.66 


68.72 


17 


8.3 


68.89 


69.06 


69.22 


69.39 


69.56 


69.72 


69.89 


70.06 


70.22 


70.39 


17 


8.4 


70.66 


70.73 


70.90 


71.06 


71.23 


71.40 


71.67 


71.74 


71.91 


72.08 


17 


8.6 


72.25 


72.42 


72.59 


72.76 


72.93 


73.10 


73.27 


73.44 


73.62 


73.79 


17 


8.6 


73.96 


74.13 


74.30 


74.48 


74.65 


74.82 


75.00 


75.17 


75.34 


75.52 


17 


8.7 


76.69 


76.86 


76.04 


76.21 


76.39 


76.56 


76.74 


76.91 


77.09 


77.26 


18 


8.8 


77.44 


77.62 


77.79 


77.97 


78.15 


78.32 


78.50 


78.68 


78.85 


79.03 


18 


8.9 


79.21 


79.39 


79.57 


79.74 


79.92 


80.10 


80.28 


80.46 


80.64 


80.82 


18 


9.0 


81.00 


81.18 


81.36 


81.54 


81.72 


81.90 


82.08 


82.26 


82.45 


82.63 


18 


9.1 


82.81 


82.99 


83.17 


83.36 


83.54 


83.72 


83.91 


84.09 


84.27 


84.46 


18 


9.2 


84.64 


84.82 


85.01 


85.19 


85.38 


85.56 


85.75 


85.93 


86.12 


86.30 


19 


9.3 


86.49 


86.68 


86.86 


87.05 


87.24 


87.42 


87.61 


87.80 


87.98 


88.17 


19 


9.4 


88.36 


88.55 


88.74 


88.92 


89.11 


89.30 


89.49 


89.68 


89.87 


90.06 


19 


9.6 


90.25 


90.44 


90.63 


90.82 


91.01 


91.20 


91.39 


91.58 


91.78 


91.97 


19 


9.6 


92.16 


92.35 


92.54 


92.74 


92.93 


93.12 


93.32 


93.61 


93.70 


93.90 


19 


9.7 


94.09 


94.28 


94.48 


94.67 


94. 8T 


95.06 


95.26 


95.45 


95.65 


95.84 


20 


9.8 


96.04 


96.24 


96.43 


96.63 


96.83 


97.02 


97.22 


97.42 


97.61 


97.81 


20 


9.9 


98.01 


98.21 


98.41 


98.60 


98.80 


99.00 


99.20 


99.40 


99.60 


99.80 


20 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 
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